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ABSTRACT 


We  present  the  results  of  a  study  on  the  current- 
induced  breakdown  of  superconduct i vity  in  a  quasi-one- 
dimensional  filament.  We  interpret  this  breakdown  in 
terms  of  the  ohase-slip  centers  (PSC's)  introduced  by 
Skocpol,  Beasley,  and  Tinkham  (SBT) ,  and  emphasize 
throughout  the  central  role  played  by  the  nonequilibrium 
quasiparticles.  Our  study  consists  of  two  parts,  one 
largely  experimental,  the  other  largely  theoretical. 

In  our  experimental  work,  we  have  investigated 
the  effects  of  a  magnetic  field  H  on  nonequilibrium 
quasiparticle  relaxation  processes  in  superconductors 
by  analyzing  the  I-V  characteristics  of  long  tin  micro¬ 
bridges  in  a  parallel  field  near  T^ (H) .  Following  SBT, 
we  take  the  differential  resistance  of  a  given  PSC  to 
be  approximately  equal  to  the  normal  resistance  of 
a  nonequilibrium  region  of  the  bridge  with  length  2/Dt  , 

K 

where  D  is  the  quasiparticle  diffusion  constant  and  t 

R 

is  the  transverse  mode  (branch  imbalance)  relaxation 

time.  The  magnitude  and  the  temperature  and  field 

dependence  of  r  as  inferred  from  our  data  aqree  well 
K 

with  those  of  the  transverse  mode  relaxation  time  in 
the  presence  of  pair-breaking,  derived  by  Schmid  and 
Schon.  The  longitudinal  mode  disequilibrium  is  inter- 


i 


preted  in  terms 


1 1 

of  local  he.it  mo  prouort.  ional  to  IV, 
and  the  result  mo  ilistoi  t  ions  in  t  lit'  1 -V  character¬ 
istics  .ire  corrected  foi  model  t  o  .mounts  of  he.it  inq. 

Most  or  our  samples  oont.im  .i  vie  liberate  weak,  spot 
with  a  depressed  entie.il  current  I  which  isolates 

a  simile  PSC  and  lowers  the  heat  dissmat  ton  ('  1  K  )  . 

e  n 

The  depression  of  1  „  or  the  weak  snot  relative'  to  that 
of  the  uniform  filament  is  incompletely  understood, 
but  is  stronqly  dependent  on  the  direction  and  magnitude 
of  the  magnetic  field. 

In  our  theoretie.il  studies,  usinu  a  oener.i  1  i zed 
two-fluid  Picture  to  doseribe  a  quas  i -one-dimens  ional 
superconductor  no.n  T  t t or  HO',  we  nrovide  a  heuristic 
derivat  ion  foi  a  set  of  eou.it  ions  uoverninq  t  tie  temporal 
and  spatial  evolution  of  the  choree  imbalance  tot 
branch  imbalance)  in  the  .mas i part i cles .  We  shew  that 
these  equ.it  ions  are  isomorphic  to  those  which  describe 
a  simple  electrical  transmission  line,  so  that  oh.iroe 
imbalance  waves  may  prepa.jate  m  tin'  superconductor 
m  analeoy  with  electrical  sionals  which  or 00.10 ate 
down  the  transmission  line.  We  propose  as  .1  model  tor 
a  PSC  in  a  supet  conduct  1  no  filament  .1  localized  Joseph- 
son  c'sc.il  later  coupled  to  the  transmission  lint'.  Apply¬ 
ing  standard  t ransrni ss ion- l 1 ne  theory  to  solve  the 
problem,  we  show  that  the  dosephson  oscillations  in  the 
center  viener.it  e  charuc  mb.i  lanuu  waves  which  prepa.jate 
t'ut  to  .1  :  1  oquonoy-dopond.ent  distance  of  ordet  ot  the 
quasipart  tele  diffusion  lenoth  )  ^  >  Pi  „  before  they 


■ 


irtwM****** 


damp  out.  The  time-averaged  behavior  of  the  model 
reduces  to  the  earlier  SBT  model.  A  novel  consequence 
of  the  model  is  a  prediction  of  intrinsic  hysteresis 
in  the  dc  current-voltaqe  relation.  The  model  also 
provides  a  convenient  framework  for  dealing  with  ac 
effects  in  PSC's,  including  resonance  and  synchronization 
in  systems  of  closely-spaced  PSC’s  and  microbridges. 
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IVKOPIVTIOM  AMP  I'RI'I. IVINARY  PISPPSSION 


The  phi’noim’tu'n  of  supet  conduct  1  vi  t  y  was  named  lot 
the  tom.it  kablo  ability  of  a  mot.il  in  this  st  .it  o  to 
carry  .1  sto.i.ly  out  tent  with  no  d i sss i pat i on ,  t.o.  w  \  t  h 
no  volt.iuo  t » 1  .It  iv*’  t  ho  current  .  It  w.is  reoo.inwed 
out  ly  on,  howovot  ,  that  tins  ability  will  break  down 
fot  some  sufficiently  lat.io  out  1  out  .  A  .liven  sample 
need  not  necessarily  >10  .liteotly  into  the’  olassto.illy 
norma  l -met  a  1 1  i  o  state  upon  exceedin']  the'  otitto.il  om  1  rnt 
howevot  .  The  case  of  vortex  t  low  in  .1  type  11  supet 
oon.luotoi  is  one  well-known  example.  The’  t  1  ur  hallm.uk 
of  supet  oon.luot  1  v  1 1  y  1  s  actually  the  existence  of  an 
older  pat ameter  with  lonu-t  anue  phase  ooliet  enee  ,  tat  hot 
than  the  infinite  conductivity  which  cave  it  its  name. 

A  supet  conduct  ot  can  .10  into  .1  none.pn  1  i  ht  1  um  state 
which  tet  ai  ns  this  ohaiaetet  while  also  exhibit  mu 
d i ss i pat  i ve  e f  f  ect  s . 

It  ts  often  the  case  that  .1  simple  o comet ry 
.treat  ly  facilitates  the  study  of  Physical  phenomena . 

A  supet conduct  in. 1  filament  is  one  such  system,  which 
is  oeornet  r  i  0.1 1  ly  simple,  but  makes  foi  .1  sm  pi  i  s  i  no  1  y 
rich  and  complex  problem.  If  the  cross  sectional 
dimensions  are  '-mill  compared  with  the  din.-buto  landau 
coherence  lenuth  (which  in  pi  act  ice  may  be’  about  a 
micron)  ,  the’  filament  can  be  treated  as  effectively 
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one-dimensional.  A  real  version  of  such  an  idealized 
sample  can  be  a  thin-film  microstrip  or  whisker  crystal, 
with  length  of  order  100  microns.  The  current-induced 
breakdown  of  superconductivity  in  such  a  sample  is  still 
not  understood  in  its  full  complexity.  We  will  deal 
with  the  problem  in  terms  of  approximations  which  are 
amenable  to  ohysical  interpretation. 

The  resistive  transition  of  such  samples  near  the 

critical  temperature  T  ,  as  seen  in  the  dc  I-V  curves, 

c 

is  characterized  by  the  appearance  of  a  series  of  voltage 
steps,  each  followed  by  a  plateau  of  approximately  con¬ 
stant  differential  resistance^  ^ ,  as  in  Fig.  1.1.  Each 
increment  of  differential  resistance  is  about  the  same 
from  one  plateau  to  the  next,  at  least  for  the  first 

few  steps.  This  basic  picture  has  been  observed  in 

1  2  3 

microstrips  of  Sn  (as  in  Fig.  1.1) ,  Al  ,  In  ,  and 

4 

other  metals  ,  and  in  whisker  crystals  of  Sn  and  Sn 
alloy^.  This  constant  increment  of  dif f ferential  resist¬ 
ance  suggests  that  what  is  happening  is  the  successive 
onset  of  basically  similar  dissipative  units  as  the 
current  is  increased,  and  that  the  total  voltage  drop 
along  the  filament  is  simply  the  sum  of  the  voltage 
drops  of  the  individual  units.  This  was  shown  directly 
by  Skocpol^,  who  fabricated  a  microstrip  with  a  series 
of  side-tabs  located  at  intervals  along  the  length  of 
the  strip,  which  acted  as  voltage  probes.  The  dissipative 
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I  ( MA) 

Fig  1.1  Current-voltage  curve  for  a  Sn  microstrip. 

Vote  the  steo  structure  with  the  incremental 
increase  in  differential  resistance.  The  inset 
shows  the  corresponding  sample  (tin  bridge  #25A 
in  Table  2.1). 
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units  wore  found  to  enter  singly,  and  tended  to  space 
themselves  relatively  far  apart.  It  makes  sense, 
therefore,  to  speak  in  terms  of  a  series  of  discrete 
spat i al ly- loca 1 i zed  un its. 

Skocpol ,  Beasley,  and  Tinkham  proposed  a  model  in 
which  each  of  these  units  is  identified  as  a  "phase-slip 
center"  (PSC)  F.ach  PSC  first  appears  at  its  own  local 
value  of  the  critical  current,  as  determined  by  the  local 
critical  temperature  and  cross-sectional  area,  and  by 
the  presence  of  other  PSO’s  nearby.  The  differential 
resistance  of  each  PSC  is  approximately  equal  to  the 
normal-state  resistance  of  a  length  2A,  where  A  is  the 
quasiparticle  diffusion  length  associated  with  a  relax¬ 
ation  time  t  which  is  of  the  order  of  magnitude  charac- 
R 

teristic  of  inelastic  electron-phonon  processes  in  the 

metal.  In  the  Sn  sample  of  Fig.  1.1,  A  is  about  6  urn 
-  n 

and  i  -  10  sec. 

K 

The  nonequilibrium  region  is  not  fully  normal, 
however,  since  extrapolating  the  straight-line  slope 
down  to  the  current  axis  gives  a  positive  value.  This 
implies  that  only  a  portion  of  the  current,  associated 
with  the  excited  quasiparticles,  contributes  resistance, 
and  the  remainder,  a  constant  fraction  ot  the  critical 
current,  is  carried  as  dissipationless  supercurrent. 

In  effect,  we  are  dealing  here  with  a  generalized  two- 
fluid  picture. 
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Those  two  "fluids"  respond  to  two  different  electro¬ 
chemical  potentials,  un  for  the  effectively  "normal" 
quasiparticles  and  u  for  the  superconduct ing  pairs. 

These  two  are  equal  in  equilibrium  far  from  the  center 
of  the  dissipative  region,  but  show  fundamentally 
different  behavior  in  the  nonequilibrium  region  near  a 
PSC.  As  indicated  in  Fig.  1.2a,  a  voltage  drop  across 
the  PSC  must  appear  as  a  chanqe  in  both  of  these  poten¬ 
tials,  but  the  change  occurs  over  two  different  length 
scales.  A  qradient  in  ug  produces  an  acceleration  of 
the  supercurrent  I  and  leads  to  a  relaxation  oscillation 
(or  phase-slip  cycle)  of  the  order  parameter  within  a 
region  about  the  size  of  the  Ginzburg-Landau  coherence 
length  (T)  .  The  time-averaged  potential  u  changes 
abruptly  within  this  region  and  is  constant  outside. 

The  phase-slip  process  causes  the  time-averaged  super¬ 
current  I  at  the  core  region  to  be  depressed  below  the 
value  far  away  (see  Fig.  1.2b) ,  so  that  if  the  total 
current  is  to  remain  constant,  the  dissipative  normal 
current  I  ,  which  is  driven  by  gradients  in  the  "normal" 
electrochemical  potential  un ,  must  make  up  the  difference. 
This  interconversion  between  normal  current  and  super¬ 
current  must  occur  over  a  scale  of  the  quasiparticle 
diffusion  length  A  (which  is  typically  larger  than  £.  ), 

the  length  required  for  the  nonequilibrium  quasiparticles 
carrying  the  normal  current  to  come  back  into  equilibrium 
with  the  pair  condensate. 


7 

In  the  limit  where  A>>£,  which  is  what  we  are 
considerinq  in  Fig.  1.2,  the  details  of  the  superconducting 
dynamics  in  the  core  region,  which  produce  the  actual 
phase-slip  process  itself,  are  not  crucial  to  the  model, 
and  we  will  for  the  most  part  ignore  them.  However,  the 
fact  that  there  are  high-frequency  Josephson  oscillations 
in  the  core  ij>  important  in  determining  the  time-depen¬ 
dent  behavior  of  the  currents  and  potentials  in  the  larger 
region  outside  the  core,  and  in  later  chapters  a  more 
detailed  model  will  be  developed  which  brings  this 
time-dependent  behavior  into  account. 

In  order  to  place  the  SBT  picture  of  the  PSC  in 
perspective,  we  will  briefly  survey  the  field  of  non¬ 
equilibrium  superconductivity . 6  A  superconductor  may 
be  described  as  consisting  of  three  interacting  systems: 
the  superconducting  condensate,  the  quasiparticle  ex¬ 
citations,  and  the  phonons.  The  condensate,  which  can 
carry  a  dissipationless  supercurrent,  can  be  character¬ 
ized  by  a  complex  order  parameter,  whose  magnitude  in 
a  BCS  superconductor  is  proportional  to  the  gap  A  in  the 
excitation  spectrum.  Quasiparticles  with  high  energy 
compared  to  the  energy  gap  A  can  be  regarded  as  either 
electrons  or  holes,  but  low  energy  quasiparticles  have 
a  mixed  character,  although  the  excitation  spectrum  is 
customarily  divided  into  an  electron-like  and  a  hole-like 
branch  according  to  the  predominant  character  of  the 


excitations.  In  thermal  equilibrium,  brought  about  by 
inelastic  collisions  between  the  quasiparticles  and  the 
phonons,  the  energy  distribution  of  the  quasiparticles 
is  given  equally  on  both  branches  by  the  usual  Fermi 
function,  so  that  there  is  no  net  charge  in  the  quasi¬ 
particles.  An  external  perturbation  which  drives  the 
system  away  from  equilibrium  may  alter  the  distributions 
of  both  the  quasiparticles  ana  the  phonons.  As  far  as 
the  quasiparticles  are  concerned,  there  are  two  general 
types  of  perturbations,  as  discussed  below. 

An  "uncharged"  perturbation  such  as  the  absorption 
of  photons  or  phonons  may  alter  the  number  and  energy 
distribution  of  the  quasiparticles,  producing  similar 
effects  on  the  distributions  of  both  electron-like  and 
hole-like  excitations.  This  will  alter  the  magnitude  of 
the  energy  gap  A  in  the  excitation  spectrum  in  accordance 
with  the  BCS  self-consistent  gap  equation,  and  relaxation 
will  occur  by  electron-phonon  processes  involving 
scattering  and  recombination  of  the  quas ioarticles . 

Such  a  perturbation  may,  of  course,  depress  the  gap, 
producing  an  effective  rise  in  the  electronic  tempera¬ 
ture.  If  the  perturbation  acts  to  redistribute  quasi¬ 
particles  away  from  the  gap  edge,  however,  the  opposite 
effect,  namely  gap  enhancement,  can  occur.7  Furthermore, 
certain  strongly  nonequilibrium  distributions  of  quasi- 
particles  can  correspond  to  multiple  solutions  of  the 


gap  equation,  permitting  the  possibility  of  an  inhomo¬ 
geneous  state  consisting  of  regions  with  differing 
9 

energy  gaps . 

In  contrast,  a  "charged"  perturbation,  such  as  the 
tunneling  injection  of  ouasipart icles  or  the  intercon¬ 
version  of  normal  and  supercurrent,  may  involve  the 
preferential  buildup  of  quasiparticles  on  either  the 

electron-like  or  the  hole-like  branch  of  the  excitation 
9 

spectrum  .  Since  this  "branch  imbalance"  also  corresponds 
to  a  net  excess  charge  in  the  quasioarticles ,  the  term 
"charge  imbalance"  has  also  been  used.  The  pair  con¬ 
densate  will  respond  to  maintain  electrical  neutrality, 
causing  a  shift  in  the  chemical  potential  of  the  con¬ 
densate.*0  A  region  with  nonzero  charge  imbalance  in 
the  quasipart icles  will  therefore  have  a  difference 
of  electrochemical  potentials  between  the  condensate 
and  the  quas ipart icles ,  which  is  observable  using  super¬ 
conducting  and  normal  voltage  probes  in  such  a  nonequi- 
librium  region. 

These  same  two  classes  of  perturbat ions  have  been 

treated  from  an  alternative  theoretical  viewpoint  by 
1 2 

Schmid  and  Schon.  ‘  They  showed  that  any  small  deviation 
from  equilibrium  can  be  resolved  into  two  characteristic 
modes,  which  under  certain  circumstances  can  be  dealt 
with  essentially  independently.  Considering  the  complex 
superconducting  order  parameter  .is  a  vector  in  the 
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complex  plane,  the  "longitudinal"  and  "transverse" 
modes  correspond  respectively  to  changes  along  the 
direction  of  this  vector  and  perpendicular  to  it. 

The  former  is  associated  with  uncharged  perturbations 
which  change  the  magnitude  of  the  gap,  while  the  latter 
is  associated  with  nonequilibrium  values  of  the  order 
parameter  phase,  and  hence  with  nonequilibrium  potentials 
and  currents  involved  with  charged  perturbations. 

We  will  be  discussing  mostly  problems  involving 
transverse  disequil  ibria,  although  with  the  understanding 
that  in  practice  a  certain  amount  of  associated  longi¬ 
tudinal  mode  disequilibrium  (which  can  often  be  approx¬ 
imated  as  simple  heating)  is  unavoidable  because  of 
dissipative  effects.  ^he  conceptually  most  direct 
experiment  in  this  area  is  the  double  tunnel  junction 
injection  experiment  (see  Fig.  1.3a)  originally  performed 
by  Clarke,^  in  which  a  quas ioart icle  current  is  injected 
into  the  middle  (superconducting)  film  of  the  sandwich 
through  one  junction,  and  the  electrochemical  potential 
of  the  quasioarticles  in  the  middle  film  is  probed  by 
the  normal-metal  electrode  of  the  other  junction.  A 
potential  difference  can  be  observed  between  the  super¬ 
conducting  lead  to  the  middle  film  and  the  normal- 
metal  electrode.  The  injection  is  spatially  uniform 
within  the  nonegui librium  region,  and  the  steady-state 
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potential  difference  is  proportional  to  the  injection 

9 

current  and  the  charge  inbalance  relaxation  time. 

The  measured  relaxation  time  is  consistent  near  T  with 

c 

12 

the  expression 


'Q* 


4kBT  T, 


(1.1) 


where  is  the  equilibrium-averaged  electron-phonon 

(inelastic)  scattering  time  of  quasiparticles  near  the 

Fermi  surface.  The  factor  of  knT/A  arises  because 

branch  imbalance  relaxation  occurs  primarily  via  states 

between  A  and  -2A.  This  situation  has  been  analyzed 

from  a  number  of  theoretical  perspectives,  with  generally 

14 

consistent  results. 

A  related  set  of  experiments  have  investigated  the 
quasipart icle-condensate  potential  difference  in  the 
vicinity  of  a  current-carrying  normal-superconducting 
boundary.  Observations  of  the  resistance  of  the  inter¬ 
mediate  state  and  of  S/M/S  sandwiches  have  demonstrated 
that  a  measurable  excess  resistance  is  associated  with 
the  superconducting  side  of  such  a  boundary. Observa¬ 
tions  by  Yu  and  Vercereau^  with  superconducting  and 
normal  probes  placed  within  a  few  microns  of  such  a 
boundary  in  a  thin-film  configuration  (see  Fig.  1.3b) 
confirmed  the  existence  of  a  potential  difference  be¬ 
tween  the  quasiparticles  and  the  condensate.  The  results 
were  consistent  with  a  spatially  uniform  condensate 
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electrochemical  potential  n  (measured  by  the  super¬ 
conduct  inq  probe)  and  an  exponential  relaxation  of  the 
quasiparticle  electrochemical  potential  p  (measured  by 
the  normal  probe)  to  the  condensate  value  with  a  charac¬ 
teristic  lenqth  scale  much  longer  than  the  C.b  coherence 
length.  Subsequent  theoretical  analyses  1 °  '  1  7  have 
concluded  that  this  exponential  spatial  variation  of  the 
quasiparticle-condensate  potential  difference  should  be 

novcrned  by  the  quas ioart icle  diffusion  lenqth  A _  = /57TT  , 

v 

where  D=v[,i;/1  is  the  normal-state  diffusion  constant. 

The  additional  resistance  contribution  from  the  super¬ 
conducting  side  of  the  S/N  boundary  is  thus  essentially 

the  normal  resistance  of  a  length  A_*  of  the  material  on 

Q" 

that  side,  although  the  situation  is  complicated  somewhat 

by  the  spatial  dependence  of  A  on  the  scale  of  the 

L  7  18 

coherence  length  near  the  S/N  interface.  ’ 

A  similar  situation  should  prevail  in  the  phase-siip 

1 9 

center,  and  recent  experiments  by  Dolan  and  Jackel  have 
provided  rather  direct  evidence  for  the  SBT  model. 

Samples  were  fabricated  which  included  arrays  of  very 
small  voltage  orobes,  closely  spaced  along  the  length  of 
a  t.hin-film  micrcstrip  (see  Fig.  1.4a).  Both  super¬ 
conducting  and  normal  probes  were  used,  which  made  it 
possible  to  make  independent  measurements  of  p  and  p 
as  a  function  of  position  along  the  strip.  As  shown  in 
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Fiq.  1.4b,  the  sDatial  behavior  of  these  t ime-averaged 
potentials  was  found  to  go  exactly  as  expected.  The 
value  of  jT  jumped  abruptly  at  the  center  of  the  dissi¬ 
pative  reqion,  while  u  relaxed  exponentially  back  to 
the  equilibrium  value  on  each  side  of  the  center  with 
a  decay  lenqth  A  which  fit  the  form  A  =  *TTi(1i>  . 

This  identifies  the  diffusion  lenqth  A  in  the  SBT 

model  with  the  lenqth  Ay* .  The  weak  diverqence  in 

Ay*  «  (TC~T)  ^  as  T *T^  was  not  apparent  in  many  of  the 

earlier  measurements  on  I-V  curves  of  PSC’s,  so  that 

SBT  made  the  initial  identification  of  A  with  the  non- 

diverqent  distance  »  P;  ,.  Some  reasons  for  the  non- 

F. 

observance  of  this  diverqence  include  the  nonuniformity 
of  the  critical  temperature,  and  the  fact  that  near  T 
the  onset  of  a  given  PSC  becomes  rounded,  and  becomes 
less  clearly  separated  from  the  onset  of  the  next  PSC. 
Furthermore,  farther  from  T  the  power  dissipated  IV 
becomes  significantly  larqer,  producing  "local  heating" 
at  the  PSC,  and  this  distorts  the  I-V  curve  by  adding 
upward  curvature  to  the  plateau  region  and  creating 
hysteresis  between  the  zero-voltage  and  the  dissipative 
branches  of  the  curve.*  Nevertheless,  within  the 
limited  temperature  range  over  which  changes  in  the 
differential  resistance  of  a  PSC  accurately  reflected 
changes  in  the  decay  length,  Dolan  and  Jackel  did  find 
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that  this  resistance  corresponded  to  the  normal  resis¬ 
tance  of  the  lenqth  equal  to  twice  this  cjuasipart icle 
diffusion  lenqth  . 

The  earlier  experiments  on  PSC's  were  done  in  zero 
applied  maqnetic  field.  Like  other  pair-breakinq  inter¬ 
actions  such  as  maqnetic  impurities  and  larqe  current 
densities,  a  maqnetic  field  parallel  to  a  thin-film 
sample  not  only  reduces  the  critical  temperature,  but 
also  induces  characteristic  chanqes  in  the  density  of 
states  and  the  excitation  spectrum."0  For  larqe  pair- 
breakinq  this  leads  to  "napless  superconductivity", 
characterized  by  the  absence  of  a  qap  in  the  excitation 
spectrum  while  maintaining  superconductinq  order. 

The  experimental  research  to  be  discussed  in 
Chapter  Two  was  undertaken  to  determine  to  what  extent 
the  SBT  model  of  the  PSC  remains  valid  in  the  presence 
of  a  parallel  maqnetic  field,  and  what  effect  the  field, 

as  an  example  of  a  pair-breakinq  perturbation,  has  on  the 

2 1 

nonequilibrium  processes  in  the  superconductor.  As 
illustrated  in  Fiq.  1.5,  we  find  that  PSC's  are  indeed 
observed  for  T  near  Tc<H)  (the  reduced  critical  tempera¬ 
ture),  and  that  the  relaxation  time  r  (t,h)  inferred 

from  the  increments  R  of  plateau  differential  resistance 

n 

depends  stronqly  on  the  intensity  of  the  field.  By 
payinq  special  attention  to  the  important  role  played  by 
heatinq  in  oroducinq  characteristic  distortions  of  the 
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Fiq.  1.5  I-v  curves  for  a  Sn  microstrip  (025A  in  Table 

2.1)  with  and  without  a  parallel  maqnetic  field. 
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I-V  curves,  we  are  able  to  minimize  and  partially 

correct  for  these  effects.  We  are  then  able  to  fit  the 

inferred  relaxation  time  to  an  expression  derived  by 
..  12 

Schmid  and  Schon  for  transverse  mode  relaxation  in 
the  presence  of  pair-breaking,  in  effect  a  generalized 


We  have  been  emohasizing  the  similarities  between 
the  PSC  and  other  situations  involving  nonequilibrium 
potentials  in  superconductors,  in  particular  quasi¬ 
particle  injection  and  the  N/S  boundary.  There  are 
important  differences,  however,  chief  among  them  being 
that  the  latter  processes  involve  a  steady  state, 
whereas  the  PSC  is  a  dynamic  t ime-deDendent  process, 
in  which  the  gap  oscillations  in  the  core  of  the  PSC 
may  be  expected  to  drive  oscillations  in  the  currents 
and  potentials.  A  connection  can  be  made  between  these 

oscillations  and  the  subject  of  collective  modes  in 
->2 

superconductors.  One  of  these  modes,  first  detected 
experimentally  in  superconducting  fluctuations  below  T  ,‘ 
consists  of  propagating  but  attenuated  waves  involving 
a  counterflow  of  normal  and  supercurrent.  As  we  demon¬ 
strate  theoretically  in  Chapter  Three,  these  waves  can 

•>  4 

be  derived  from  a  generalized  two-fluid  picture"  as  an 
ac  generalization  of  the  steady-state  charge  imbalance 
effects  discussed  earlie*:.  We  show  further  that  these 


charqe  imbalance  waves  are  analogous,  for  a  one-dimen¬ 
sional  superconductor  near  T  ,  to  electrical  signals 

25 

propagating  down  a  simt'le  transmission  line.  This 
transmission  line  provides  a  convenient  picture  in  which 
to  view  the  whole  range  of  experiments  involving  charge 
imbalance  effects  in  superconductors,  and  we  make  some 
predictions  of  additional  effects  which  have  not  yet 
been  observed  experimentally. 

We  proceed  with  our  theoretical  study  in  Chapter 

Four,  with  a  model  for  the  PSC  in  which  the  essential 

features  of  the  FBT  model  are  generalized  to  account 

for  ac  effects,  and  the  temporal  and  spatial  evolution 

of  the  currents  and  potentials  can  be  understood  in 

terms  of  charge  imbalance  waves  driven  by  Josephson 

2  ^ 

oscillations  in  the  core  of  the  PSC.  This  can  be 
represented  as  an  ideal  Josenhson  element  of  negligible 
extent  (corresponding  to  the  core)  coupled  to  the 
transmission- 1 i ne  equivalent  circuit  mentioned  above. 

We  then  take  advantage  of  standard  transmission-line 
theory  to  calculate  a  st raight foward  solution  to  the 
problem.  In  addition  to  providing  a  consistent  picture 
of  the  ac  properties  of  PSCs,  this  solution  also  predicts 
the  existence  of  intrinsic  hysteresis  in  the  dc  1 -V  curves, 
quite  apart  from  hysteresis  due  to  heating,  which  we 
have  not  included  in  this  model. 


As  in  any  ,vt i ve  field  of  research,  the  present 
work  does  not  answer  all  the  questions  or  follow  up  on 
all  possible  leads.  In  Chapter  Five  we  discuss  some 
additional  topics,  which,  although  incomplete,  may  be 
a  guide  for  future  research.  These  include,  on  the 
experimental  side,  a  more  detailed  discussion  of  some 
effects  related  to  localized  weak  spots  in  magnetic 
fields,  in  particular  critical  currents  and  low-vo 1 taae 
behavior.  On  the  theoretical  side,  we  discuss  some 
preliminary  calculations  concerning  the  application 
of  our  PSC  model  to  the  phenomena  of  synchronization 
and  resonance  in  series  of  two  or  more  PSC's. 

To  summarize  the  structure  of  the  present  report , 
we  continue  in  Chapter  'T'wo  with  a  presentation  and 
discussion  of  our  experimental  measurements  on  PSC's 
in  parallel  magnetic  fields.  Chapter  Three  contains 
a  general  theoretical  analysis  of  charge  imbalance 
dynamics ,  which  is  used  in  Chapter  Four  to  develop 
our  model  for  the  PSC.  Some  additional  tooics,  both 
experimental  and  theoretical,  are  discussed  in  Chapter 
Five,  and  in  Chapter  Six  we  summarize  our  current 
understanding  of  noneou i 1 ibr ium  phenomena  in  super¬ 
conducting  phase-slip  centers. 


CHAPTER  TWO:  PHASE-SLIP  CENTERS  IN  MAGNETIC  FIELDS 


2 . 1  Introduct ion 

We  present  in  this  chapter  a  review  of  our  experi¬ 
mental  studies  on  lonq  tin  superconducting  microbridges 

in  the  presence  of  parallel  magnetic  fields,  following  up 

1  2  6 

on  the  work  of  Skocpol  et_  a_l_.  '  on  such  samples  in  zero 

field.  Much  of  this  work  has  been  reported  in  previous 
works  of  the  author. 

We  show  in  Fig.  2.1  a  set  of  I-V  curves  of  a  tin 
microstrip  in  a  parallel  magnetic  field,  for  a  set  of 
temperatures  near  the  reduced  critical  temperature  T  (H) . 
As  we  noted  in  Chapter  I,  the  general  appearance  of 
these  curves  is  similar  to  that  in  zero  field.  The 
voltage  "step"  may  be  somewhat  rounded,  but  it  is  followed 
in  the  usual  way  by  a  region  of  approximately  constant 
differential  resistance.  In  this  sample  we  have  used  a 
weak  spot  in  the  center  of  the  bridge  to  effectively 
isolate  the  curve  for  a  single  PSC,  from  which  we 
extract  our  data.  The  differential  resistance  of  the 
PSC  increases  at  high  voltages,  which  as  we  will  show 
later  is  due  to  heating.  It  is  also  larger  at  low 
voltages  because  of  the  rounded  onset  of  the  plateau 
region.  We  therefore  choose  the  minimum  value  of  dv/dl 
(the  slope  of  the  dashed  lines  tangent  to  the  curves  of 
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Fig.  2.1  V ( I )  for  Tin  microbridge  #15A  in  a  parallel 

field  of  1200  Oe,  for  a  range  of  temperatures  near  Tc(H) 

Curve  A  -  T=3.709  K;  B  3.684;  C  3.649;  D  3.628;  E  3.595; 

F  3.568;  G  3.535;  H  3.505;  I  3.472;  J  3.434.  The  slopes 

of  the  dashed  lines  tangent  to  the  curves  are  taken  as 

the  characteristic  resistance  R  of  the  PSC.  Curve  C  - 

n 

Rn=  0.83  $1;  D  .793;  E  .727;  F  .678;  G  .648;  H  .677; 

I  .685;  J  .683. 
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Fig.  2.1)  as  the  characteristic  resistance  R^  of  the 

nonequilibrium  region  of  the  bridge.  Following  SBT,  we 

then  infer  the  relaxation  time  t  using  the  relation 

K 

R  =  |(2A)  =  ^Vd7~  (2.1) 

n  A  A  R 

where  p  is  the  normal-state  resistivity,  A  is  the  cross- 
sectional  area,  A ( T , H )  is  the  experimentally  inferred 
quasiparticle  diffusion  length,  and  D  is  the  normal-state 
diffusion  constant. 

Before  discussing  the  data  we  obtained,  we  will 
review  in  Section  2.2  the  fabrication  and  characteriza¬ 
tion  of  our  thin-film  microbridges  and  our  experimental 
procedure.  In  Section  2.3  we  oresent  our  measurements 
of  the  temperature  and  magnetic  field  dependence  of  the 
inferred  transverse  relaxation  time,  and  compare  them  to 
theory.  In  Section  2.4,  we  discuss  heating,  on  both  a 
microscopic  and  a  phenomenological  level,  and  explain 
how  distortions  in  the  I-V  characteristics  produced  by 
heating  can  be  minimized  and  partially  corrected  for. 

In  Section  2 . 5  we  discuss  how  the  use  of  a  deliberate 
weak  spot  to  isolate  a  single  PSC  and  reduce  heating 
allowed  vis  to  obtain  more  accurate  data  over  a  wider 
range  of  temperatures  and  fields.  Finally,  in  Section 
2.6,  we  summarize  our  current  understanding  of  phase-slip 
centers  and  nonequilibrium  relaxation  times  in  parallel 
magnetic  fields. 


2 . 2  Samples  and  Hxper iment a  1  rrocodu res 

The  procedures  for  sample  preparation  were,  except 

where  otherwise  noted,  very  similar  to  those  used  by 
26 

Skocpol  for  his  measurements  on  Sn  microbndqes  in 
zero  magnetic  field,  and  we  will  only  review  them  here. 
The  samples  were  prepared  on  glass  or  crystal  sapphire 
substrates,  typically  1  *  V  in.  in  size,  to  which  leads 
were  first  attatched  with  small  patches  of  In  or  In-Sn 
alloy.  A  film  of  Sn  was  then  vapor -deposited  on  a  cold 
(77k)  substrate  in  a  vacuum  of  better  than  10  *'  torr, 

O 

and  the  thickness  (of  order  1000  A)  was  measured  with  a 
calibrated  quartz  crystal  monitor  to  an  accuracy  of 
better  than  10%.  By  varying  the  deposition  rate,  the 
purity  of  the  tin  film  and  hence  the  residual  resistivity 
could  be  varied  substantially .  A  typical  rate  for  a 

O 

reasonably  clean  film  was  100  A  sec.  In  some  samples, 

the  Sn  film  was  applied  on  top  of  a  thin  undercoat  of 

Sn  which  had  been  oxidized  by  a  glow  discharge  in  0,, 

in  order  to  achieve  better  adhesion  of  the  film  to  the 
27 

subst rate . 

The  critical  temperature  of  pure  bulk  Sn  is  1.7.'  K, 

but  that  of  thin  films  is  generally  found  to  be  higher. 

This  enhanced  T  lias  been  shown  to  be  associated  with 

increased  elect ron-phonon  coupling  at  thin-film  surfaces 

2  8 

and  crystallite  boundaries."  The  values  of  T ,  for  our 
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samples  (see  Table  2.1)  were  typically  in  the  range 
from  3.8  to  3.9  K,  with  those  for  some  extremely  thin 
samples  up  to  4.0  K. 

With  the  substrate  mounted  on  an  adjustable  micro¬ 
scope  stage  and  a  diamond  knife  held  above,  a  small 
region  of  the  film  together  with  the  edge  of  the  knife 
were  brought  into  focus  under  a  microscope.  The  knife 
was  carefully  lowered  onto  the  surface,  and  adjacent 
regions  of  the  film  were  scraped  away,  leaving  the  long, 
thin  rectangular  microbridge  connecting  the  two  halves 
of  the  film.  Dimensions  ranged  from  1  to  5  pm  in  width 
and  from  15  to  150  pm  in  length.  The  parameters  for  a 
number  of  bridges  are  listed  in  Table  2.1. 

In  general,  our  samples  were  not  quite  the  ideal 
quasi-one-dimensional  filaments  that  our  discussion 

might  indicate.  In  the  absence  of  a  magnetic  field, 

29 

the  GL  coherence  length  takes  the  form 

C  (T)  =  0.85  (C  «.)^(1-T/T  )~*  (2.2) 

o  c 

for  a  dirty  superconductor  near  T  ,  where  l  is  the 
mean  free  path  and  £  =  hv_/TrA(0)  is  the  BCS  coherence 

O  r 

length.  A  similar  formula  using  the  renormalized  Tc(H) 
will  apply  in  the  presence  of  a  parallel  magnetic  field. 
(This  is  only  valid  if  the  transition  to  the  normal  state 
is  second  order  in  the  presence  of  the  field.)  Choosing 

O 

typical  parameters  l  =  700  A  and  Tc~T  =  40  mK,  one  obtains 
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(T)  1  iim.  Although  the  thickness  d  of  our  samples 

was  always  much  less  than  this,  the  width  was  often 

somewhat  greater.  This  was  in  part  because  submicron- 

width  bridges  are  more  difficult  both  to  fabricate  and 

to  maintain  (because  of  burnout  problems)  than  the  wider 

ones.  In  the  limit  of  very  large  width  (w'-'"r,  (T)  )  ,  the 

current-induced  breakdown  of  superconductivity  has  been 

shown  to  occur  via  the  nucleation  of  vortices  at  the 

edges  of  the  film,  which  then  move  across  the  film.^0 

It  seems  that  the  transition  between  narrow  and  wide 

strips  is  sufficiently  broad  that  the  SBT  Picture  of 

a  PSC  still  holds  approximately  in  the  range  w  i  £  (T)  .  ^ 

A  further  assumption  is  that  the  applied  field 

penetrates  the  film  uniformly.  This  is  valid  if  the 

thickness  d  is  much  smaller  than  the  magnetic  penetration 

32 

depth  in  zero  field 

MT)  ~  \T(0)  (£  /l)^  [1-('T’/t  )4]~*  (2.3) 

L  o  c 

o 

where  \  (0)  (=  340  A  for  Sn)  is  the  London  penetration 

depth  at  T=0 .  In  the  opposite  limit,  if  d"V5  \,  the 
film  undergoes  a  first  order  phase  transition  at  the 
critical  field  Hc||  (T)  ,  i.e.  the  gap  goes  discontinuously 
to  zero  at  the  transition.  All  our  samples  are  thin 
enough  to  have  a  second  order  transition,  so  that  the 
parallel  critical  field  is  given  by ^ 


N 
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Hc||(T)  =  2s6  Hcb(T)  \<T)/d  (2.4) 

which  is  always  larger  than  the  bulk  critical  field  (T) 
For  tin,  Hcb(0)  *  308  qauss ,  and  we  have  studied  samples 
with  the  extrapolated  (0)  from  about  700  qauss  to  over 
4  kiloqauss.  In  those  samples  where  d''\,  the  gap  in  the 
presence  of  the  field  satisfies  the  relation*1 


nart-way  through  the  thickness  of  the  film  ("nick". 

Fig.  2.2b).  This  was  a  rather  touchy  procedure,  and 
required  considerable  trial  and  error  as  well  as  frequent 
trios  between  the  low-oowered  microscope  used  for 
fabrication  and  the  hiqh-powered  one  (('00')  used  for 
careful  samole  examination.  The  width  of  one  such  cut 


was  determined  from  a  scanning  electron  micrograph  to 

O 

be  about  1000  A,  far  less  than  A  under  any  conditions. 


& 
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(a) 


Typical  geometries  for  bridges  with  weak  spots, 
a)  Overall  view,  b)  and  c)  are  magnified  side 
views  of  the  "nick"  and  "groove"  types  of  weak 
spots,  respectively. 


Fig.  2.2 
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This  type  of  weak  spot  tended  to  heal  itself  in  a  period 

of  order  a  few  weeks  at  room  temperature.  We  conjecture 

that  the  mechanism  for  this  may  be  thermally  activated 

surface-self  diffusion,  in  part  because  the  healing  of 

small  scratches  on  a  metal  surface  has  actually  been 

34 

used  to  study  surface  self-diffusion  in  metals.  In 
other  cases,  the  "weak"  spot  was  produced  by  depositing 

O 

the  film  over  a  groove  in  the  substrate  (width  <  5000  A)  , 

and  then  fabricating  the  bridge  transverse  to  this  groove 

("groove",  Fig.  2.2c).  The  local  reduction  of  the 

critical  current  at  the  groove  was  a  product  of  its 

interaction  with  the  external  magnetic  field,  as  will  be 

discussed  in  Section  2.5,  rather  than  a  simple  decrease 

in  geometrical  cross-section.  The  fabrication  of  this 

groove  was  made  with  the  sharp  pointed  tip  of  a  diamond 

knife,  in  a  similar  way  to  analogous  grooves  made  by 

27 

M.  Octavio  for  his  variable-thickness  bridges.  This 
form  of  weak  spot  seemed  to  remain  relatively  stable, 
presumably  because  of  the  adhesion  to  the  substrate 
at  the  lower  surface. 

In  mounting  the  sample  onto  the  cryostat,  care  had 

to  be  taken  to  prevent  burnout  by  electrical  transients, 

although  the  problem  was  generally  not  as  severe  as  that 

27 

experienced  by  Octavio  with  very  small  microbridges. 
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Burnout  involved  melting  and  drawing  back  of  the  tin  in 
the  weakest  region  of  the  bridge,  and  was  normally 
clearly  visible  under  the  microscope.  To  assist  in 
preventing  it,  a  conducting  path  parallel  to  the  bridge 
was  left  connected  until  after  the  sample  leads  were 
soldered  to  the  voltage  and  current  leads  on  the  cryo¬ 
stat.  The  soldering  iron  was  unplugged  for  this  operation, 
since  otherwise  inductive  coupling  to  the  heating  element 
produced  ac  voltages  of  order  10  mV  in  the  hot  tip. 

After  scratching  away  the  conducting  shunt,  the  continuity 
of  the  bridge  was  checked  using  a  1  uA  source  and  a 
microammeter,  rather  than  an  ohmmeter  which  might  cause 
burnout.  Most  burnouts  tended  to  occur  at  room  temperature, 
but  transients  in  the  electrical  measuring  system  (e.g. 
due  to  turning  equipment  on  or  off)  have  also  caused 
problems  on  occasion. 

^he  sample  was  mounted  vertically  in  a  horizontal 
magnetic  field,  and  was  either  directly  immersed  in  the 
liquid  helium  bath,  with  pressure  regulation  to  fix  the 
temperature,  or  inside  a  vacuum  can,  with  a  heater  and 
electronic  temperature  regulation  (see  Appendix  1  for 
the  details  of  the  temperature  regulation  circuit). 

Bath  temperatures  could  be  maintained  in  the  range  from 
4.2  K  down  to  about  1.3  K.  The  temperature  was  measured 
with  a  calibrated  doped-Ge  resistance  thermometer  to  a 


relative  accuracy  of  about  1  mK.  Farlier  measurements 

taken  with  a  carbon  resistance  thermometer  were  found  to 

be  insufficiently  reproducible.  The  magnet,  with 

conventional  pole-pieces  located  outside  the  dewar, 

was  feedback-stabilized  with  relative  accuracy  of  about 

1  gauss.  Fields  of  up  to  several  kiloaauss  were  used. 

At  this  level  magnetoresistance  in  the  Ge  thermometer 
2 

(which  goes  as  H  )  produced  errors  of  at  most  several  mK , 

which  could  be  ignored.  The  magnet  could  be  rotated 

around  the  dewar  on  a  track,  so  that  one  could  align  the 

field  parallel  to  the  plane  of  the  thin-film  sample. 

This  was  accomplished  to  within  a  quarter  of  a  degree 

by  maximizing  I  at  fixed  T  and  large,  fixed  H.  This 

maximum  reflects  the  fact  that  the  critical  field  of 

a  film  is  a  sharply-peaked  function  of  angle,  with  a 

maximum  for  the  parallel  orientation  and  a  minimum  for 
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perpendicular  orientation.  It  is  only  for  the  parallel 
orientation  that  the  applied  field  causes  the  gap  to  be 
depressed  continuously  and  uniformly  to  zero;  a  perpen¬ 
dicular  component  will  cause  the  entry  of  vortices  or 
flux  bundles. 

In  taking  the  I-V  curves,  the  sample  was  current- 
biased  with  a  battery-powered  source,  and  the  voltage 
was  measured  with  a  Keith  ley  nanovoltmeter  (or  at  times, 
a  PAR-113).  In  some  cases,  a  PAR  MR-8  lock-in  amplifier 


was  used  to  measure  dV/dl  directly  as  a  function  of 
current.  The  curves  were  nlotted  on  an  X-Y  recorder 
as  the  bias  current  was  swept  by  hand.  The  schematic 
circuits  for  these  situations  are  shown  in  Fig.  2.3. 

In  order  to  reduce  rf  interference,  the  leads  con¬ 
sisted  of  twisted  wire  oairs,  and  a  set  of  Erie  rf 
filters  were  also  used.  The  measurements  were  not , 
however,  taken  in  an  rf-shielded  room.  Mechanical 
stability  of  the  magnet  dewar  was  also  important,  since 
one  source  of  low-frequency  interference  was  inductive 
pick-up  due  to  vibrations  of  the  cryostat  in  large 
magnetic  fields. 


DC/AC  BIAS 
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2 . 3  Behavior  of  the  Transverse  Relaxation  Time  x_R 

As  we  have  noted  earlier,  the  I-V  characteristics 

of  a  tin  microstrip  in  a  parallel  magnetic  field,  for  T 

near  Tc(H)  *  are  qualitatively  the  same  as  those  for  H=0 

near  Tc  (see  Fig.  1.5).  Quantitatively,  the  major 

effect  of  increasing  the  magnetic  field  is  the  decrease 

of  R  ,  the  differential  resistance  of  a  single  PSC  in 
n 

the  plateau  region,  together  with  the  corresponding 

increase  in  the  number  of  PSC's  which  can  "fit"  along 

the  length  of  the  bridge.  There  is  also  a  variation 

of  Rn  with  temperature  for  constant  H  (see  Fig.  2.1), 

which  consists  of  a  weak  divergence  near  Tc(H)  and  an 

enhancement  due  to  heating  farther  away.  Still,  one  can 

take  the  minimum  value  of  R  for  a  given  field  H  as  a 

measure  of  the  general  field  dependence  of  R^  and  the 

relaxation  time  Tn  inferred  from  it.  In  Fig.  2.4,  these 

values  of  x  are  plotted  against  the  reduced  magnetic 

field  H/Hc||  (0)  for  the  first  PSC  of  a  series  of  different 

samples.  Typically,  the  length  scale  decreases  from  about 

5  urn  for  H=0  to  less  than  1  urn  for  large  fields,  and  the 

-9 

corresponding  times  decrease  from  about  10  sec  to 
almost  10  ^  sec. 

The  experimental  values  of  the  relaxation  time  were 
inferred  from  the  measured  values  of  the  differential 
resistance  R  usina  the  relation 

n 
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Fig.  2.4  Relaxation  time  t  versus  reduced  magnetic  1 

H 

field  H/H^i 10) .  Each  point  corresponds  to  the  minimum 

value  of  R  of  the  first  PSC  for  T  near  T  (H) .  This 
n  c 

shows  the  aeneral  magnetic  field  dependence  of  r  , 

R 

neglectinq  the  divergence  at  T  (H) .  The  experimental  1 

c 

points  are  compared  with  three  characteristic  times, 

which  are  defined  as  in  the  text  with  t„(T  )  =  3  "  10-10 

E  c 

sec .  | 

! 


N 


Samples 
3 


•  15  A 
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Here  d  is  the  film  thickness,  w  is  the  bridge  width,  and 
p  is  the  normal-state  remnant  resistivity.  This  resist¬ 
ivity  was  determined  either  from  the  resistivity  ratio  of 
the  sample  between  room  temperature  and  4.2  K  or  from 
the  bridge  geometry  and  its  4.2  K  resistance.  The  geometry 
was  measured  from  optical  or  scanning  electron  micrographs. 

The  diffusion  constant  D=v„t/3  was  determined  using  the 

r 

-11  2 

relation  p£  =  1.0  *  10  Q-cm  together  with  the  value 
8  3  6 

vF  =  0.7  x  10  cm/sec.  We  estimate  that  taking  into 
account  the  various  uncertainties  in  all  these  measure¬ 
ments,  the  proportionality  constant  between  tr  and  FT 
should  be  known  for  a  given  sample  to  within  a  factor 
of  two. 

In  zero  field,  SBT  made  the  tentative  identification 

t0  -  t_(T  ),  where  x_  is  the  electron-ohonon  collision 
R  E  c  E 

time.  In  the  absence  of  a  comprehensive  understanding 
of  the  problem,  we  initially  expected  that  the  effect  of 
a  magnetic  field  would  be  to  replace  Tc  by  the  reduced 
TC(H),  yielding  tr  =  t  [T  (H)].  Since  the  phonons 
"freeze  out"  at  low  temperatures,  this  time  gets  sub¬ 
stantially  longer:  t  T  3 .  This,  however,  is  in 

Jj 

strong  disagreement  with  the  experimental  data  (see  Fig. 
2.4);  the  basic  trend  cioes  in  the  other  direction. 
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There  is,  however,  another  characteristic  time  in 
the  problem.  A  parallel  magnetic  field  in  a  dirty  super¬ 
conductor  (£<<£  )  is  one  member  of  the  class  of  pair- 
o 

breaking  perturbations,  which  in  the  limit  of  strong 

20 

perturbation  leads  to  gapless  superconductivity . 

The  effect  of  any  of  these  perturbations  on  the  gap  and 
the  excitation  spectrum  can  be  characterized  in  terms  of 
the  pair-breaking  energy  2ct  ,  or  equivalently  the  pair¬ 
breaking  time  Ts=fi/(2a),  which  for  a  magnetic  field 

20 

parallel  to  a  thin  film  take  the  form 

For  large  pair-breaking,  superconductivity  will  be 
totally  destroyed  (i.e.  the  critical  temperature  driven 
to  zero)  when  this  pair-breaking  energy  is  equal  to  the 
energy  gap  in  the  absence  of  pair-breaking: 

2  ot  ( 0 )  =  A  (T=0  ,H=0 )  =  1.76  kDT  (2.8) 

c  Be 

The  zero-temperature  critical  field  can  therefore  be 
expressed  as 


t*  ( r,  \  _  -s 

h^ll  ~  — nSTa's - 


so  that  we  can  write 


Ti 

A (0,0) 


(2.9) 


(2.10) 


For  small  pair-breaking,  the  depression  of  the  critical 
temperature  is  linear  in  the  pair-breaking  energy, 
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T  - 
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.(H) 
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4k 


(2.11) 


B 


so  that  the  nai r-break inq  tine  can  also  be  written 


_  ft 

s  "  ^kB(Tc-Tc(!l)  ] 


(2.12) 


(cf.  the  Ginzburg-Landau  time  i  =  uft/8k. (T  -T)  ) 

(»Li  Be 

and  the  temperature  dependence  of  the  critical  field 
takes  the  form 


|,c||(-)  =  1.2  ll.i  (0)  (1-T/r.) 


(2.13) 


Fq .  (2.13)  is  approximately  valid  down  to  about  %TC  • 

In  practice,  we  use  Fn.  (2.13)  to  determine  l!^  (0)  , 
by  fitting  to  the  measured  values  of  T  (H) .  Since  we 
want  this  to  correspond  to  the  critical  temperature  of 
the  uniform  film  and  not  that  of  the  weak  spot  (which  is 
somewhat  depressed  in  many  of  our  samples),  this  ^.(’l) 
is  determined  by  the  onset  of  an  essentially  normal 
linear  I-V  characteristic. 

If  other  pair-breaking  mechanisms  are  oresent , 
the  total  pair-breaking  energy  is  the  sum  of  the  indi¬ 
vidual  ones.  Collisions  with  phonons  also  break  pairs, 

12 

so  that  the  comnlete  nair-break ing  time  is 


i  +  - 

TS  2  1  F 


(2.14) 


which  crosses  over  from  2 1 ,,  in  zero  field  to  i  in  large 

F  s 


fields.  As  one  can  see  in  Fig.  2.4,  1 


n-1 


is  too  small 


and  has  the  wrong  field  dependence  to  fit  the  data  for  i 
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But  the  qeomotric  mean  between  and  r  turns  out 

F 

to  have  the  riqht  order  of  magnitude,  and  furthermore, 

seems  to  be  iustified  theoretically  in  the  work  of 

l'>  37-39 

Schmid  and  Schon  and  others 

The  relaxation  time  in  the  presence  of  pair-breaking 

derived  by  Schmid  and  Schon  for  a  homogeneous  transverse 

mode  disequilibrium  (which  can  be  produced  by  uniform 

quasiparticle  injection  into  a  film)  is 


4k  T 
(T) 

tf. 

n-r  ' 

1  4- 

A 

1 

3 

C> 

,2r 

X  • 

2A?rF 

where  A=  A(T,H)  and  tp,  =  Tr(T).  For  our  experiments 
the  final  factor  is  close  to  unity,  so  that 


tr(T,H) 


4y 

TT  A 


2r 


4kRT 

trA 


IT  A 


Q  ,  for  t  g  >  ^  t  R 


(2.16) 


<ys>  ,  for  >S-<'E 


Note  that  this  reduces  correctly  to  for  H=0  (the 
above  relaxation  time  is  in  effect  a  generalized  t  ) , 
while  in  the  hiqh-field  limit  (typically  H-Hc|j  (0)  /10) 
one  gets  the  geometric  mean  between  t„  and  F  * ,  with  a 

r. 

divergence  as  A  ^  at  T  (H) .  Similar  expressions  to 

c 

Eg.  (2.15)  have  been  derived  by  Artemenko  et  al_. '  , 

3S  39 

Shelankov,  and  ^nt  in-Wohlman  and  Orbach  . 

The  temperature  and  field  behavior  of  R^  for  two 
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different  samples  are  shown  in  Figs.  2.5  and  2.6. 

The  values  of  Rn  are  plotted  as  a  function  of  bath 
temperature  T  for  various  constant  values  of  parallel 
magnetic  field  H.  In  Fig.  2.5,  we  also  include  values 
of  R  corrected  for  heating  effects,  as  will  be  discussed 
in  Section  2.4.  Also  plotted  (the  solid  lines)  are  the 
corresponding  theoretical  values  from  Fa.  (2.16),  taking 

te(T)  =  te(Tc) (Tc/T) 3  (2.17) 

A  (T  ,H)  =  A(T)  [1-H2/H^||  (T)  ]  ^  (2.18) 

Rn  =  ( 2  p/wd)  [Dt^T) (2.19) 

and  treating  r  (T  )  as  a  fitting  parameter  for  all  the 
data  from  the  particular  sample.  We  find  rather  good 
agreement,  considering  all  the  uncertainties.  This 
agreement  is  achieved  using  values  of  t  (T  )  in  the* 
range  from  1  to  5  «  10  18  sec  for  all  our  samples. 
Variation  within  this  range  does  not  seem  to  correlate 
with  any  particular  samnle  parameter,  indicating  that 
it  may  be  due  to  uncertainties  in  some  of  the  experimental 
constants  (e.g.  sample  geometry  and  mean  free  path)  in 

/  rp\ 

Eq.  (2.19)  which  relate  and  .  However,  all  of 

the  values  of  t  (T  )  which  we  have  inferred  are  in 
I’j  c 

reasonable  accord  with  other  experimental  and  theoretical 
estimates  of  the  same  constant  for  tin . 18  '  ' 36 ' 40 

For  example,  Clarke  and  Paterson36  obtained  1.4  x  10 


- . 
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T(K) 

Fig.  2.5  Differential  resistance  of  the  first  PSC 

in  tin  bridge  #15?,  for  various  constant  values  of  H  ||J, 

as  a  function  of  T.  ...  T’he  solid  curves  represent  the 

bath 

theoretical  predictions  of  Fqs.  (2.17)  -  (2.19),  with 
Tj,(Tc)  =  3  *  10  ^  sec.  The  solid  symbols  include  correc¬ 
tions  for  heating-induced  distortion  (Section  2.4). 
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1 


4  0  -10+  l1 

sec  *  and  Kaplan  et  al.  calculated  2.8  x  10  sec. 

The  complete  form  of  F.q.  (2.16)  is  difficult  to 
understand  in  detail,  but  one  can  gain  a  partial  under¬ 
standing  by  considering  the  physical  processes  involved. 

The  basic  picture  seems  to  be  that  when  t  <<tp,  pair- 
breaking  interactions  allow  for  very  rapid  branch-crossing 
for  low  energy  quasiparticles.  If,  however,  a  signifi¬ 
cant  amount  of  the  branch  imbalance  is  being  created  at 
high  energies,  the  excess  quasiparticles  must  first  be 
relaxed  by  inelastic  scattering  into  the  low  energy 
region,  so  that  the  overall  relaxation  rate  will  involve 
both  t  and  x_.  This  interpretation  has  been  emphasized 

S  fc. 

by  Shelankov,*®  by  Clarke  and  Schmid, and  by  Entin- 

39  I 

V’ohlman  and  Orbach. 

.  .  4  2 

Within  this  picture,  Clarke  and  Schmid  have 

suggested  a  heuristic  method  for  deriving  the  general 

form  of  Eq.  (2.16)  in  the  limit  t  <<i_.  Consider  a  BCS 

s  E 

superconductor ,  with  quasiparticle  scattering  events 

-  - -  x 

*  Their  result  was  expressed  in  the  form 

tq(T)  =  1.0  x  lo-10  sec  [ A ( 0 ) /A (T) ]  . 

4 

x  They  computed  a  different  characteristic  electron- 

_  Q 

phonon  time,  t  =  2.3  x  10  sec,  which  is  related 
o 

to  t  by  the  relation4*  t_ (T  )  =  t  /7£(3)  =  t  /8.4  . 
r»  t,  c  o  O 


3 
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across  the  Fermi  surface  induced  by  the  pair-breakinq 

interaction,  ''’he  transition  rate  is  proportional  to 

4 1 

the  coherence  factor 

( uu 1  m')2  =  %(1  +  -  fp-r)  (2.20) 

where  f.  is  the  normal  electronic  enerqy  referred  to  the 
condensate  chemical  potential  and  r  =  / f,  ’  +  A •'  is  the 
quasiparticle  excitation  energy  (a  pair-breaking  inter¬ 
action  due  to  maonetic  fields  or  paramagnetic  impurities 

4  3 

is  a  "Type  II"  perturbation  ) .  For  scattering  straight 

across  from  one  branch  to  the  other,  \  '  =  -f. ,  and  the 

coherence  factor  is  A2/E:.  This  is  reasonable,  since  a 

branch-crossing  transition  is  forbidden  in  a  normal 

metal  where  A=0.  The  branch-crossing  transition  at  a 

given  energy  will  proceed  at  a  rate  —  A-’/E',  and  if  this 

s 

is  faster  than  the  rate  at  which  additional  nonequilibrium 
quasiparticles  can  scatter  down  from  higher  enerqies  , 
the  latter  rate  will  determine  the  net  rate  of  branch 
imbalance  relaxation.  The  cross-over  between  these  two 
rates  occurs  at  an  energy 

r-0  *  a/T^TT  (2.2D 

and  if  we  assume  that  the  noneouilibrium  quasiparticles 
are  created  with  an  initially  thermal  distribution, 
a  fraction  ''Eo/k^T  of  the  inelastic  scattering  events 
will  involve  scattering  into  the  lower  energy  region 
and  consentient  branch  imbalance  relaxation.  The  net 


A  problem  with  this  argument  is  that  it  is  only 

(T) 

valid  when  F  is  less  that  k_T,  i.e.  >  t„.  This 

O  9  K  h 

condition  is  seldom  satisfied  and  usually  strongly 

violated  in  our  experiments.  It  is,  however,  the  same 

condition  under  which  Schmid  and  Schon  originally  derived 

Eq.  (2.15).  The  oreceeding  argument  seems  to  suggest 

that  if  E  >  k  T,  then  the  overall  relaxation  rate  will 
O  D 

1^2  30 

be  of  order  —  <f)~z'  *  Nevertheless,  Shelankov  has 

s  3 

proposed  conditions  whereby  Eg.  (2.16)  should  be  valid 
over  a  wide  range  of  fields,  which  seem  to  be  approx¬ 
imately  fulfilled  in  our  experiments,  and  Schmid  and 

a  A 

Schon  have  made  similar  suggestions.  ‘  Furthermore, 

4  5 

Lemberger  and  Clarke  have  recently  taken  measurements 
on  the  branch  imbalance  induced  by  tunnel  injection  into 
superconductors  doped  with  paramagnetic  impurities,  and 
have  confirmed  Eq.  (2.16),  again  largely  outside  the 
regime  in  which  the  theoretical  basis  is  most  certain. 
Taken  together  with  these  other  considerations,  our 
experimental  evidence  helps  establish  a  broader  validity 
of  Eq.  (2.16)  for  the  transverse  relaxation  time. 


2 . 4  Heating  Ff  f ects 
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We  have  discussed  how  the  transverse  node  dis¬ 
equilibrium  determines  the  general  form  of  the  I-V  curves. 
The  longitudinal  mode  is  also  excited  in  a  PSC,  however, 
since  the  electric  field  deposits  energy  into  the 
quasiparticles .  In  general  this  produces  a  complicated 
disequilibrium  of  the  quasioarticles  and  the  phonons, 
but  to  a  first  approximation  most  of  its  properties  can 
be  accounted  for  in  terms  of  an  increased  effective 
temperature  of  the  electrons,  i.e.  local  heating. 
Ironically,  heating  effects  tend  to  be  nreater  at  lower 
temperatures,  since  the  critical  current  and  hence  the 
characteristic  dissipation  level  increases.  For  a 
sufficient  amount  of  dissipation,  the  raised  effective 

temperature  at  the  center  of  the  PSC  rises  to  T  and  then 

c 

surpasses  it,  at  which  point  the  dissipative  region 


develops  into  an  expanding  normal  self-heating  hotspot. 

The  properties  of  this  hotspot,  and  more  generally 

the  effects  of  heating  in  long  microbr idnes ,  were 

1  4  6 

discussed  by  Skocpol  et  a_l.  '  in  terms  of  a  simnle 
one-dimensional  macroscopic  steady-state  heat  equation 
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-K  P  *  J^bath’  * 


(2.23) 

where  d  is  the  film  thickness,  t  is  the  heat  transfer 
coefficient  per  unit  area  (the  reciprocal  of  the  thermal 
boundary  resistance) ,  K  is  the  thermal  conductivity  of 
the  metal,  and  JF.  is  the  Joule  heat  ner  unit  volume. 
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Conceptually,  this  equation  balances  heat  production 
against  conduction  along  the  bridge  and  surface  heat 
transfer  across  the  interface  with  the  substrate  or  He 
bath.  The  length  scale  is  characterized  by  a  "thermal 
healing  length" 


n  =  (Kd/ap  (2.24) 

which  represents  the  decay  length  of  the  temperature 
near  a  localized  heat  source. 

We  are  assuming  in  Eq.  (2.23)  that  the  energy 
enters  the  electronic  system,  is  transported  along  the 
bridge  by  quasiparticles  (the  lattice  thermal  conductivity 
can  be  neglected  near  the  critical  temperature  of  a 
superconductor) ,  and  is  transferred  out  of  the  film 
via  phonons.  Since  the  electrons  and  phonons  in  the 
film  can  be  characterized  by  some  temperature  T,  which 
is  larger  than  the  bath  temperature  the  bottle¬ 

neck  in  the  enerav  transfer  is  not  between  the  electrons 
and  the  phonons  in  the  film,  but  between  the  phonons  in 
the  film  and  those  in  the  thermal  bath.  This  is  supported 
experimentally  by  the  fact  that  when  the  He  bath  passes 
through  the  lambda-ooint  transition,  there  is  an  abrupt 

change  in  the  effective  electron  temperature,  because  of 

4  6  * 

the  greatly  enhanced  phonon  transfer  to  superfluid  He.  ' 


47 

*  There  is  some  evidence  in  the  Russian  literature  that 
under  certain  circumstances  (below  the  lambda-point,  and 
for  a  wide  range  of  temperatures  for  extremely  narrow  films) 
the  heating  can  be  described  in  terms  of  electron  "super¬ 
heating"  relative  to  the  lattice. 
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Near  the  critical  temperature,  the  electronic  thermal 

4  8 

conductivity  is  aporoximately  that  in  the  normal  metal  , 
since  only  the  quasiparticles  and  not  the  superconducting 
pairs  can  transport  entropy.  We  can  therefore  write, 
using  standard  kinetic  theory. 


K  =  C  .D 
el 


(2.25) 


where  D  and  C  .  are  the  diffusion  constant  and  heat 
el 

capacity  of  the  normal  metal  at  T  .  The  thermal  conduc¬ 
tance  to  the  bath  can  be  expressed  in  terms  of  the  phonon 
proDerties 


—  =  C  /t 
d  ph'  T  esc 


(2.26) 


where  C  ,  is  the  phonon  heat  caoacity  and  i  is  the 
ph  1  esc 

characteristic  escape  time  of  a  phonon  out  of  the  film. 
The  thermal  healing  length  of  Eq.  (2.24)  can  then  be 
written  as 


n  =  /Dr  C  ./C  .  =  /Dr  ,, 

esc  el  ph  eff 


(2.27) 


where  r  ^  is  the  effective  electronic  cooling  time. 
Furthermore,  by  detailed  balance  and  energy  conservation 
we  have 

Cph/TPh  =  CelAE  (2'28) 

where  tDh  is  the  inelastic  phonon-electon  collision  time 
(i.e.  the  average  time  for  a  given  phonon  to  scatter 
inelastically  off  any  electron) .  Therefore  we  can  write 


r  ef  f  r E  ^  resc'  1 oh^ 


(2.29) 
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This  is  completely  analogous  to  previous  studies  of 

nonequi librium  quasiparticle  recombination  times  at  low 
49 

temperatures.  The  "hot"  electrons  diffuse  a  distance 

vfF_  before  cominq  to  equilibrium  with  the  ohonons  in  the 

h 

film.  In  the  limit  that  the  ohonons  in  the  film  are  main¬ 
tained  in  equilibrium  at  the  bath  temperature,  t  ,,  =  t„  .  * 
1  e  f  f  E 

In  the  more  usual  case  of  "ohonon  traooing",  however, 

>>x  ,  ,  and  a  hot  ohonon  is  much  more  likely  to  inter- 
esc  ph  1 

act  with  the  electrons  to  maintain  the  heated  electron 
distribution  than  to  escape  the  film,  producing  an 
enhancement  by  the  factor  r  /i  ,  in  the  relevant 
quasioarticle  relaxation  time.  In  general,  whenever 
ohonon  trapping  is  dominant,  the  spatial  extent  of  heating 
gives  no  information  on  the  intrinsic  quasiparticle 
relaxation  time,  but  rather  provides  a  measurement  of  the 
phonon  escape  time.  Mote  that  a  similar  argument  does 
not  apply  to  the  case  of  the  transverse  mode,  since 
nonequilibrium  phonons,  being  uncharged,  cannot  create 
or  maintain  a  branch  imbalance. 

*  This  should  be  distinguished  from  the  longitudinal 
relaxation  time  for  transient  processes  jtj  the  super¬ 
conductor,  which  near  T  takes  the  form 

(L)  _  7r  3  kgT  ° 

R  -  7c,  (3)  A 

The  difference  between  this  divergent  time  and  the 
steady-state  relaxation  time,  which  goes  over  continuously 
to  the  normal  time  t  at  Tc,  is  due  to  the  fact  that  in 
the  former  case  the  gap  itself  is  changing  in  time.  That 
these  two  should  be  different  seems  not  to  be  generally 
appreciated. 
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For  the  SST  model  of  the  PSC,  where  the  electric 

field  E  decays  exponentially  over  the  distance  A,  Eq. 

2.23  can  be  solved  explicitly,  yielding  a  temperature 

1  46 

rise  6t  =  CVI  in  the  center  of  the  PSC  ,  where  the 


heating  coefficient  C  is 


2 ( n+A) wa 


(2.30) 


This  is  essentially  what  one  would  expect  for  a  power 

input  VI  distributed  over  a  distance  2A  producing  a 

temperature  rise  extending  out  over  the  larger  range 

2A  +  2 n  .*  A  rough  estimate  of  the  heating  coefficient 

can  be  made  as  follows.  For  small  values  of  the  gap, 

most  of  the  heat  conduction  is  carried  by  quasiparticles, 

and  the  thermal  conductivity  is  approximately  that  of 

the  normal  metal,  which  can  be  determined  from  the 

resistivity  p  using  the  Wiedemann-Franz  law  (K=iT2kgT/3e2  p)  . 

The  thermal  boundary  conductance  a  scales  as  T^  since 

it  is  proportional  to  the  density  of  phonons,  so  that 

n  =  /Kd /ct  should  scale  as  T  Values  of  a(Tc)  of  the 

2 

order  of  2  watts/cm  K  have  been  inferred  from  the  I-V 

characteristic  of  long  tin  microbridges  in  the  self- 

46 

heating  hotspot  regime  ,  where  dissipation  sustains 


the  local  temperature  above  T  . 


*  We  have  assumed  here  that  the  heat  generated  goes  as  JE 
®  exp(-x/A).  It  may  be,  on  the  other  hand,  that  the  heat 
from  the  normal  component  is  generated  as  J  E  <*  exp(-2x/A), 
which  will  produce  a  heated  region  -  2g+A,  and  the  heat 
from  the  condensate  component  is  generated  in  the  core 
region,  producing  a  heating  region  ~2g.  In  any  case,  this 
treatment  is  only  approximate,  and  furthermore ,n  is  gener¬ 
ally  >>A  for  our  data. 
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We  will  assume  that  this  heating  analysis  remains 
essentially  valid  near  the  reduced  T  (H)  in  a  magnetic 
field.  The  oresence  of  the  magnetic  field  will  not 
affect  the  phonon  properties,  and  the  magnitude  of  the 
field  is  not  large  enough  to  produce  any  significant 
magnetoresistive  effects  in  the  relatively  dirty  thin- 
film  samples. 

If  heating  is  imoortant.  of  course,  then  one  is  no 
longer  justified  in  simply  substituting  the  bath  tempera¬ 
ture  into  Eq .  (2.16)  for  the  transverse  relaxation  time. 

If  n  is  much  larger  than  A  ,  however,  then  the  temperature 
rise  at  the  PSC  is  aoproximately  constant  over  the  scale 
of  A  ,  and  one  can  use  the  heating  coefficient  C  to 
determine  the  raised  temperature  for  substitution  into 
Eq.  (2.16).  From  the  above  estimates  of  n,  we  find  that 
n  is  approximately  twice  A  for  H=0 ,  whereas  for  larger 
H  and  lower  7,  n  increases  while  A  decreases  substantially. 
To  a  first  aporoximation ,  therefore,  one  can  take  the 
temperature  rise  to  be  constant  within  the  PSC. 

In  the  SBT  model,  the  isothermal  characteristic  of 
a  PSC  in  the  plateau  region  is 

V  =  Rn(I-Ia)  =  Rn(I-fIc)  (2.31) 

where  f  is  a  fraction  (v%)  which  is  at  least  approximately 
independent  of  T.  "’hen,  since  heating  causes  I  to 
decrease,  it  causes  V  to  increase  if  I  remains  fixed. 


1 1:  R  diverges  .it  T  (H)  ,  heating  causes  R^  t  o  increase, 
thus  furthering  the  same  trend.  Therefore,  since  the 
heat  dissipation  IV  increases  with  increasing  I,  the 
slope  dV/dl  will  increase  as  I  increases.  In  our 
presentation  of  the  I-V  curves,  with  I  along  the  horizontal 
axis,  the  curve  turns  upward  in  a  region  where  it  would 
otherwise  be  straight.  Another  characteristic  indication 
of  heating  is  the  steepening  at  the  onset  of  the  PSC  and 
the  ultimate  appearance  of  hysteresis  between  the  zero- 
voltage  and  dissipative  states  for  current-biased 
measurements.  Again,  this  is  because  1^  in  the  dissipative 
state  can  he  significantly  below  that  corresponding  to 
the  bath  temperature.  These  two  trends  can  be  seen  in 
Fio.  2.7,  where  a  set  of  simulated  I-V  character i st i cs 
are  shown  with  and  without  heating  effects  included. 

For  this  simulated  model  of  heating,  we  take  the  set 
of  T -V'  curves  determined  by  the  equation 


V 


R  I  _ 

«nTr/d-ic) ) 


'  1  '  Tc 


where 


I  ■  I  (l -T/T  ) 
c  CO  c 


3/2 


(2.32) 


(2.31) 


and 


R  -  R  (l-V 
n  o  c 


(2.34) 

These  isothermal  curves  (the  solid  lines  in  Fig.  2.7) 


m 


ar«'  very  steep  near 


and  for  large  1  approach  asymp- 


Kii],  2.7  Simulated  I -V  curves  calculated  front  Kqs . 

(2.32)  -  (2.34)  with  and  without  heating  effects  included, 
using  parameters  I  ,o  -  22  mA,  R  0 .  r>5  and  C  - 
’  '  10*  K  w.  These  curves  are  modeled  after  the  experi¬ 
mental  curves  of  Fig.  2.8.  The  inset  shows  how  the 
minimum  differential  tesi stance  of  the  l -V  curve  varies 

bath  ’ 


w  i  t  It  T 


with  and  without  heating. 


RMRPMi 


■ 
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totically  the  straight  line  V  =  R^tl-Vl^).  The  analytic 

form  of  Fq.  (2.32)  follows  from  a  modified  resist ively 

50 

shunted  junction  (RSJ)  model,  where  one  takes  the 

supercurrent-phase  relation  to  be  piocewise  linear  with 

an  average  value  ovei  the  cycle  of  I  =  (see  Appendix 

II  for  details).  A  current-phase  relation  of  this 

general  form  was  found  for  long  bridges  by  DC  measure- 
50 

ments.  Many  other  reasonable  cur rent -phase  relations 

yield  qualitatively  similar  i-v  curves.  The  above 

temperature  dependence  of  the  critical  current  is  the 

standard  mean-field  behavior,  and  is  the  correct  form 

for  a  long  microbridge.  The  given  behavior  for  is 

in  accord  with  a  divergence  of  A  *  in  the  relaxation 

time  within  the  Sut  model.  The  values  of  the  coefficients 

R  and  I  in  Fig.  2.7  were  fixed  bv  modeling  this 
o  co  '  J 

system  after  a  particular  set  of  experimental  data. 

Bridge  ft  1 5 A  for  H=0,  presented  for  comparison  in  Fig.  2.8. 

Ileatinn  is  introduced  into  the  simulation  by  assuming 
that  the  temperature  which  enters  into  determining  1  and 
Rn  is  raised  above  the  bath  temperature  by  the  quantity 
iST  =  CVI,  where  the  single  value  C  =  2  x  10^  K/W,  consist¬ 
ent  with  experimental  estimates  for  Bridge  015a,  is  used 
for  all  the  curves  in  Fig.  2.7.  The  dashed  lines  repre¬ 
sent  the  set  of  curves  with  heating  included,  correspond¬ 
ing  to  the  same  set  of  bath  temperatures  as  the  isothermal 


L 
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Fig.  2.8  Experimental  I-V  curves  for  Sn  bridge  #15A 
for  H  =  0,  for  a  range  of  temperatures  near  T  .  The 
temperatures,  together  with  the  corresponding  measured 
values  of  R^ ,  are:  curve  a  -  T=3.389  K,  R  =1.83  fi; 

B  3.867  K,  2.05  C  3.849  K,  2.49  Q;  D  3.828  K,  2.4  Q 
E  3.806  K,  2.55  Q;  F  3.782  K,  2.55  Q.  These  curves 
are  to  be  compared  with  the  corresponding  simulated 
I-V  curves  in  Fig.  2.7. 
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curves.  Mote  that  the  heating-induced  distortion 

increases  quite  strongly  the  farther  one  goes  from  T  . 

This  follows  from  the  fact  that  both  V  and  I  scale  with 

lc  in  the  isothermal  curve,  so  that  the  heat  dissipation 
2 

VI  goes  as  I  for  comparable  positions  along  the  curve. 

If  heating  were  negligible,  the  best  measurement  of 
the  resistance  R  could  be  made  at  the  largest  current 
before  the  second  PSC  entered,  since  the  curve  is  asymp¬ 
totic  to  the  straight  line  with  slope  R  .  In  the  presence 

n 

of  heating,  with  the  tendency  toward  increasing  values 
of  dv/dl  for  large  currents,  we  take  the  slope  at  the 
largest  current  before  heating  becomes  evident,  i.e.  we 
take  the  minimum  value  of  dV/dl.  /'s  one  gets  farther  from 
Tc  and  heating  becomes  more  significant,  however,  this 
becomes  a  less  and  less  reliable  method  of  measuring  the 
appropriate  slope.  Since  the  intrinsic  isothermal  curve 
is  concave  downward  and  heating  tends  to  make  this  concave 
upward  for  large  currents,  there  can  often  be  a  deceptive¬ 
ly  straight  region  of  the  curve  around  dV/dl^^,  which 
however  is  significantly  greater  in  slope  than  R  . 

This  is  evident  from  the  inset  of  Fig.  2.7,  which  shows 

the  major  effect  of  heating  upon  dv/dl  .  .  In  the  case 

mm 

of  the  simulation  here,  the  intrinsic  resistance 

(the  solid  curve)  is  decreasing  as  one  moves  away  from 
Tc,  while  beyond  a  certain  Point  the  minimum  value  of 
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dV/dl  with  heating  is  increasing.  The  sane  thing  occurs 
in  our  experimental  data,  and  when  heating  is  very 
significant,  this  can  obscure  the  divergence  in  the 
relaxation  time  as  T  -*•  Tc(H). 

The  hysteresis  in  the  critical  currents  is  also 
quite  similar  between  Figs.  2.7  and  2.8.  In  part,  this 
may  be  because  the  value  of  C  was  chosen  so  as  to  match 
the  onsets  of  hysteresis  in  the  two  cases.  The  develop¬ 
ment  of  hysteresis  for  lower  temperatures  is  also  in  good 
agreement,  however,  implying  that  this  hysteresis  may  be 
largely  attributed  to  heating.  On  the  other  hand,  a  few 
reservations  about  hysteresis  may  be  in  order.  For  one 
thing,  constant  current-biasing  in  a  noisy  environment 
may  lead  to  premature  switching  to  the  dissipative  state 
or  back  again  (see  e.g.  curve  F  in  Fig.  2.8),  leading  to 
a  possible  underestimate  of  the  actual  hysteresis. 
Secondly,  the  apparent  onset  of  hysteresis  in  our 
samples"^  tended  to  occur  at  about  30  to  50  mK  from  T^ 

(in  zero  field),  whereas  measurements  on  similar  samoles 
of  Skocpol,51  taken  in  an  rf-shielded  room,  indicated  an 
onset  considerably  closer,  typically  5  -  10  mK  from  T  . 

It  may  be  that  some  hysteresis  in  our  data  is  hidden  by 
rf  noise-rounding,  in  spite  of  the  rf  filters  used. 

(Note,  in  this  regard,  the  rounded  onsets  of  the  non- 
hysteretic  curves  in  Fig.  2.8.)  In  that  case,  an  alter¬ 
native  explanation  of  hysteresis  close  to  T^  might  be 


required  (a  possible  source  of  intrinsic  hysteresis  is 
discussed  in  Chapter  IV) . 

In  practice,  we  do  not  know  whether  any  simple  set 
of  formulas  such  as  Eqs.  (2.32)  -  (2.34)  can  actually 
generate  our  experimental  I-V  curves.  We  were  able, 
however,  to  reconstruct  isothermal  curves  from  our  raw 
data,  at  least  when  the  heating-induced  distortion  was 
not  too  great  (in  particular,  necessary  information  is 
lost  in  hysteresis) .  The  procedure  begins  by  considering 
a  set  of  I-V  curves  for  constant  H  and  a  series  of  bath 


temperatures  near  Tc<H) ,  as  in  Fig.  2.9a.  Given  a  value 
for  the  heating  coefficient  C,  the  raised  local  tempera¬ 
ture  Tioca;L  =  Tbath  +  can  b®  comPute<3  at  various 

points  along  each  curve  of  the  set.  One  can  then  connect 
the  points  corresponding  to  the  same  value  of  T^oca^» 
interpolating  where  necessary,  to  obtain  a  family  of 
isothermal  curves.  The  curves  in  Fig.  2.9a  were  pro¬ 
cessed  in  this  way  to  yield  those  of  Fig.  2.9b,  using 
for  simplicity  a  single  value  of  C  for  all  the  curves. 

The  magnitude  of  C  was  chosen  so  as  to  bring  the  greatest 
reduction  in  the  amount  of  heating-induced  curvature  for 


the  whole  set  of  curves. 


The  values  of  C  obtained  in  this  way  for  several 
values  of  H  are  plotted  in  Fig.  2.10  as  a  function  of 


Tbath”  T*ie  h°ri-zontal  bars  represent  the  range  of 


heating"  relative  to  the  lattice. 
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Fig.  2.9  I-V  curves  for  bridge  *15A  for  H  =  2000  Oe 
(  ||  J)  before  and  after  processing  to  remove  heating- 
induced  distortion,  a)  Raw  I-V  curves  for  the  first 
PSC.  b)  Reconstructed  isothermal  curves  corresponding 
to  the  bath  temperatures  in  (a) .  See  text  for  clarifi¬ 
cation.  Curves  for  bath  temperatures  N  and  O  are  not 
reconstructed  because  this  would  require  extrapolation 
byond  the  measured  data.  Bath  temperatures:  A,  3.309  K: 


B  ,  3.284; 

c, 

3.250; 

n,  3.212; 

F. ,  3.149; 

F,  3.086; 

G,  3.037; 

H, 

2.977; 

I,  2.939; 

J,  2.852; 

K,  2.775; 

L,  2.670; 

M. 

2.622; 

M,  2.546; 

0,  2.501. 

appreciated . 
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temperatures  over  which  the  fit  was  made,  and  the 

vertical  bars  are  the  estimated  uncertainties.  These 

values  are  compared  to  independent  estimates  of  C 

obtained  from  F.o .  (2.30)  and  experimental  parameters. 

The  value  of  a(Tc)  is  known  only  roughly,  from  data  by 
46 

Skocpol  in  the  self-heating  hotspot  regime.  In  Fig. 

2 

2.10,  curves  are  olotted  using  (m  )  =2  and  3  W/cm  K 

c 

for  T  >  T^  (where  heat  is  transferred  primarily  to  the 

substrate  rather  than  directly  to  the  normal  He  bath) 

2 

and  8  W/cm  K  below  the  lambda  point,  where  heat  transfer 
out  of  the  film  suddenly  increases  due  to  the  superfluid 
He.  Since  we  are  comparing  with  data  in  the  oresence  of 
magnetic  fields  near  Tc(H),  the  lower  temperatures 
correspond  to  higher  fields,  where  and  C  becomes 

effectively  a  function  only  of  T.  Above  about  3.4  K, 
however,  A(H)  is  sufficiently  large  compared  to  n  to 
cause  a  significant  dependence  of  C  on  magnetic  field. 

The  computed  curves  there  are  not  unique,  but  use 
typical  values  of  A ( H ) . 

Overall,  considering  the  roughness  of  the  approach, 

the  agreement  between  these  two  sets  of  estimates  for  C 

seems  to  be  quite  satisfactory,  providing  support  for  the 

1  4  6 

simple  heating  model  of  ^kocpol  et  a_l.  '  The  evidence 
seems  to  indicate  that  the  "electron  heating"  is  dominated 
by  phonon-trapping  effects  (at  least  for  T  N  T^),  so  that 


-  i- - — ~  i-c^xun  -  c.\\-rn,  auu  cut:  ;icau 

from  the  condensate  component  is  generated  in  the  core 
region,  producing  a  heating  region  ~2r\.  In  any  case,  this 
treatment  is  only  approximate,  and  furthermore , n  is  gener¬ 
ally  >>A  for  our  data. 
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Fig.  2.10  Heating  coefficient  C  versus  bath  temperature 
for  bridge  $15A.  The  experimental  points,  marked  by  (-4—*, 
represent  constant  values  chosen  to  minimice  apparent 
heatinci-induced  distortion  for  a  set  of  1-V  curves 
over  a  range  of  temperatures,  with  the  indicated  un¬ 
certainty.  The  lines  represent  semi-empirical  predictions 

of  Fa.  (2.30)  ,  usina  tvpical  values  of  \(m  )  and  tvnical 

c  *  * 

variation  of  A(H) .  At  the  onset  of  superfluidity  in  the 
He  bath  (2.18  K)  ,  the  appropriate  value  of  Mm  )  changes 
abruptly,  due  to  the  additional  heat- t ransf er  ability 
of  the  superfluid. 
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we  are  probably  unable  to  obtain  anv  information  on  the 
effect  of  the  magnetic  field  on  the  intrinsic  longitudinal 
relaxation  time  from  the  variation  of  the  thermal  healing 
length  with  field.  In  contrast  to  the  transverse  mode 
case,  however,  we  would  not  expect  the  field  to  have  a 
dramatic  effect,  since  near  T  (H)  the  intrinsic  steady- 
state  relaxation  time  should  be  approximately  the 
electron-phonon  collision  time  t  (T) ,  which  is  essen¬ 
tially  independent  of  magnetic  field. 

Using  the  curves  corrected  for  heating  effects,  we 

(T' 

have  found  that  the  variation  of  i  for  constant  H 

K 

can  be  described  fully  by  a  divergence  as  T^fH)  is 

aporoached;  the  increase  of  dV/dl  (T)  in  the  raw  data 

mm 

far  from  T  (H)  is  due  to  heating.  This  can  be  seen, 
for  example,  in  the  data  in  Fig.  2.5  for  H  =  1200  Oe. 

Using  the  data  corrected  for  heating,  we  are  able  to 
provide  a  better  fit  to  tire  Schmid  and  Schon  formula 


F.q .  (2.16). 


Heating  is  a  fundamental  problem  intrinsic  to 
uniform  one-dimensional  superconducting  filaments. 

An  estimate  of  the  degree  ot  heating-induced  distortion 
in  zero  field  can  be  found  in  the  dimensionless  parameter 


Y  = 


■ST 

T  -T 
c 


I  R 
c  n 


2  (n  +  .\)w  r(T  -T) 
c 


K 


V3  _A_ 

T  -T  n+A 
c 


I 


(2.35) 
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which  represents  the  relative  temperature  rise  at  the 

2 

characteristic  dissipation  level  I  R  .  Here  o  is  the 

c  n 

normal-state  conductivity  and  J  is  the  critical  current 

'  c 

density,  which  for  H=0  in  the  dirty  limit  near  T  takes 

c 

the  form  (see  ref.  52  and  Fqs .  (3.28)  and  (3.29)  in 
Chapter  III) 


,  _  a  2  h 

c  nsevc  2eiQ  ^  TTcTtT 


=  2.837 


kB(VT) 


(2.36) 


eC(T) 

Then,  using  the  Widemann-Franz  law  and  estimating 
A/(n+A)  -  1/3,  we  obtain  the  very  simnle  expression 


v  *  Tef f  ^T1 


(2.37) 


where  Teff  is  the  effective  cooling  time  discussed 


earlier . 


This  heatina  parameter  is  independent  of  the 
sample  width,  the  electron  mean  free  path,  and  even  the 
thickness  of  the  film  (provided  it  is  small  enough  to 
optimize  phonon  escape).  Furthermore,  since  Te££  -  tp, 
the  coefficient  Te££A(0)/fi  is  very  large,  about  1000 
for  Sn.  Even  if  T  is  within  lfl  of  Tc>  which  is  tyoical 
for  our  data,  y  a,  0.1,  corresponding  to  a  significant 
amount  of  heating  and  consequent  distortion.  For  Al, 
with  a  value  of  r  two  orders  of  magnitude  larger  than 
that  for  Fn,  heating  should  be  an  even  more  serious 
problem  at  comparable  reduced  temperatures.  In  the  clean 
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limit  <"?.  ,  i  is  reduced  by  a  factor  /i  ,  hut  it  is 

difficult  to  obtain  thin-film  samples  clean  enough  for 

this  to  apply  (although  it  may  in  clean  whisker  crystals) 

The  presence  of  a  parallel  magnetic  field  will  reduce 

heating  somewhat  because  of  the  lower  character ist ic 

resistance  R  ■>  \,  but  heating  will  still  be  quite 

significant.  The  method  we  have  used  to  reduce  the 

effects  ot  heating  to  manageable  proportions  is  to 

produce  a  localised  reduction  in  the  critical  current 

1  ,  which  will  reduce  the  heating  parameter  i  by  a 

factor  ot  l  We  will  discuss  the  use  ot  these  local  ir.ei 

c 

weak  spot s  further  in  the  next  section. 
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2.5  Localized  Weak  Spots 


The  critical  currents  of  successive  PSC's  in  a 

bridge  which  is  nominally  uniform  can  vary  somewhat,  due 

to  slight  variations  in  or  in  cross-section  along  the 
1  21 

bridge,  '  or  to  the  influence  of  one  PSC  on  others 
1,53  54 

nearby.  '  Fven  so,  it  is  often  difficult  in  such 
uniform  bridges  to  make  certain  that  in  following  the 

PSC  with  the  smallest  I  through  large  variations  in  T 

c 

and  H,  one  is  actually  following  the  same  single  unit. 

Furthermore,  near  T  (H)  the  second  PSC  often  enters 

c 

before  the  first  PSC  has  reached  the  plateau  region  in 
its  I-V  curve,  making  an  accurate  determination  of  R 
difficult.  For  these  reasons,  a  sample  having  a  single 
deliberate  "weak  soot"  associated  with  a  locally  reduced 


value  of  I  is  highly  advantageous.  In  addition,  since 

2 

the  heat  dissipated  at  the  PSC  is  of  order  1  ‘R  ,  de- 
r  c  n 

creasing  Tc  reduces  the  heating-induced  distortion  of 
the  I-V  curve  substantially.  Since  we  are  interested 
in  probing  the  nonequilibrium  properties  of  the  uniform 


region  outside  the  weak  spot,  however,  the  length  of  this 
weak  soot  should  be  small  compared  with  A  and  the  resist¬ 
ance  it  contributes  should  be  small  compared  with  R  .  * 

n 


*  Tf  the  weak  spot  is  made  by  a  const rict ion ,  it  will 
add  some  additional  resistance  even  if  it  is  short.  The 
requirement  that  this  additional  resistance  be  small 
relative  to  R  puts  a  practical  limit  on  the  degree  to 
which  the  critical  current  may  be  reduced  by  the  con¬ 
striction  weakening  illustrated  in  Fig.  2.2. 


. 
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For  a  uniform  bridge  without  a  weak  spot,  the  I-V 

curves  are  relatively  independent  of  the  orientation  of 

the  magnetic  field  relative  to  the  direction  of  the 

current  density  J,  as  long  as  it  is  within  the  plane  of 

the  film.  The  critical  current  I  of  a  weak  soot, 

c 

however,  is  depressed  far  more  in  a  large  magnetic 
field  H  ||  J  (and  perpendicular  to  the  cut)  than  it  is 
for  the  same  field  with  H  [  J  (and  parallel  to  the  cut) . 
For  H  [  J,  the  critical  current  of  the  weak  spot  seems  to 
behave  similarly  as  a  function  of  H  and  T  to  the  corre¬ 
sponding  Ic's  of  the  other  PSC’s.  The  relative  "weakness 
of  the  weak  spot  is  then  presumably  due  to  a  simple 
decrease  in  the  local  cross-sectional  area.  For  H  ||  J, 
on  the  other  hand,  not  only  is  I  of  the  weak  spot 
smaller  than  that  of  the  other  PSC's,  but  it  also  goes 
to  zero  at  a  lower  temperature.  This  seems  to  be  par¬ 
ticularly  true  for  the  "groove"  type  of  weak  spot  (Fig. 
2.2c)  where  the  cross  section  at  the  weak  spot  is  not 
really  reduced  at  all.  We  note  that  in  this  orientation, 
the  field  has  a  component  perpendicular  to  the  surface  at 
the  weak  spot.  Since  perpendicular  critical  fields  are 
substantially  smaller  than  parallel  critical  fields,  a 
depression  of  the  critical  temperature  T^fH)  at  the 
localized  weak  spot  is  not  unreasonable.  We  will  discuss 
some  further  aspects  of  the  critical  currents  of  these 
magnetically  weakened  regions  in  Section  5.2. 
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Due  to  the  localized  depression  of  Tc(H),  there  is  ■ 

a  significant  range  of  temperatures  for  large  H  ||j  where 
Ic  of  the  weak  spot  has  already  gone  to  zero,  whereas 
that  of  the  second  PSC  has  not.  We  are  therefore  able 

(T) 

to  extract  tr  from  Rn  of  the  first  PSC  over  this 

entire  range,  and  we  find  that  nothing  unusual  happens 

when  Ic  of  the  weak  spot  (and  presumably  also  the  local 

* 

value  of  the  gap  A)  goes  to  zero.  The  divergence  of 
(T) 

tr  seems  to  occur  at  the  undepressed  value  of  T  (H) 

as  determined  by  the  higher  PSC's.  This  is  shown  in 

-4 

Fig.  2.11a,  where  Rn  is  plotted  versus  T  to  show  that 

_v 

the  divergence  of  Rn  as  [l-T/Tc(H)]  occurs  at  a  tempera¬ 
ture  150  mK  above  the  temperature  where  I  went  to  zero. 

This  can  be  understood  by  noting  that  the  magnetic 
depression  of  I  of  the  weak  spot  is  a  highly  localized 
phenomenon.  Outside  this  region,  yet  still  within  the 
range  of  the  quasiparticle  diffusion  length,  the  appro¬ 
priate  critical  temperature  becomes  that  corresponding  to 
the  uniform  film.  The  same  is  true,  although  to  a  much 

lesser  extent,  for  a  uniform  bridge  in  zero  magnetic  'Ss. 

field  (see  Fig.  2.11b).  Even  for  the  uniform  bridge, 

*  When  the  center  region  has  gone  normal,  what  we  have 

is  more  correctly  an  SNS  junction  than  a  phase-slip 

18 

center.  Some  recent  measurements  on  large-area  SNS 
junctions  confirm  that  the  excess  resistance  contributed 
by  the  superconductor  corresponds  near  Tc  to  the  length 
A  =  /DTq  on  either  side  of  the  normal  region. 


. - . .  i  ii . ^ ‘  ~ 
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there  are  significant  local  variations  of  T  ,  typically 

1  19 

about  30  mK  for  Sn ,  '  which  cause  some  of  the  variation 

in  I  of  the  PSC 1 s .  Since  I  is  determined  by  local 
c  c  1 

characteristics  over  a  region  of  the  order  of  the 
coherence  length,  the  "nucleation  center"  of  a  PSC 
may  be  expected  to  have  a  somewhat  smaller  value  of  Tc 
than  that  of  some  of  the  surrounding  regions  within  a 
distance  of  the  quasiparticle  diffusion  length  A. 

The  data  presented  in  Figs.  2.5  and  2.6  come  from 
samples  with  magnetically  weakened  spots.  In  samples  not 
containing  such  a  feature,  the  divergence  of  Rn  near  T  (H) 
is  often  largely  hidden  by  the  presence  of  the  other  PSC's 
and  by  heating.  T’he  nominally  uniform  sample  in  Fig. 

O 

2.11b  was  less  than  200  A  thick  and  was  well-adhered  to 
the  sapphire  substrate,  both  features  which  should 
optimize  the  conduction  of  heat  out  of  the  film.  Yet 
data  indicating  the  divergent  behavior  covered  a  range 
of  only  about  50  mK ,  before  heating-induced  distortion 
started  becoming  important.  This  distortion  can  be 
approximately  removed  by  the  methods  described  earlier 
to  extend  this  range  somewhat,  but  it  remains  quite 
limited.  Compare  this  to  Fig.  2.11a,  where  the  corre¬ 
sponding  data  for  a  bridge  with  a  magnetically  weakened 
spot  (not  particularly  optimized  for  heat  removal)  shows 
the  divergent  behavior  over  a  range  of  300  mK.  It  is 
perhaps  ironic,  but  nontheless  true,  that  it  is  necessary 
to  use  a  weak  soot  to  probe  the  properties  of  a  uniform  film. 
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2 . 6  Conclusions 

In  summary,  we  have  examined  the  I-V  curves  of  long 

tin  microbridges  in  parallel  magnetic  fields  and  have 

successfully  interpreted  them  within  the  phase-slip 

center  model  of  Skocpol,  Beasley,  and  Tinkham.^  The 

differential  resistance  in  the  plateau  region  of  the  PSC, 

after  correcting  for  dissipative  effects  using  a  simple 

heating  model,  is  set  equal  to  the  normal  resistance  Rn 

of  a  region  of  the  bridge  with  length  2/Dt  ,  where  D 

K 

is  the  quasiparticle  diffusion  constant  and  t  is  the 

K 

appropriate  relaxation  time.  The  temperature  and  field 

dependence  of  Tn  thus  inferred  from  our  data  agrees  well 
K 

with  that  of  the  transverse  mode  relaxation  time  in  the 

••  12 

presence  of  pair-breaking  derived  by  Schmid  and  Schon 
[Eq.  (2.15)].  The  values  of  t  (T  )  which  orovide  the 
best  fit  with  the  data  fall  between  1  and  5  *  10  10  sec, 
and  much  of  this  scatter  can  be  attributed  to  an  un¬ 
certain  knowledge  of  the  experimental  factors  relating 
the  measured  resistance  to  the  corresponding  time. 

In  any  case,  these  figures  are  generally  consistent 
with  other  estimates  obtained  for  this  constant  in  tin, 
both  experimentally  and  theoretically.  The  elucidation 
of  the  role  of  heating  in  modifying  the  I-V  characteristics 
has  been  very  important  in  gaining  a  more  complete  under¬ 
standing  of  them.  By  using  a  semi-empirical  heat  transfer 


1 
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model  with  local  heating,  corresponding  to  the  regime  of  t 

phonon  trapping,  we  have  been  able  to  account  for  most 
of  these  heating  induced  distortions.  Finally,  we  re¬ 
emphasize  the  importance  of  the  deliberately  inserted 
weak  spots,  which  caused  significant  depression  of  I 
for  large  H  ||J,  in  allowing  us  to  extend  substantially 
the  temperature  range  near  T^(H)  where  data  relatively 
undistorted  by  heating  could  be  obtained. 


CHAPTFR  THREE:  CHARGE  IMBALANCE  DYNAMICS  IN 


NONEQUI LIBRIUM  SUPERCONDUCTORS 

3 . 1  Introduction 

We  now  step  back  from  phase-slip  centers  to  examine 
more  generally  the  dynamics  of  charge  imbalance  in  non¬ 
equilibrium  superconductors.  Using  a  generalized  two- 
fluid  approach  to  BCS  superconductors  near  T  ,  we  derive 
in  Section  3.2  a  charge  imbalance  wave  equation,  which 
provides  a  unified  picture  capable  of  dealing  with  both 
steady-state  branch  imbalance  relaxation  and  ligh-fre- 
quency  collective  mode  propagation.  We  show  in  Section 
3.3  that  the  associated  equations  for  the  nonequilibrium 
currents  and  potentials  are  isomorphic  to  equations 
which  describe  electrical  signals  on  a  simple  classical 
transmission  line,  and  oroceed  to  use  analogies  with 
this  well-known  case  to  illustrate  directly  and  graphic¬ 
ally  the  physical  significance  of  the  nonequilibrium 
processes.  In  Section  3.4  we  demonstrate  the  facility 
of  this  transmission-line  picture,  as  well  as  some  of 
its  limitations,  by  discussing  a  number  of  previous 
experiments  involving  nonequilibrium  potentials  in 
superconductors,  and  make  some  predictions  for  experi¬ 
ments  yet  to  be  performed.  We  conclude  Chapter  III  by 
summarizing  those  elements  of  the  analysis  which  will  be 
useful  in  developing  our  model  for  the  phase-slip  center 
in  Chapter  IV. 
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3.2  Derivation  of  Charge  Inbalance  Waves 

Our  approach  to  charge  imbalance  dynamics  follows 

in  part  the  treatments  of  Waldram3*3  and  Pethick  and 
24 

Smith.  While  many  of  our  results  have  been  obtained 
previously,  our  heuristic  derivation  of  the  charge 
imbalance  wave  equation  clarifies  the  physical  origin 
of  the  waves  and  brings  out  certain  aspects  not  previous¬ 
ly  emphasized.  The  quantitative  results  are  strictly 
limited  to  a  regime  very  close  to  Tc  in  the  absence  of 
pair  breaking,  but  we  expect  that  the  qualitative  picture 
has  a  wider  validity. 

The  electronic  charge  density  in  a  BCS  superconductor 
is 

Qe  =  2e^[ukfk  +  vk(1-fk)]  (3.1) 

k 

where  the  sum  is  over  all  electronic  k-states  (with 
the  factor  of  two  for  the  spins)  and  f^  is  the  quasi¬ 
particle  distribution.  The  ion  cores,  of  course, 
contribute  an  eoual  and  opposite  charge  to  maintain 
overall  charge  neutrality.  The  factors  u£  and  v^  can 
be  interpreted  as  the  probabilities  that  a  given  quasi¬ 
particle  state  represents  an  electron  and  a  hole  resDec- 
tively,  and  are  given  by 


where 


E 


k 


1  -  Vs  ■  *<i  +  W 


(3.2) 


1  IIMM— 


1 


N 


(3.3) 
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Here  is  the  quasiparticle  excitation  enerqy ,  positive 

on  both  the  electron-like  (f.^0)  and  the  hole-like  (£.,  <0) 

branches  of  the  excitation  spectrum,  A  is  the  3CS  gap 

parameter,  aid  Cj.  is  the  electron  kinetic  energy  referred 

to  the  chemical  Dotential  c s  of  the  superconducting 

condensate  (which  may  be  different  from  the  Fermi  energy 

tp^O).  For  simplicity,  we  have  dropoed  the  terms  involving 

5  5 

the  superfluid  velocity  ,  even  thouqh  this  approximation 
may  be  valid  only  in  situations  where  the  supercurrent 
is  small  compared  to  the  critical  current. 

We  can  express  Qg  as  the  sum  of  the  two  components 

0*  =  2e£f  (uk2  -  vk2)  =  2^q.f  (3.4) 

k  k 

and 

Qc  =  2eSvk2  •  (3.5) 

k 

Q*  is  the  net  charge  density  in  the  quasiparticle 
system,  commonly  called  the  branch  imbalance  or  charge 

It 

imbalance.  The  quantitv 

qk  =  e{uk2  "  vk2)  =  e  r’k/rk 
in  Eq.  (3.4)  can  be  regarded  as  an  "effective  charge" 

*  Notation  and  units  vary;  our  0*  has  units  of  charge 

density,  that  of  refs.  8  and  14,  number  density  (=o*/e) , 

and  that  of  ref.  10,  energy  (=0*/2N ( 0) e) .  The  quantities 

SO3-'6'  and  50  in  ref.  24  are  related  to  our  Q *  by 
n  n 

.  4k  T 

6Q„'  =  =  °*/e  • 

n  ti  A  n 
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for  a  given  quasiparticle.  The  remaining  part  Qc  of 

the  total  electronic  charge  density  can  be  attributed 

to  the  condensate,  but  should  not  be  confused  with  the 

2 

superfluid  density  ns  <*  A  . 

When  the  electronic  system  is  in  equilibrium, 

£s=0,  and  f  k  is  simply  the  Fermi  function  f°,  given  by 

f°(ek)  =  [1  +  exo(/c£  +AVk?T)]'1  (3.7) 

where  kg  is  Boltzmann's  constant.  (We  will  be  considering 

only  cases  where  the  gap  A  and  the  temperature  T  are 

constant,  so  that  we  can  safely  suppress  the  functional 

dependence  on  them.)  Since  qk  is  an  odd  function  of 

Ck  and  f°(ek)  an  even  one,  the  charge  imbalance  Q*  =  0 

and  all  the  electronic  charge  resides  in  the  condensate. 

For  a  general  disequilibrium,  both  es  and  Q*  may  be 

nonzero,  but  overall  charge  neutrality  should  still  be 
5  6 

maintained.  Therefore  we  have  (see  Fig.  3.1a) 


or  qk: 


Q* 


2£  <qkfk  +  fkqk> 

K 


(3.9) 
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Fig.  3.1  The  effects  of  a  shift  in  the  condensate 

chemical  potential  eg  on  the  condensate  and  on  the 

quasiparticles.  The  dashed  lines  indicate  equilibrium, 

the  solid  lines  disequilibrium,  and  the  dotted  line 

2 

m  (c)  local  equilibrium.  a)  gives  the  occupation 

of  electrons  in  the  condensate.  b)  A  shift  in  e 

s 

shifts  the  excitation  spectrum,  but  in  the  absence  of 

inelastic  collisions,  a  given  quasiparticle  remains 

at  the  same  value  of  but  changes  its  excitation  energy 

and  its  effective  charge  .  c)  If  the  nonequilibrium 

quasiparticles  can  be  characterized  by  a  quasi-equilibrium 

distribution  f°(e.  -e  )  around  a  "normal"  chemical 

k  n 

potential  e  ,  this  is  actually  closer  to  the  global 
equilibrium  distribution  f°(ek)  than  to  the  local 
equilibrium  distribution  f°(e.-e  ).  The  extent  of  the 

K  S 

disequilibrium  in  the  figures  here  is  exaggerated  for 
clarity. 


R5 


A  change  in  the  distribution  function  f  ^  can  occur  either 
by  direct  injection  of  quasiparticles,  quas ioarticle 
diffusion  from  adjacent  regions,  or  collisions.  Thus 
inelastic  electron-phonon  scattering  between  states 
of  different  g.  will  change  0*.  In  addition,  if  s 

K  S 

changes,  the  effective  charge  q^  will  change  for  all 
the  quasiparticles  (see  Fig.  3.1b),  which  has  an  even 
greater  effect  on  Q*. 

In  a  steady-state  injection  experiment,  such  as 
those  by  Clarke  and  colleagues,  ^  t  ^  and  lienee  the  q^'s 
are  constant.  The  resulting  steady-state  charge  imbalance 
will  then  be  pronort iona 1  to  the  time  \ ^  which  character¬ 
izes  relaxation  by  changes  in  the  distribution  function 

alone,  which  has  been  shown  to  take  the  following  form 
12 

near  T  : 


4knT 


0* 


n  A  F. 


'  c 


(3.10) 


where  is  the  inelastic  electron-phonon  collision 

time.  '’’he  k  T/A  enhancement  factor  reflects  the  fact 

that  a  transition  between  states  of  different  o,  can 

k 

occur  only  if  either  the  initial  or  the  final  state 
lies  within  approximately  A  of  the  gap  edge. 

In  a  pulsed  experiment,  however,  the  change  in  the 
effective  charge  q.  plays  an  essential  role.  Since 


is?  ‘W  ■  -e‘v'2'V 


(3.11) 


F.q.  (3.9)  becomes 
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Q*  =  2j]qkfk  -  2eeg  ^  fkA2/ER3  ,  (3.12) 

k  k 

so  that  using  Eq.  (3.8)  we  obtain 

n*  -  uE  Vki/[i-  jirrr?  VX3'  •  <3-131 

Since  the  denominator  in  Eq.  (3.13)  is  oositive  (for 
fk  <  ^)  the  response  of  the  superfluid  in  maintaining 
electroneutrality  acts  to  enhance  the  change  in  Q*  due 
to  the  change  in  the  distribution  function  alone. 

★ 

In  order  to  evaluate  this  enhancement  factor , 
we  write  the  nonequilibrium  distribution  function  in 
the  form 


fk  "  f°(e 
where  f°(ek  - 
function  [Eq. 
potential  e  , 


k  ‘  cs>  +  6fk 


4k 


+  6f‘ 


(3.14) 


eg)  represents  the  equilibrium  distribution 
(3.7)]  centered  about  the  shifted  chemical 
and  the  form  on  the  right  is  appropriate 


for  low  energies  near  7  .  We  take  the  small  deviation 

c 

from  equilibrium  6fk  to  be  odd  in  £k  so  as  to  produce  a 

charge  imbalance.  A  distribution  of  this  form  implies 

that  the  gap  A  has  its  equilibrium  value,  so  that  in  the 

12 

terminology  of  Schmid  and  Schon,  the  disequilibrium  is 
purely  transverse,  the  longitudinal  mode  (which  deals 
with  changes  in  the  magnitude  of  the  gap)  not  being 
excited.  If  we  substitute  the  above  distribution  into 
the  denominator  of  Eq.  (3.13),  and  note  that  6fk  cancels 


*  Although  the  present  treatment  is  restricted  to  the 
regime  close  to  T  ,  the  denominator  in  Eq.  (3.13)  can  be 
evaluated  more  generally,  in  terms  of  the  integral  Z(T) 
defined  in  ref.  14,  to  yield  an  enhancement  factor  of 
1/fl-Z (T) ) . 
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out  since  it  is  odd  in  £k  we  obtain 


/*  r 

,  I  a  i  1  k  .  A  ‘  .  . ,  n  A  .  _  n  A  , , 

1  ~  J.J k (y  "  (  4kB-r  "  4kBT  •  (3' 


Therefore,  Eq .  (3.13)  becomes 

•Ik  r 

°*  *  -TT  '2E  Vx1 


(3.16) 


so  that  the  total  chance  in  0*  is  enhanced  over  the 
explicit  change  (due  to  the  changes  in  f^)  by  the  factor 
4kBT/nA. 

In  general  the  factor  in  brackets  in  Eq .  (3.16) 

will  be  the  sum  of  terms  for  the  various  mechanisms  by 

which  the  quasiparticle  distribution  function  can  chance, 

e.g.  inelastic  scattering,  quasiparticle  diffusion  and 

tunnel  injection.  For  the  simplest  case ,  relaxation  of  a 

homogeneous  charce  imbalance  by  inelastic  scattering, 

we  can  use  a  relaxation  time  approximation  with  the 

single  time  (although  the  actual  relaxation  time 

4  0 

is  expected  to  be  energy  dependent  ) : 

■  -  vlfk '  f°"k '  •  (3-17’ 

Note  that  recharacterizes  the  rate  at  which  fj.  relaxes 

toward  the  "local  equilibrium"  distribution  f°(i^  -  <  ) , 

rather  than  towards  the  "global  equilibrium"  distribution 

f°(ck) .  Substituting  into  Eq.  (3.16),  we  find  that 

f°(r.  -  r  )  contributes  zero  charge  imbalance,  and 

K  S 


4kBT 


1-Q*/t0J  -  -OVt, 


(3.18) 


Thus  a  t imc -resol ved  observation  of  charge  imbalance 
decay  following  pulse  injection  of  quasiparticles  should 
yield  a  relaxation  time  r  ,  rather  than  the  value  r  . 

1.  Q  * 

which  is  appropriate  for  steady-state  experiments  where 
t'  is  held  constant.  To  our  knowledqe  this  prediction 
has  not  been  previously  emphasized,  and  remains  to  be 
verified  experimentally. 

Thus  far  the  only  assumptions  about  the  disequilibrium 

have  been  that  it  is  odd  in  ' ^  (transverse)  and  small 

compared  to  the  eoui librium  populations.  Let  us  go 

further  and  assume  that  f.  may  be  described  by  f°(>,  -  r  ) , 

k  •  *  '  k  n 

where  r  is  different  from  >:  and  can  be  viewed  as  the 
n  s 

chemical  potential  of  the  quasiparticles.  In  effect,  we 


are  assuming  that  the  quasiparticles  come  to  equilibrium 
among  themselves  before  they  come  to  equilibrium  with 


the  nair  condensate.  Expanding  the  effective  charge  q^ 

in  a  Taylor  series  in  -r  (which  we  take  to  be  small 

s  n 

compared  to  A)  and  dropping  higher  order  terms,  we  have 


"k  *  qk  , *  <•„*•„>  — 


'kl>n  -  -j|, 

k  n 


(3.19) 


(3.20) 
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o 
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spends  to  local  e.’iitlibt  mm  with  the  condensat  e  .About  the 

shitted  ohemiu.Al  potential.  '.’his  is  .a  it't  l<vt  ion  of  the 

point  made  earlier,  that  ne.Ai  sm.A  l  l  changes  m  the 

o 

d  i  st  i  i  but  ion  function  m  oduoed  enh.Anoed  chanues  in  0*  .And 

henoe  In  wh.At  follows,  stnce  we  .no  retatninn  terms 

on  l  S'  of  le.Admn  order  in  k  „T,  we  will  take  :  0.  ’.'his 

H  n 

approx  im.Af  ion  should  remain  v.i  1  i  d  even  when  t  ^  ounnot 
striotly  be  expressed  in  terms  ot  .a  Fermi  tunot  ion. 

In  order  to  deal  with  sp.At  i.Ally  i nhomoueneous 
d  i  soon  ill  hr  i  .a  ,  we  define  eleot  roohemio.Al  pv'tenti.Als 
n  for  the  super  flu  id  condensate  .And  i.  for  the  quasi- 
particles,  .as  t he  sum  of  the  chemical  and  elect  t  teal 
not  ent  i  -a  l  s  : 


u 
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n 


t  +  Ov' 
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t  ►  e 
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We  then  use  the  standard  Ginzburg-Landau  expressions 
relating  ug  and  the  supercurrent  density  Jg  to  the 
phase  8  of  the  superconducting  order  parameter, 

Us  =  -%fie  (3.24) 

n  e 

=  n  evo  =  ~  (fi^8  -  2eA)  (3.25) 

s  s  s  2m 

to  obtain  the  following  relation  for  the  time  evolution 

of  3  : 
s 

T  i  =  -  -  oA  =  oE  +  D^Q*  (3.26) 

os  e  s 

Here  2  =  —  v <f>  -  A  is  the  electric  field, 

o  =  ne^i/m  =  2N(0)e^D  (3.27) 


is  the  normal-state  electronic  conductivity,  x  is  the 

2 

momentum-transport  relaxation  time,  and  D  =  Vpi/3  is  the 
normal  electronic  diffusion  constant.  The  supercurrent 
response  time  x  is  one  of  the  fundamental  characteristic 
times  in  our  treatment  of  the  problem,  and  is  defined  by 
_  4ttX^o  mo  _  n 

To  -  2  „  2  “  n  T  (3.28) 

c  nge  s 

This  can  be  evaluated  using  the  standard  normalizations 

for  the  superelectron  density  n  ,  or  equivalently  the 

s 

29 

formulas  for  the  magnetic  penetration  depth  \ ,  in  the 
clean  and  dirty  limits  near  Tc: 


Clean  limit 


(3.29) 


Dirty  limit 
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2 

We  have  assumed  here  that  n  «  A  is  constant  and  holds 

s 

its  equilibrium  value. 

For  the  normal  current  density  A  we  use  the 

n 

usual  relation 

J  =  of.  =  -  -?u  “  cA  (3.30) 

n  e  n 

but  with  the  understanding  that  it  is  a  good  approximation 
only  near  T  .  This  is  because  superconduct i ng  quasi¬ 
particles  are  not  simply  either  electrons  or  holes,  as 
they  are  in  a  normal  metal,  nor  even  particles  with  a 
well-defined  charge  q^,  but  rather  coherent  superpositions 
of  electrons  and  holes.  However,  auasiparticles  with 
energies  much  greater  than  A  (which  form  the  majority 

in  the  limit  T  +  T  )  can  be  considered,  at  least  operation 

24 

ally,  as  either  electrons  or  holes.  Pethick  and  Smith 

have  compared  the  normal  current  .1  =5-5  with  the 

ns 

current  carried  by  the  flow  of  quas ipart icles  (which 
they  call  5  ^)  ,  and  have  demonstrated  that  the  difference 
goes  to  zero  as  A/k  T,  as  do  the  corrections  to  Fq .  (3.30) 

D 

'’’he  physical  significance  of  t0  becomes  clear  if  we 
subtract  F.q.  (3.30)  from  Eq .  (3.26)  to  obtain 

’.A  -  A  * Dfo*  •  l3-n’ 

In  the  equilibrium  case  where  Q*=0 ,  a  sudden  change  in 
the  total  current  density  5=5  +  5  will  initially  be 

carried  as  normal  current.  This  will  decay  into  super¬ 
current  exponentially  with  time  tq,  so  that  to  is  the 
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equilibrium  response  time  of  the  supercurrent  to  a 
change  in  the  total  current. 

The  presence  of  means  that  charge  can  leave  the 
quasiparticle  system  in  a  given  region  of  space  not 
only  through  inelastic  scattering,  but  also  through  the 
net  inflow  or  outflow  of  quasiparticle  current  as  deter¬ 
mined  by  the  term  ?-J  .  Because  this  acts  to  produce  a 

n 

direct  change  in  the  distribution  function,  it  enters  as 
a  term  inside  the  brackets  in  Eq.  (3.16),  and  is  subject 
to  the  same  enhancement  factor  as  that  due  to  to  inelastic 
charge  imbalance  relaxation.  We  thus  have 


4k  T 


(3.32) 


Similarly  a  quasioarticle  injection  term  would  also  enter 
inside  the  brackets.  Wow,  if  we  take  the  divergence  of 
Eq.  (3.31),  add  tq  times  the  time  derivative  of  Eq.  (3.32), 
and  assume  *-3=0  from  electroneutrality,  we  obtain  the 
following  partial  differential  equation  involving  only 
the  charge  imbalance  O* : 

,2, 


Dtq* v  Q*  =  tqteQ*  +  (tq+te)Q*  +  Q* 


(3.33) 


This  is  the  charge  imbalance  wave  equation  we  have 
been  seeking.  In  the  low  freauency  limit  w<<x  x_  ^ 

O  hi 

it  describes  charge  imbalance  decay  with  a  characteristic 
decay  length 


V  =  ^ 


(3.34) 


This  length,  the  quasiparticle  diffusion  length  corresponding 
to  the  branch  imbalance  relaxation  time  Tnt,  has  been  used 
to  characterize  the  penetration  of  an  electric  field 
into  a  superconductor  in  the  presence  of  a  current  flow 
through  a  normal-superconduct ing  boundary . ^ ^ 7  ' 

From  here  on,  we  will  denote  this  length  simply  by  A, 
but  it  should  be  clearly  distinguished  from  the  related 
length  /Dr"  ,  for  which  the  same  notation  was  used  by  SBT. 

In  the  high  frequency  limit  (.<<1^  ^ ,  tq  ^) ,  Eq . 

(3.33)  describes  waves  of  charge  imbalance,  propagating 
with  velocity 


v 


a//Ve 


(3.35) 


Physically,  this  is  the  same  phenomenon  as  the  propagating 

collective  mode  observed  experimentally  by  Carlson  and 
2  3 

Goldman4-  in  the  tunneling  fluctuation  spectra  of 
superconducting  films  below  T  .  We  can  evaluate  Fq . 

(3.35)  by  substituting  in  for  x  from  Eo .  (3.29): 

O' 


8kBT  T  3, 

Vpl-yfj  (1-  £  )]* 

c 


0.53  v„  (1-  %  )*,  Clean  limit 
F  ‘c 

(3.36) 


(  ( 2DA/T0  ^  -  0.61  v  (£/t  )  ^  (1-  ~  )*  ,  Dirty  limit 


where  v_  is  the  Fermi  velocity,  i=v _r  is  the  mean  free 

r  r 

path,  and  o=hv_/r (0)  is  the  BCS  coherence  length. 

The  above  expressions  for  the  velocity  are  identical  to 

those  derived  for  this  propagating  mode  within  other 

22 

theoretical  frameworks. 
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For  the  more  general  case  where  neither  the  high 
nor  the  low  freauency  limits  is  appropriate.  Eg.  (3.33) 
describes  damped,  dispersive  waves  of  charge  imbalance 


with  the  dispersion  relation 


2,  2 


57,58 


-A  k  =  (1  +  iut  )  (1  +  imt  ) 

O  r. 


(3.37) 

for  waves  which  vary  as  exp[i(^t  -  k-r)]. 

Thus,  the  two  phenomena  of  dc  electric  field 
penetration  into  a  superconductor  and  the  high  fre¬ 
quency  propagating  collective  mode  are  both  special 
cases  of  the  more  general  phenomenon  of  charge  imbalance 
dynamics  in  a  superconductor,  governed  by  Eqs .  (3.31) 
to  (3.33).  We  will  show  in  the  next  section  that  these 
equations  in  turn  can  be  represented  by  a  physically 
transparent  transmission-line  equivalent  circuit. 
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3.3  A  Transmission- Line  Equivalent  for  Charge  Imbalance 
Dynamics 

When  expressed  in  one  dimension,  the  charge  imbalance 
wave  equation  (3.33)  is  well  known  historically  as  the 
telegraph  equation,  in  the  form 


3 2  E 

3^ 


=  LC  E  +  (RC  +  GL)  E  +  RG  E 


(3.38) 


This  describes  the  propagation  of  electromagnetic  waves 

down  a  cable  with  a  series  inductance  and  resistance 

per  unit  length  L  and  R,  and  a  shunt  capacitance  and 

leakage  conductance  per  unit  length  C  and  G.  An  early 

59 

treatment  of  the  problem  was  by  Kelvin  in  1855  to 
explain  the  wave-carrying  properties  of  submarine 
telegraph  cables.  In  general,  the  waves  are  not  only 
attenuated,  but  are  dispersed  as  well;  the  phase  velocity 
depends  on  frequency.  Since  dispersion  is  an  undesirable 
feature  when  sending  a  signal  down  a  cable,  it  was  noted 
early  on  that  in  the  special  case  that  RC=GL,  the  wave, 
although  still  attenuated,  is  dispersionless  with 
velocity  1//LC  for  all  frequencies. 

The  correspondence  between  these  telegraph  waves  and 
charge  imbalance  waves  can  be  made  more  explicit. 

Defining  the  gauge-invariant  nonequilibrium  potential  $ 


*  =  *  +  2i  6  =  'cs/e 


we  can  rewrite  Eqs .  (3.31) 


0* 

2N(0)e2 
and  (3.32) 


(3.39) 

in  the  form  (for 


-  ---•  ~ 


^6 


a  one-dimensional  wire  with  cross-sect ional  area  A) : 
31 


3  x 


+  G4>  +  C4*  *  0 


■  +  RI  -  LI  =  0 
3x  n  s 


where  the  parameters  (each  per  unit  lenqth) 


(3.40) 

(3.41) 


are : 


R 


L 


1_ 
a  A 
r 

o 
a  A 


G  =  2H  (0)  e  A/t 


C  =  2  M  ( 0 )  e  ■"  A 


Q* 

it  A 

4knT 


a  A 
2 


A 

aA 

7*  F 


(3.42) 


'Hiis  system  of  equations,  strictly  valid  only  in  the  limit 

A<<k  T,  is  represented  by  the  transmission-line  equivalent 

in  Fig.  3.2.  The  resistive  line  represents  the  channel 

for  the  normal  current  I  and  its  voltaqe  u  /e,  and  the 

n  n 

inductive  line  represents  the  condensate  channel  with  Ig 
and  l‘g/e  (we  can  iqnore  the  vector  potential  A  for  a 
quasi-one-dimens ional  superconductor).  The  voltaqe 
between  the  lines  /e  is  the  nonequilibrium 

potential  41.  Eqs.  (3.40)  and  (3.41)  then  follow  directly 
from  Kirchoff’s  laws. 

The  circuit  parameters  R,  L,  C,  and  G  have  a  simple 
and  appeal inq  physical  interpretation.  The  relaxation 
of  41  is  governed  by  the  time  constant  t  =  C/G,  while 

r. 

the  ratio  L/R  is  the  characteristic  response  time  i  for 

o 

the  supercurrent  when  <>=0.  The  resistance  per  unit  length 
R  is  simply  the  normal  resistance  per  unit  lenqth  of  the 
filament,  while  the  inductance  per  unit  length  L  is  the 
usual  London  theory  kinetic  inductance  of  the  super- 


r 


‘>7 


(*)UI  '  3/(X)Ut/ 
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current.  To  assist  in  the  interpretation  of  C  and  G, 


we  introduce  the  quantity'' 


^  =  4CT 

D 


(3.43) 


which  is  a  measure  of  the  departure  of  the  distribution 
function  from  global  equilibrium,  as  distinguished  from 
Q*,  which  measures  the  departure  from  local  equilibrium. 
Neglecting  conversion  processes,  we  obtain  from  Eq.  (3.32) 
the  continuity  equation  Q  =  -dJ  /dx.  Thus,  0  is  the 
quantity  which  is  analogous  to  the  "free  charge  density" 
in  the  transmission  line.  T’he  capacitance  per  unit  length 
C  is  then  equal  to  QnA/4> ,  which  differs  from  Q*A/4>  by 
the  factor  7rA/4kQT.  This  factor  can  be  interpreted  as 
an  effective  "dielectric  constant"  (although  it  is  less 
than  one  rather  than  greater  than  one  as  is  the  case  for 
the  usual  dielectric) ,  representing  the  fact  that  a 
change  of  the  "free  charge  density"  induces  a  much 
larger  change  of  Q* ,  due  to  the  adjustment  of  required 
to  maintain  overall  charge  neutrality.  The  conductance 
per  unit  length  G  gives  the  proportionality  between  the 
rate  at  which  QnA  changes  by  conversion  processes 
(=  ~Q*A/t  *  =  -Q  A/t  )  and  the  nonequilibrium  potential  4>. 


*  The  geometrical  inductance  remains  undefined  because 
the  return  path  for  the  current  is  not  included  in  this 
circuit.  For  our  usual  case  of  constant  total  current 


it  can  be  ignored. 


The  significance  of  the  equivalent  circuit  oarameters 
can  also  be  seen  from  the  energy  balance  equation  which 
follows  from  F.qs.  (3.40)  and  (3.41): 

in:R  +  $~g  +  ^-[kLI^2  +  %C^2 1 


x—  (I  u  /e  +  I  u  /e) 
Ox  n  n  s  s 


(3.44) 


The  first  tv/o  terms  represent  energy  dissipated  by  Joule 
heating  in  the  normal  resistance,  and  by  charge  imbalance 
decay.  Energy  is  stored  "inductively"  in  the  motion  of 
the  supercurrent  and  "capacitat ively "  in  the  field  of 
the  nonequilibrium  potential,  and  can  be  carried  into 
or  out  of  a  region  of  space  by  the  flow  of  normal  or 
supercurrent  at  their  corresponding  potentials. 

The  propagation  characteristics  of  the  transmission 
line  depend  on  the  relative  magnitudes  of  the  two  times 
t  and  r_.  When  t  =r„,  we  have  the  special  dispersion- 
less  case  mentioned  earlier  (L/R  =  C/G) ,  where  Fq .  (3.37) 


simplifies  to 


k  =  %  +  iki  -  Hr +  * 


(3.45) 


so  that  the  phase  velocity  and  decay  length  are 


w  A  ,  - 1  . 

v  '  k-  =  7  :  ki  *  A  • 

R  o 


(3.46) 


independent  of  frequency.  Otherwise  they  are  bounded  by 
the  low  and  high  frequency  limits  given  by 


A/ [ ( t  +t  ) /2J  =  v ( 0)  <  v  <  v(~)  =  A//t  T 
o  t.  or. 


(3.47) 


A  =  kT_1(0)  >  kT-1 


kI“1(»)  =  A/l0TF/(^r^r  (3'48) 
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For  illustration  wo  show  in  Fiq.  3.3  the  real  and  imaginary 

parts  of  k(m)  for  two  values  of  the  ratio  t  „/  r  Q . 

We  can  compare  these  characteristic  times  for  some 

representative  superconductors  by  evaluating  Eq.  (3.29) 

for  t  in  the  dirty  limit: 
o 


2k_T  h 
R  c 


5.2  "  10 


13 


sec 


(T  -T)  IK1 


(3.49) 


(r  here  has  the  same  temperature  dependence  as  = 

o  r  C.L 

nfi/  (  8k^  ( T  -T)  1  but  is  smaller  by  a  factor  of  n  / 1 4 (  3 ) 

=  5.79).  Thus  for  a  typical  experimental  temperature 

about  50  mK  from  T  ,  r  is  of  order  10  ^  sec.  The 

c  o 

electron-phonon  time  t  is  a  normal-metal  property 
and  hence  is  nonsinqular  at  Tc ,  so  that  we  can  take  it 
as  essentially  temperature- indendependent  for  any  qiven 


material  near  T  .  It  varies  from  2.3  '  10 
c 


■11 


Pb  ^  through  1.5  "  10  sec  in  Sn  to  1.3  v  10 


sec  in 
-8 


4 1 

sec  in  A1 ,  so  that  t 


r  is  the  more  usual  relation- 
o 


ship  except  very  close  to  T^.  The  dispersion  relations 
shown  in  Fiq.  3a  and  3b  correspond  to  En  and  Al  respective¬ 
ly  at  reduced  temperature  T/Tc  =  0.99. 

The  characteristic  impedance  7.  (m)  of  the  transmission 
line  can  be  found  by  standard  methods.  I f  we  constrain 

the  total  current  1=1  +1  to  be  constant  in  time, 

s  n 

then  the  charge  imbalance  wave  equation  (3.33)  is  obeyed 
not  only  by  4> ,  but  also  by  I  .  T  f  we  now  consider  a  wave 


■ 


Fig.  3.3  The  dispersion  relation  k(w)  =  k  +  k_  for 

H  1 

the  charge  imbalance  waves  [Eq.  (3.37)]  for  the  cases 

a)  t E  =  11  tq  (Sn  at  T  =  0.99  T  )  ,  and 

b)  tr  =  2R0  tq  (A1  at  T  =  0.99  T  ). 

For  a  given  ,  the  phase  velocity  is  w/k  and  the 

-1  ^ 

decay  length  is  kT 


102 


propagatinq  to  the  riqht  as  expliwt  -  ik(io)x)  ,  we  can 
Fourier  transform  Fqs.  (3.40)  and  (3.41)  in  space  and 
time  to  obtain 

ZU)  =  *(k,w) /In(k,w)  *  /( R+iwL)  /(G+iwCT 

-  /(l  +  iuao)/(l  +  iw,  F)  (3.50) 

For  a  wave  traveling  to  the  left,  on  the  other  hand, 

4>  =  -ZI  ,  so  that  in  general  one  is  justified  in  using 
n 

Fq .  (3.50)  only  for  a  semi-infinite  line  with  no  re¬ 

flections.  This  is  discussed  in  more  detail  in  Appendix 
III,  as  is  the  extension  to  short  lengths  of  transmission 
line.  For  the  special  case  r  =t„,  Z(w)  is  always  equal 

O  r* 

to  A/oA,  the  resistance  of  a  lenqth  A  of  the  line.  If 

t  ^t_,  it  is  equal  to  A/aA  for  low  frequencies  and 
o  t. 

(A/i'A)  >TyTr  for  high  frequencies,  but  will  have  a 
O  E 

reactive  component  for  intermediate  frequeicies.  The 
use  of  this  quantity  Z(u>)  will  greatly  facilitate  our 
discussion  of  the  phase-slip  center. 
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3 . 4  Discuss  ion 

The  transmiss ion- 1 ine  equivalent  which  we  have 
developed  to  describe  superconduct i no  charge  imbalance 
dynamics  Drovides  a  coherent  picture  in  which  to  view 
the  whole  range  of  phenomena  involving  a  nonequilibrium 
potential  in  a  superconductor.  In  order  to  highlight 
both  its  usefulness  and  its  limitations  we  will  briefly 
consider  its  applicabi 1 i ty  to  a  number  of  experiments 
that  have  already  been  described  in  the  literature. 

As  we  discussed  in  Chapter  I ,  the  concept  of 

9 

branch  imbalance  was  first  introduced  to  explain  the 

measurements  by  Clarke'1  of  differences  between  ;i  and 

n 

u  induced  in  a  superconducting  film  by  spatially  uniform 
quasiparticle  tunnel  injection.  Within  the  regime  near 
where  our  transmission-line  picture  is  most  valid, 
we  can  represent  the  physical  processes  involved  by  the 
equivalent  circuit  of  Fig.  3.4a  (cf.  Fig.  1.3a).  A 
constant  current  I  tunnels  from  a  normal  film  through  an 
oxide  barrier  into  a  superconducting  film  below.  This 
current  creates  "free  charge"  in  the  normal  channel  at 
a  rate  (Q^K  n ^  =  1/ ( ( A)  ,  charging  up  the  capacitors. 

They  discharge  .at  the  same  rate  through  the  leakage 
conductance  G,  so  that  the  current  I  leaves  the  nonequi¬ 
librium  region  as  pure  supercurrent  in  the  condensate 
channel.  The  induced  nonequilibrium  potential  is  therefore 


♦  = 


•  lT0./I2N(0)e2tAl 


k 


MHHWHMba 


1/  (GO 


(3.51) 


Use  of  the  transmission-line  eouivalent 
to  model  a)  uniform  cjuas ipart icle  injection 
and  b)  current  flow  through  a  normal-suoer- 
conductinct  boundary. 
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where  i  is  the  length  of  the  nonequi 1 ibr i um  region, 
much  larger  than  A  (to  minimize  edge  effects).  The  form 
on  the  right  of  Fq.  (3.51)  is  essentially  the  same  as 
that  originally  derived  by  Tinkham  and  Clarke  to 
explain  the  nonequilibrium  potential,  measured  using 
superconducting  and  normal  probes  in  the  nonequilibrium 
region.  We  reemphasize  here  that  if  the  current  is 
suddenly  turned  off,  the  potential  4>  will  decay  exponen¬ 
tially  with  a  time  constant  C/G  =  i„,  not 

I  vj 

One  can  also  see  from  this  picture  why  the  S  and  N 
probes  measure  their  respective  electrochemical  potentials, 
within  the  region  near  where  the  picture  is  valid. 

The  normal  detector  junction  is  similar  to  the  normal 
injector  junction,  except  that  there  is  negligible 
current  oassing  through  (assuming  a  high-imoedance 
voltmeter) .  therefore  there  is  no  voltage  drop  across 
the  oxide  barrier,  and  the  potential  of  the  normal 
probe  is  equal  to  u  .  A  superconducting  probe  may  be 
either  a  continuation  of  the  superconducting  film  out 
of  the  nonequilibrium  region  (as  it  was  for  the  Clarke 
experiment)  or  a  superconductor-superconductor  Joseohson 

1 o 

tunnel  junction  (as  in  the  Dolan-Jackel  experiment  ) . 

In  the  former  case,  ug  will  be  constant  in  space  in  a 
steady-state  situation,  provided  that  the  critical 
current  of  the  superconducting  film  is  nowhere  exceeded. 


lOh 


In  the  latter  case  we  have  essentially  the  same  condition. 

There  mav  be  a  difference  in  the  values  of  u  on  the 
'  n 

two  sides  of  the  S-S  junction,  so  that  there  is  quasi¬ 
particle  current  passing  through  the  hairier,  but  so 
long  as  a  counter  flow  of  supercurrent  is  possible,  u 
will  be  the  same  on  both  sides. 

Another  important  phenomenon  in  which  charge 
imbalance  plays  a  crucial  role  is  the  excess  resistance 
on  the  superconducting  side  of  a  current-carry incj 
superconducting-normal  interface.  We  can  represent  this 

(aqain  for  T  near  T  )  bv  the  equivalent  circuit  in  Fig. 

c 

3.4b,  in  which  the  normal  current  connects  up  to  the 
"normal"  channel  on  the  suoerconduct ing  side.  Over  a 
length  equal  to  A,  the  current  leaks  over  to  the  con¬ 
densate  line  by  charge  imbalance  decay,  and  the  excess 
resistance  is  the  normal  resistance  of  the  length  A. 
Furthermore,  the  equivalent  circuit  gives  rather  directly 
the  exponentially  decaying  nonequilibrium  potential 
near  the  boundary,  which  as  discussed  in  Chapter  I  was 
observed  by  Yu  and  Mercereau  ^  using  small  normal  and 
superconducting  probes  (see  Fig.  1.3b). 

As  in  our  model  for  the  PSC,  the  pictures  in  Fig. 

3.4  are  strictly  valid  only  very  near  T„.  For  lower 
temperatures  the  fact  that  quasiparticles  near  the  gap 


have  both  elect ron- 1  ike  and  hole-1  ike  characteristics 
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becomes  imnortant.  A  fraction  of  the  injected  current 
of  order  A/k  T  enters  the  condensate  line  directly  in 

u  C 

Fig.  3.4a/  reducing  the  contribution  to  the  charge 
g  13 

imbalance.  '  Similarly,  electrons  with  energy  <  A 

approaching  the  N-S  boundary  from  the  normal  side  in 

Fig.  3.4b  are  reflected  as  holes  via  Andreev  reflection, 

with  superconducting  pairs  to  carry  the  current  in  the 

superconductor.  This  corresponds  to  modifying  Fig.  3.4b 

so  that  some  fraction  vA/kpTc  of  the  current  from  the 

normal  metal  is  injected  directly  into  the  condensate 
37  18 

line.  '  A  microscopic  analysis  which  calculates 
this  fraction  of  current  Andreev-reflected  at  the 
boundary  has  been  used  successfully  to  explain  the 

observed  boundary  resistance  over  a  wide  range  of 

.  18 
temperatures . 

Within  the  restricted  range  in  which  the  pictures 
are  valid,  we  can  vise  them  to  consider  ac  effects  as 
well.  For  instance,  for  an  ac  current,  the  decay  length 
of  the  nonequilibrium  potential  becomes  shorter,  as  was 
noted  in  Eq.  (3.48).  Similarly,  an  ac  injection  current 
localized  at  a  point  alone  a  superconducting  strip  will 
qenerate  charge  imbalance  waves  which  propagate  down  the 
strip. 

We  wish  to  emphasize  that  this  physically  appealing 
transmission-line  picture  is  exactly  equivalent  in 
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physical  content  to  some  of  the  recent  generalized 

24  57 

two-fluid  theories  of  nonequilibrium  superconductors.  ' 
These  theories  also  include  only  the  leading  order  term 
in  A/kQT  and  are  therefore  valid  only  near  T  .  They 
also  ignore  all  pair-breaking  effects  of  the  superflow, 
and  so  are  restricted  to  currents  well  below  the  critical 
current.  Finally,  the  discussion  is  limited  to  a  de¬ 
scription  of  the  transverse  mode,  thus  ignoring  all 
other  effects  which  would  lead  to  changes  in  the  gap 
parameter  A.  Nevertheless,  within  their  range  of  validity, 
the  theories,  and  hence  the  transmission-line  equivalent, 
are  expected  to  contain  the  essential  features  of  the 
spatial  and  temporal  dependence  of  the  potentials  and 
currents  in  nonequilibrium  superconductors.  Even  outside 
the  range  of  quantitative  validity,  the  same  qualitative 
pictures  might  well  continue  to  provide  a  useful  guide 
to  the  understanding  of  superconducting  dynamics. 


CHAPTER  FOl'R :  A  NEW  MODEL  FOR  THF.  PHASE-SLIP  CENTER 


4 . 1  Elements  of  the  Mode  1 

The  current-induced  breakdown  of  superconductivity 
in  a  quas i-one-dimens ional  filament  has  been  shown  to 
occur  via  the  nucleation  of  spatially-localized  ohase- 
slip  centers . ^ We  will  consider  a  single  isolated 

PSC  (in  zero  magnetic  field),  centered  about  x=0 ,  subject 

•  1 
to  a  dc  current  bias  (1=0).  We  will  follow  SBT  in 

dividing  the  problem  into  two  oarts  on  two  different 
length  scales.  The  oscillations  in  the  magnitude  of  the 
gap,  which  are  necessary  for  the  chase-slip  process,  are 
assumed  to  be  restricted  to  a  core  region  with  size  £  f, 
which  is  small  compared  to  the  length  scale  /V  character¬ 
izing  variations  in  the  charge  imbalance.  We  further 
assume  that  it  is  possible  to  model  this  core  region  by  a 
discrete  Josephson  oscillator  of  zero  length  and  zero 
resistance,  with  a  current-phase  relation  of  the  form 

I  ( x=0 , t)  =  If  ( A  0 ( t ) )  (4.1) 

o  C 

where  AO  is  the  phase  difference  across  the  core,  f  is 
some  function  with  period  2 tt  and  amplitude  1,  and  I  is 
the  critical  current  of  the  core.  ’’’he  form  of  this 
function  is  undetermined  in  our  model,  although  we 
would  not  expect  the  purely  sinusoidal  dependence  since 
this  is  appropriate  only  for  weak  links  small  compared 
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l  10 

I  is  to  bo  taken  .is  oqu.il  to  ot  sliqht  lv 
c 

smaller  than  t  ho  oi  it  ical  current  ot  t  ho  adiaoent  tc'oions 
of  t  ho  t  i  lament  .  Our  approx  im.it  ions  consistent  ly  nouloot 
t  ho  pa  ir-break inq  ot  toots  associated  with  t  ho  presence  ot 
noar-critical  supot current s  in  those  ton  ions.  however, 
wo  boliovo  that  the  essential  features  ot  the  model  at o 
qua  lit  at i vo ly  cot  root . 

Since  dosephson  oscillat  ions  are  dt  i  von  by  differences 
in  the  condensate  e  loot  rochenti  ca  1  potential,  we  can  write 
the  voltauo  V(t)  dtivinu  these'  osc  t  1  1  at  ions  as 

V(t)  -  v  At'  ^  -In  ( x  -  0 f ,  t )  -  v.  (x  O', til 

“  (0 f ,  t)  -  4'  (0~ ,  t)  ,  (4..'1 

where  the  condition  that  i,  is  continuous  throuoh  the 

n 

ot  w in  follows  from  our  assumption  of  nealiuible  timet  ion 
resistance.  by  symmetry  the  magnitude  of  4  must  be  equal 
on  the  two  sides  ot  the  oriuin,  and  since'  V 1 1  '  is  not 
identically  zero  the  values  of  4>  on  the  two  sides  must 
have  opposite  siuns.  Therefore,  we  have 

V(t)  =  2$(0+,t)  2I  d,.'  o'"  *  4(0*,, id 

-  2J  du'  o'"'t  Z  (...)  Int0  ,„0  ,  (4.3) 

where  the  expression  in  terms  ot  .*  p.O  as  in  V‘o .  (1.S01 
can  be  iustitied  by  reuardinq  the  dosophson  oscillatoi 
as  a  source  of  waves  ptopae.it  inu  away  from  t  ho  ot  i.nn. 

We  are  as  sum  i  nq  here  that  there  is  not  anothei  l'Sr  01  a 
reflective  discontinuity  m  the  transmission  line  ne.11  bv 


in 


on  the  scale  of  A,  since  then  waves  would  also  be 
propagating  towards  the  oriqin  and  Eq .  (3.50)  would 

not  be  directly  applicable. 

The  above  physical  picture  is  equivalent  within  our 

transmission-line  representation  to  placing  an  ideal 

Josephson  clement  characterized  by  Fq .  (4.1)  in  the 

inductive  (condensate)  line  at  x=0,  as  shown  in  Fig. 

4.1a.  As  in  standard  transmission-line  theory,  the  effect 

of  the  semi - i n f i n i te  line  on  either  side  of  the  Josenhson 

oscillator  can  also  be  obtained  by  replacing  the  line 

with  its  character i stic  impedance  Z(..>)  as  in  Fig.  4.1b. 

The  impedance  shunting  the  Josephson  element  is  then  2Z(u>), 

and  Eq.  (4.3)  follows  directly.  For  the  current-biased 

case  being  considered  here,  all  the  current  far  from 

the  Josephson  element  flows  as  supercurrent .  In  Fig. 

4.1b,  this  is  modeled  by  connecting  the  external  current 

source  to  the  condensate  side  of  the  impedance  Z(m). 

(In  the  more  general  case  of  an  ac  total  current,  the 

current  far  from  the  Josephson  element  would  be  split 

between  I  and  I  in  accordance  with  the  ratio  between 
s  n 

R  and  iu'T.,  so  that  the  modified  Fig.  4.1b  would  include 
connections  to  both  sides  of  Z(m).  A  simple  general¬ 
ization  a l lows  us  also  to  consider  cases  where  the 
resistance  of  the  junction  region  itself  RT  is  not 
small  compared  to  Z (w) ,  e.g.  for  small  in-line  tunnel 
junctions,  as  in  the  simplified  schematic  in  Fig.  4c.] 


Fig.  4.1 
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Equivalent  circuits  for  a  phase-slip  center. 

a)  The  ideal  Josephson  element  is  located  on  the 
inductive  (condensate)  line  of  the  transmission  line 
of  Fig .  3.2. 

b)  The  equivalent  of  (a)  as  far  as  the  currents  and 
potentials  at  the  center  are  concerned,  where  Z(u>) 

is  the  complex  f reouencv -dependent  characteristic  im¬ 
pedance  of  the  transmission  line  [Fq.  (3.50)1. 

c)  A  generalization  of  (b)  to  allow  for  cases  where 
tire  intrinsic  resistance  RT  of  the  Josephson  element 
is  not  small  compared  to  7.  (w). 


1 


114 


The  high-vol tage  dc  behavior  of  the  PSC  tares  a 

particularly  simple  form.  At  the  center  of  the  PSC 

most  of  the  current  I'"I  must  be  normal  current  I  , 

c  n 

which  in  terms  of  Fig.  4.1b  flows  through  the  shunt 

impedance.  The  oscillations  in  I  of  order  I  are  then 

n  c 

relatively  small,  so  that  the  Josenhson  element  is  subject 
to  an  almost  constant  voltage  bias.  Therefore  the  phase 
difference  AO  across  the  element  increases  nearly  linearly 
with  time,  and  we  can  express  the  current-voltage  relation 
in  the  form 

V  =  2Z(0)  [I  -  ls(0,t)  ]  =  ^(1  -  Ic<f")  (4.4) 


where 


<f>  = 


2^  J  f ( AO) d ( AO) 


(4.5) 


is  the  average  value  of  the  periodic  function  f  and  the 
horizontal  bars  represent  averages  over  time.  In  the 
high-voltage  regime,  then,  the  model  predicts  a  constant 
differential  resistance 


^  =  2Z(0)  =  ^4  =  R 
dl  a  A  eff 


(4.0) 


and  an  extrapolated  excess  current  1  =  <-f'T  ,  in 

ex  c 

agreement  with  SBT .  The  high-voltage  I-V  relation  of 
Fq.  (4.4)  can  also  be  derived  simply  by  replacing  the 
impedances  Z(m)  in  Fig.  4.1b  by  their  dc  limits.  This 
is  equivalent  to  a  simple  resistively  shunted  junction 
(RSJ)  model  with  shunt  resistance  R  although  here  we 


2.ir . :■ 
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are  allowing  a  more  general  supercurrent-phase  relation 
that  the  sinusoidal  one  usually  associated  with  the  RSJ 
model.  Note,  however,  that  R  cc  of  the  PSC  is  not  the 
intrinsic  junction  resistance,  which  we  have  assumed  to 
be  negligible,  but  the  resistance  of  the  nonequilibrium 
region  which  extends  a  distance  A  on  either  side  of  the 
core . 


This  effective  resistance  R  ^  suggests  a  character¬ 
istic  voltage  vc=IcReff'  at  which  the  fundamental 
Josephson  frequency  is 


=  2e 

*cReff 


4el  A 
c 


(4.7) 


For  much  lower  voltages,  the  Josephson  oscillations 
become  much  more  pulse-like  and  are  dominated  by  higher 
harmonics  of  2eV/fi  uo  to  v  uic.  For  larger  voltages, 
on  the  other  hand,  the  Josephson  frequency  is  larger 
than  ui  ,  but  the  oscillations  in  this  regime  become 
negligibly  small  oerturbations  of  the  dc  average,  and 
do  not  significantly  affect  the  FSC  behavior. 

It  is  of  interest  to  compare  wc  ^  to  the  character¬ 
istic  times  tq  and  t  .  This  is  facilitated  by  expressing 

Reff  *n  terms  t^ie  Ginzburg-Landau  critical  current 
52 

of  the  strip 


IC  (strip)  =  Angevc  =  J  77F  .  (4.8) 

o 

in  the  form 
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so  that 


(D  T  = 
C  O 


3/3  C  I  (strip) 


(4.10) 


If  A>>£,  as  our  model  assumes,  and  Ic  is  only  slightly 


less  than  I  (strip),  it  follows  that  w  t  >>1 
c  co 


In  the 


usual  regime  t  >  t  ,  w  r  >>1  as  well.  In  subsequent 

r.  O  C  L 

sections  we  will  make  considerable  use  of  approximations 

involving  this  "high-frecmency  limit"  where  w  r  ,  co  r_>>l. 

CO  C  fcj 

In  Fig.  4.2  we  show  how  *  comoares  with  t  and  Tr 

c  o  t. 

as  a  function  of  temperature  near  T  ,  for  parameters 


appropriate  for  aluminum,  tin,  and  lead  in  the  dirty 
limit.  The  temperature  and  material  dependence  of  w 


-1 


takes  the  form 


-1 


1.8  x  10 


•18 


Tglsec]  Tc4[K] 


=££  (T  -  T)-5/4[K] 

4  r  r»  i  C 


(4.11) 


For  this  expression  we  have  taken  I  =1  ( strip) ,  with 
the  understanding  that  a  reduced  value  of  Ic  will 
increase  the  value  of  wc  ^  proportionately.  It  is 
evident  from  Fig.  4.2  that  of  the  three  materials,  Al 
satisfies  the  approximations  in  the  high-frequency  limit 


*  For  the  case  of  an  ideal  Josephson  weak  link  (tunnel 

junction  or  small  metallic  constriction)  with  critical 

current  I  and  intrinsic  junction  resistance  RT,  it  has 

been  shown"^  that  I  R  =  h/(2er  ).  Therefore,  the 

C  J  o 

neqlect  of  R_  in  R  c,  may  be  valid  only  when  w  t  >>1. 

J  0LL  CO 
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acteristic  times  Tp,  t  ,  and  ^  1  as  a 
tion  of  temperature  near  Tc  for  Al ,  Sn, 

Pb  in  the  dirty  limit,  according  to  Eqs. 
9)  and  (4.11).  The  values  of  Tp  are  from 
.  41  (Al) ,  36  (Sn) ,  and  40  (Pb) . 


I 


-i  I  '*V 
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over  the  widest  range  of  temperatures.  They  are  approx¬ 
imately  true  for  Sn  beyond  about  10  mK  from  T^,  and 
generally  untrue  for  Pb.  In  the  clean  limit,  tq  is 

increased  by  a  factor  ~l/E,  and  to  *  by  '(C/g  )  ^  while 

o  c  o 

te  remains  the  sane,  which  have  the  effect  of  increasing 
the  temperature  range  over  which  the  approximations  in 
the  high-frequency  limit  remain  valid. 

Returning  to  the  behavior  of  the  phase-slip  center, 
we  will  concentrate  our  attention  on  the  right-hand  side 
xs0,  since  the  solution  for  the  other  side  follows  from 
symmetry  considerations.  Since  the  disequilibrium  is 
localized  near  the  origin  and  the  total  current  is 
constant  in  time,  the  current  becomes  pure  supercurrent 
far  from  the  PSC.  The  potentials  out  in  this  region  are 
then  constant  in  space,  and  we  are  free  to  set 


U n  ( ‘^ )  =  Ug(w)  =  0 

We  then  have 

Ug(x,t)  =  eL  C  Igdx  =  -eL  f  Indx 
J  x  J  x 


(4.12) 


L  •  ,  .. 

-  -  5  wn(x,t)  = 

This,  together  with  the  fact 
all  times,  implies  that 

ug(x)  =  0  ,  x  >  0 

Since  the  dc  decay  length  is 


~Tov'n 

that  vi n  is  bounded  over 


(4.14) 

A,  the  other  time  averages 


are : 
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I  (x)  =  I  -  I  (x)  =  I  (x=0)e~x/A  (4.15) 

•  •  3  il 

Un(x)/e  =  <F(x)  =  ln(x=0)^-  e~x/A  (4.16) 

This  behavior  is  plotted  in  Fig.  4.3  together  with  the 
corresponding  behavior  on  the  left  side  of  the  PSC. 

Note  that  the  currents  are  symmetric  with  respect  to 
the  origin,  while  the  nonequilibrium  potential  is  anti¬ 
symmetric.  This  time-averaged  behavior  of  the  currents 
and  potentials  is  the  same  as  that  in  the  SBT  model, 
which  was  confirmed  experimentally  by  Dolan  and  Jackel. 
In  the  next  section  we  examine  the  associated  time- 
dependent  behavior  within  our  transmission-line  model 
for  the  phase-slip  center. 


19 
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4 . 2  Waves  in  the  Vicinity  of  the  PSC 

We  first  direct  our  attention  to  the  oscillations 

at  the  origin,  and  look  for  a  solution  with  I  oer iodic 

n 

in  time  with  a  fundamental  frequency  1  =  2eV/f .  Expanding 

in  a  Fourier  series  with  coefficients  c  ,  w'e  have 

m 

09 

I  (x-0,t)  =Vc  (  ')  eln-fc  (4.17) 

n  /  j  m 

where 

co  =  in(0)  -  V/Reff  (4.19) 

Applying  the  results  of  Appendix  II  [Fqs.  (A2.18)  and 
(A2.20)]  for  a  Josephson  element  shunted  by  an  impedance 
2z(,o)  (see  Fig.  3.1b)  ,  we  obtain 

q I  Jn/n 

cm(')  - - 2^  /  dt  e~im  tf<At>(t)),  m/0 ,  (4.19) 

"0 

where 

A?(t)  =  It  +  s  (m ')  -r-4  (einv:t-l)  .  (4.20) 

h  ‘  ■>  imu 

m^O 

Given  any  relation  f(.V)  ,  these  equations  provide  a 

formal  solution  to  the  problem,  which  can  be  evaluated 

on  a  computer  by  iteration,  although  the  existence  of 

a  solution  is  not  necessarily  guaranteed.  We  will  deal 

later  with  a  particular  f(A0)  for  which  these  equations 

can  he  solved  analytically. 

Given  the  set  of  these  Fourier  coefficients  c  , 

m 

we  can  express  the  snace  and  time  dependence  of  all  the 
physical  quantities  of  interest  in  the  problem.  Con¬ 
sidering  the  right-hand  side  x^O,  we  have 
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I  (x,t)  =  I  -  I  (x , t )  Tc  el(mnt  "  kmx)  ,  (4.21) 
n  s  i  -d  m 


where 


m 


k  =  k  +  ik 
m  m_  m_ 

R  i 


k  (mil) 

is  the  wave  vector  from  Eq .  (3.37)  ,  and 


k  <0  ;  k  -  0  for  m  0  , 

ml  V 


(4.22) 


(4.23) 


which  corresponds  to  decaying  waves  propagating  to  the 

right.  The  notentials  4> ,  y  ,  and  y  can  likewise  be 

n  s 

expressed  in  terms  of  the  c  's: 

m 


r,  .  M  „  l  ( mil  t  -  k  x 

Z  (m.. )  c  e  m 

m 


4>  (x ,  t)  =  y  ^ 

m 

yr»  OO 

I  (X'  ,t)dx 
X 

eV~^cm  i(mHt 
cA/  ,,ik  e 

m 


(4.24! 


m 


(4.25) 


m 


us(x't)  =  *Toyn(x't) 
t  e 


o  V'  imft  i(mQt  -  k  x)  . 

Fa  IF“  cm  e  m  (4-26) 

m 

m 


The  dependence  for  negative  x  can  be  obtained  by  symmetry 
the  currents  1^  and  Ig  are  symmetric  on  either  side  of 
the  origin,  the  potential  4>  is  antisymmetric,  and 


yn(_x)  '  yn(0)  =  “  tun(x)  "  un(0)] 


(4.27) 


’Js(-X)  “  un(0)  =  ~  [us(x) 


un(0) 


The  apparent  asymmetry  between  the  behavior  of  y  (and  yg) 
at  x  =  t«>  arises  from  our  arbitrary  boundary  condition 


that  o  (*»)  “  ns  («)  *  0.  Since  only  potential  differences 
are  important,  this  is  of  no  concern. 

The  condition  of  constant  current  implies  that  the 
voltage  V  across  the  PSC,  as  measured  by  probes  far 
away  from  the  nonequilibrium  reqion,  must  include  an 
ac  component : 

V.  *  hi  (-"■)  ~  un("’)l/e  =  2iin(0)/e 

-  V  t  -4  /  c  elmi2t//(l  +  iimli  )  ( l  +  inu2r_)  .  (-1.28) 
n  A  /  >  m  o  F 

m/0 

This  is  annloqous  to  the  ac  voltaqe  which  appears  in  the 
RSJ  model,  and  follows  from  inteqrating  the  (spatially 
symmetric)  oscillating  normal  current. 

The  linear  relation 

f  (AO)  =  A0/2n  ,  0  <  AO  <  2n  ,  (-1.29) 

which  is  one  branch  of  a  periodic  sawtooth  curve  passing 
between  zero  and  one,  has  the  unique  advantage  of  yield  inn 
an  analytic  solution  for  the  Fourier  coefficients. 

Alt  hough  its  instantaneous  phase-slip  at  the  end  of 
each  cycle  may  be  physically  unrealistic,  similar  dc 
supercurrent-phase  relations  have  been  measured  and 

50 

have  been  anp 1 ied  to  explain  properties  of  nncrobr idges . 

If  we  assume  here  that  the  system  remains  on  a  single 
branch  for  the  entire  period  (i.e.  there  is  a  single 
phase-slip  event  at  the  beginning  of  each  period),  then 
the  integral  in  Fq.  (*1.19)  can  be  evaluated  explicitly 
(see  Apoendix  T  T )  to  yield 


I  .'4 


Fig.  4.4  Tin e  deoendence  within  a  single  Josephson 

oeriod  tT  of  the  oscillations  in  the  center  of  (a)  I  (0,t)/l 
J  n 

and  ( b )  l(0,t)/I  *1  ,  ,  using  the  linear  f  (AO)  and  parameters 
c  ef f 

te/to  =  11,  u>ctq  =  2.2,  v  «  Ic?eff,  and  I  =  1.5  Ic,  appro¬ 
priate  for  Sn  at  t/t  =  0.99.  Both  I  and  41  are  dis- 

c  n 

continuous  at  the  beginning  and  end  of  each  cycle  for  the 
linear  f  ( AO) . 
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Fig.  4.5  Spatial  dependence  (for  x  >  0)  of  (a)  In/Ic 
and  (b)  at  t  =  0.25  for  the  same 

parameters  as  in  Fig.  4.4.  The  dashed  curves 
are  the  dc  averages. 


. 
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The  discontinuous  nature  of  these  jumps  is  again  an 
artifact  of  the  linear  current-phase  relation,  but  here 
they  make  the  progression  of  the  waves  easily  visible. 
Their  wavelength  L  is  given  in  the  high-frequency 


approximation  as 


L  =  vr  = 
•J 


/x  t  2eV 
o  E 


(4.31) 


which  works  out  to  =  0.85  A  for  the  case  of  Fig.  4.5. 
Increasing  the  voltage  would  decrease  the  wavelength 
proportionately.  Since  the  high-frequency  limit  is  only 
approximately  valid  here,  the  phase  velocity  depends 
somewhat  on  frequency,  producing  a  waveform  whose  shape 
evolves  slightly  as  it  propagates  down  the  line.  The 
corresponding  behavior  of  un(x,t)  and  us(x,t)  is  shown 
in  Fig.  4.6. 

In  Fig.  4.7  we  present  the  nonequi librium  potential 
0 (x>0 , t=Tj/4)  which  would  arise  from  the  classic  sinu¬ 
soidal  f(A0),  again  using  the  same  parameters  as  in 
Figs.  4.4  -  4.6.  This  solution  was  obtained  by  iteration 
on  a  computer.  Here  the  waves  have  no  discontinuous 
jumps,  since  f ( A 9 )  is  continuous.  Nevertheless,  the 
essential  picture  of  propagating  and  attenuating  waves 
remains  the  same. 


As  we  noted  in  Fa.  (3.43) ,  when  t  the  ac  decay 

O  E 

length  A  becomes  less  than  the  dc  decay  length  A. 

3C 

In  the  high-frequency  limit,  for  the  case  where  TE>>T0' 


o 

Fig.  4.6  Spatial  dependence  of  the  potentials  un/eV 
and  ug/eV  for  the  same  set  of  parameters  as  in  Fig.  4.5. 
The  behavior  for  x  <  0  is  included  to  make  clear  the 
symmetry  of  the  solutions.  The  range  of  variation  of 
u(-°°)/ev  is  indicated  by  arrows  on  the  left. 


The  spatial  dependence  of  $/IcR  ^  for  the 
same  parameters  as  in  Fig.  4.5  (except  for 
I  =  1.0  Ic) ,  but  using  a  sinusoidal  f ( A 8 ) . 
The  dashed  curve  is  the  dc  average. 
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we  have 

A 


ac  *  2A  /V’TE  ' 


(4.32) 

which  is  about  0.6  A  tor  Pins.  4.5  -  4.7.  The  difference 


between  two  two  lengths  is  more  dramatic  for  aluminum, 

where  is  about  80  times  longer  than  it  is  for  tin. 

This  is  demonstrated  in  Fig.  4.8,  where  I  (xs 0 , t« t  /2) 

n  J 

is  plotted  for  Al  with  T/T  =  0.90,  V  =  I  R  ,./2,  and 

c  c  eff 

I  =  I  (strip) .  In  this  case  we  have  A  -  A/8,  although 
c  c  ac 

we  emphasize  that  in  terms  of  absolute  length  A  is 

J  ac 

independent  of  r£  (in  the  high-frequency  limit)  and  has 
not  decreased;  rather  it  is  A  «  /FT  which  has  increased. 


Fig.  4.8 


SDatial  deoendence  of  I  ( x ,  t  ~  0 .5t  » )  / 1  for 

n  J  c 

the  linear  f(A6)  and  oarameters  t_  =  280  t  , 

r.  o 

cox  =  11,  V  *  0.5  I  R  and  I  =  1.0  I  , 
co  c  eir  c 

appropriate  to  Al  at  T/Tc  =  0.99.  Note  that 
the  horizontal  scale  is  expanded  by  a  factor 
of  two  compared  to  Fig.  4.5. 
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4 .  i  I-V  Ch  a  r  ac  t  or i s  t i cs  for  the  Phase-Si  ip  Center 

The  I-V  relation  for  a  general  f(.v3)  can  bo  expressed 
in  the  form 

I  -  I  (x=0 , t=0)  +  I  (0,0) 
n  s 

=  7  c  (f!)  +  I  [  f  ( A  0  ( t  =  0  + )  )  +  f  ( A  0  ( t  ~  0  ~ )  )  ]  2  .  (4.33) 

m*-<» 

Since  .0  =  2oV  Ti  this  allows  an  explicit  expression  for  I 

as  a  function  of  V.  (The  second  term  allows  for  a 

discontinuity  of  f (At4)  at  t  =  0.)  Curve  b  in  Fin.  4.o 

shows  I (V)  for  the  same  set  of  parameters  used  in  Fins. 

4.4  -  4.6,  namelv  the  linear  f  ( A  )  with  „•  2.2  and 

c  o 

* 

t..  i  =11,  cor  respond  i  no  to  tin  at  0 . 0  o  T  . 

!  O  c 

Fven  without  solving  explicit lv  for  the  c  's, 
however,  we  can  obtain  approximate  expressions  for 
I (V)  which  permit  a  qualitative  understanding  of  its 
behavior,  by  comparison  with  the  resistively  shunted 
junction  (FFJ)  model.  Since  the  terms  "PSJ"  lias  tra¬ 
ditionally  implied  a  sinusoidal  f(AP)  ,  we  will  refer  tv' 

*  Curve  b  ends  abruptly  at  a  point  in  the  negative 

differential  resistance  region.  For  lower  voltages,  the 

formal  solution  obtained  by  using  the  Fourier  coefficients 

of  F.q .  (4.30)  corresponds  to  an  oscillation  with  AO  N  2- 

(and  I  s  I  )  over  a  portion  of  the  cvcle.  It  we  attempt 
sc 

to  solve  Fqs .  (4.10)  and  (4.20)  by  computer  iteration 

using  the  actual  sawtooth  f (AO) ,  the  c  's  do  not  converge 
to  a  self-consistent  solution,  implying  that  no  such 
solution  exists  in  this  region. 


limit  [Eq.  (4.36)],  (b)  the  actual  solution  for  wcxo  =  2.2  [Eq.  (4.33)],  and  (c)  the 
low-frequency  or  GRSJ  solution  (Eq.  (A2.13)]. 
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a  "generalized  RSJ"  model  (GRSJ)  which  is  unrestricted 
in  this  regard.  (See  Appendix  II  for  a  detailed  expo¬ 
sition  of  the  GRP J  model.)  Curve  c  of  Fig.  4.9  shows 
the  GRPJ  curve  [Fq.  (A2.13)]  for  the  linear  f  (At?)  and 

a  shunt  resistance  R  ,,  =  2A/oA.  In  the  limit  of  hiah 

et  r 

voltages,  curves  b  and  c  agree  [compare  Fas.  (4.4)  and 
(A2.9)].  However,  in  the  range  shown  in  Fig.  4.9, 
curve  c  deviates  substantially  from  the  exact  result 
of  curve  b.  This  occurs  because  of  the  frequency 
dependence  of 


2  ff  =  2Z  (co)  =  R  ►TT+iu-T0)/(l+iwT  )  (4.34) 

which  is  the  shunt  impedance  equivalent  to  the  trans¬ 
mission  lines  on  either  side  of  the  core,  as  in  Fig. 

4.1b.  For  the  soecial  case  t  =  t„,  Z  cc  =  R  cc  for 
*  o  E  ef  f  ef  f 

all  frequencies  and  the  GRSJ  results  Eg.  (A?. 13)  is  the 

exact  solution  for  our  model. 

For  f  t  ,  a  simple  approximation  that  does  much 

better,  and  becomes  exact  in  the  high-frequency  limit 

to  t  ,  io  t  >>1,  is  to  consider  the  dc  resistance  to  be 
CO  C  F 

R  cc  but  the  ac  resistance  to  be 
ef  f 


R  =  R 
ac 


*ef  f'r\7~x  F  ’ 

The  solution  is  then  aiven  bv 


(4.35) 


I  =  I 


GRS 


SJ<*'  Rac>  +  R- 


-)  , 


(4.36) 


eff  ac 

that  is,  by  the  inverse  function  of  v  e  (I,R  )  with  a 

or.Du  «C 


Fig.  4.10  I-V  characteristics  in  the  high-frequency  limit, 
from  Eq.  (4.36),  for  the  linear  f(A6)  and  various 


values  of  t_/t  . 

F  o 


The  curves  are  scaled  so  that 


all  have  the  same  critical  current  and  high- 
voltage  differential  resistance.  The  values 
of  Tj?/T0  corresponding  to  the  curves  from  left 
to  right  are  1024,  256,  64,  16,  4,  1,  i,  and  1/16 
The  dashed  line  is  the  high-voltage  asymptote. 


(i.e.  it  has  a  larqer  I (v) ]  ,  while  in  the  more  usual 

case  i,,  >  i  ,  it  lies  to  the  left  .  This  is  true  for 
E  o 

any  current-phase  relation,  and  since  the  GtSJ  curve 
rises  vertically  at  I  I  ,  i  ,,  N  1  orovides  the  condition 
for  hysteresis. 

This  hysteresis  can  be  qualitatively  understood 
in  terms  of  the'  fact  that  the  effective  ac  resistance 
shuntinq  the  junction  is  smaller  than  t he  dc  resistance 
[ Fq .  (4.35)|  in  the  reqime  1  >  1  in  which  hysteresis 
is  present.  The  amplitude  of  voltage  oscillations 
across  the  junction  is  therefore  reduced  relative  to  the 
dc  avoraqo  value,  so  that  the  voltaqe  drivinq  the 
dosephson  oscillation  is  closer  to  beinq  a  constant 
voltaqe  bias.  The  resulting  I -V  curve  then  lies  closer 
to  the  vol taqe-biased  line  which  forms  the  hiqh-voltaqe 
asymptote,  and  as  is  evident  from  Fiq.  4.10  approaches 
this  line  at  lower  and  lower  voltages  as  t  /i^  is 
increased.  The  maximum  amount  of  hysteresis,  therefore, 
corresponds  to  the  curve  remaining  on  the  voltage-biased 
line  all  the  way  down  to  /.oro  voltaqe: 


m 


f"I  as  i _/ i 
c  F  o 


where  I  is  the  current 
m 


hysteret ical 1 y  back  into 


*  "  ,  (4 . 30) 

at  which  the  system  switches 
the  zero-voltage  state.  Even 


when  the  high-frequency  limit  does  not  strictly  hold. 


1  Mt 


.is  in  curve  b  ot  Fig,  4  .  *> ,  •  •*  i  t  (>m,i  i  n«  the  oondit  ion 

I  o 

tot  the  presence  ot  hysteiesis,  since  the  maun  \  t  ude  ot 
?■  .f(m  J  is  still  small  cm  than  R  ....  .  The'  amount  ot 
hysteresis  is  somewhat  t  educed ,  hut  this  is  not  a 
universal  quantity  even  in  » he  high-frequency  limit, 
depending  on  the  details  ot  the  curt  out -phase  relation. 

We  note  that  the  same  argument  can  also  explain  oualt 
tat ively  the  well-known  hysteresis  in  t he  oapaeitivol^ 
shunted  RSd  model ,  when'  the  et  toot  ivo  impedance 
1' ,  (1  *  itoR^i’.j)  also  decreases  in  magnitude'  with 
i nc re as  inn  t  r equency . 

In  an  actual  oxpet  iment  ,  one  tyoioally  measures 
voltage  as  current  is  swept  out  toi  a  s<'i  les  ot  tempeia 
tines  near  T,.  In  Fig.  4.11,  we  plot  theoretical  curves 
from  our  model  in  an  analogous  way,  showing  the  onset  and 
development  ot  hysteresis.  They  correspond  to  t  tie 
linear  current -phase  relation  in  the  h i uh- f t equency  limit, 
and  were  obtained  from  curves  such  as  those  in  Fig.  4.10 
by  using  the  temperature  scaling  iclat  ions 

i»/t  ,  «  ( 1 -T/T  ) 

F.  O  C 

t  ' 

1  "  ( 1 - T/T  1  l  4  .  PM 

c  c 

T  R  ,  r  "  (1  -T  T  T'4 

C  off  o 

Within  the'  17  i  von  npptoxinmt  ions  t  hose  nu  vos  ,11  r  unto 
pendent  of  UMforidl.  The'  onset  of  hystotosis  \  x%  \  \ 

r.  O 

would  correspond  to  T  -T  4  mK  foi  Sn,  nnd  0.0S  mK  tot  M 

c 
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Fig.  4.11  I-V  curves  for  the  linear  current-phase 
relation  in  the  high-frequency  limit  for  various  values 
of  te/t0»  with  the  temperature  dependence  reintroduced 
in  accordance  with  Eq .  (4.39).  The  arrows  pointing 

upward  represent  hysteretic  switching  at  the  critical 
current  I  into  the  dissipative  state,  and  the  downward¬ 
pointing  arrows  represent  switching  back  into  the  zero- 

voltage  state  at  the  current  I  .  The  values  of  t„/t 

m  E  o 

corresponding  to  the  curves  from  left  to  right  are 
k i  1 ,  2 ,  3 ,  4 ,  and  5 . 
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The  development  of  hysteresis  in  Fig.  4.11  is 

qualitatively  similar  to  that  seen  in  many  experiments. 

Quantitative  comparisons,  however,  are  complicated  by 

a  number  of  factors  which  we  have  not  considered  in  our 

model,  such  as  noise-induced  smearing  of  I-V  curves 

close  to  T  and  heatinq  effects  farther  from  T  .  In 
c  c 

particular,  heating  has  been  shown  to  account  reasonably 

well  for  hysteresis  in  some  cases.^'^'^'t'''  The 

intrinsic  hysteresis  which  our  model  predicts  should  be 

viewed  not  as  an  alternative  explanation  of  hysteresis, 

but  as  a  supplementary  mechanism  which  should  become 

important  near  T  ,  where  the  heating-induced  hysteresis 

disappears  more  rapidly. 

We  can,  however,  make  some  rough  preliminary 

comparisons.  In  our  model  for  intrinsic  hysteresis, 

the  normalized  return  current  I  /I  saturates  at  a 

m  c 

fixed  value  <fN  with  decreasinci  temperature,  whereas 

that  associated  with  heatinq  decreases  towards  zero  as 

the  PSC  eventually  becomes  a  sel  f-heat.i  ng  normal  hotspot. 

In  thin-film  microstrips  of  tin,  the  latter  behavior 

is  typically  seen.  In  clean  tin  whisker  crystals,  however, 

where  heating  effects  should  be  somewhat  less  severe, 

64 

Tidecks  and  Meyer  have  noted  a  saturation  effect  similar 
to  that  predicted  by  our  model.  The  onset  of  hysteresis 
near  Tc  is  another  basis  for  comparison.  This  onset  has 
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been  observed  in  the  range  from  1.5  mK  to  10  mK  from  Tc 
for  both  tin  whiskers^  and  tin  microstrips ,  roughly 
consistent  with  the  condition  t  =  x_  (about  4  mK  in 

O  r i 

Fig.  4.2),  although  rather  large  samole-to-sample 

variations  remain  to  be  explained.  Comparable  data  for  tin 

microstrips  in  an  electrically  unshielded  environment  do  not 

21  51 

show  hysteresis  until  significantly  farther  from  T  ,  ' 

presumably  because  of  rf  noise  rounding.  Because  of 
these  effects,  a  careful  set  of  measurements  on  hysteresis 
near  Tc  in  well-characterized  and  well-shielded  samples 
are  needed  to  provide  a  more  detailed  and  definitive 
comparison  with  our  theory. 

In  summary,  we  have  developed  a  new  model  for  the 
phase-slip  center,  which  is  based  on  coupling  Josephson 
oscillations  in  the  core  of  the  PSC  to  the  transmission¬ 
line  equivalent  developed  in  Chapter  III  to  describe 
charge  imbalance  dynamics.  When  averaged  over  time, 
our  model  reduces  to  the  SBT  picture  of  the  spatial 
dependence  of  the  potentials  inside  a  PSC,  which  has 
been  experimentally  confirmed  by  Dolan  and  Jackel. 

Another  important  consequence  of  our  model  is  the 
predicted  existence  of  intrinsic  hysteresis  in  the 
dc  I-V  curves,  and  preliminary  comparisons  with  experiment 
are  promising  but  complicated  by  heating  and  other 
factors.  We  have  also  made  preliminary  investigations 
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into  the  possible  application  of  our  model  to  the 
interaction  of  closely-spaced  PSC's,  with  some  encouraging 
results  (see  Section  5.1).  in  conclusion,  our  model 
provides  a  qualitative  understanding  of  ac  effects 
in  PSC's,  even  thouqh  its  quantitative  validity  remains 
uncertain . 


CHAPTER  FIVE:  OTHER  TOPICS 

The  previous  chapters  could  not  of  course  cover 
all  aspects  of  nonequilibrium  phenomena  in  superconducting 
phase-slip  centers.  In  the  present  chapter  we  round  off 
the  main  body  of  this  report  with  a  discussion  of  some 
additional  topics  which  were  dealt  with  in  the  course 
of  our  investigations.  Although  the  treatment  is  not 
as  complete  as  that  in  the  earlier  chapters,  it  should 
help  provide  guidance  for  future  research  in  the  area. 

5 . 1  Detection  of  Charge  Imbalance  Waves  and 

Interactions  Between  Phase-Slip  Centers 

How  does  one  go  about  detecting  high-frequency 

waves  in  a  superconductor ?  Two  ideas  come  to  mind. 

2  3 

One  is  the  Car lson-Goldman  experiment,  where  a  super¬ 
conducting  tunnel  junction  was  made  to  serve  as  a 
frequency  and  wave-vector  selective  detector  of  thermal ly- 
activated  collective  waves  by  application  of  a  dc  voltage 
and  an  external  magnetic  field.  The  other  is  the 
detection  of  microwaves  by  the  appearance  of  Josephson 
steps  in  a  Josephson  junction.  These  two  methods  are 
fundamentally  rather  similar,  and  in  the  present  case 
one  can  use  one  Josephson  element  (e.o.  a  small  tunnel 
junction  probe)  in  a  similar  way  as  a  selective  detector 
of  waves  generated  by  another  Josephson  element  (e.q.  a 
phase-slip  center  in  a  long  filament). 
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Consider  Fig.  5.1a,  where  a  superconducting  microstrip 
contains  a  localized  weak  spot,  and  nearby  on  the  scale 
of  the  ac  decay  lenqth  (of  order  the  auasiparticle 
diffusion  length  A)  is  located  a  small  Josephson  tunnel 
junction.  If  the  current  flowing  in  the  strip  exceeds 
the  critical  current  of  the  weak  spot,  then  a  PSC  will 
nucleate  there,  and  will  generate  charge  imbalance 
waves  which  extend  out  beyond  the  tunnel  junction  probe. 

In  particular,  the  condensate  electrochemical  potential 
Ug  will  oscillate  with  fundamental  frequency  ft  =  2eV/fi 
(Eq.  (4.26)],  where  V  is  the  total  dc  voltage  across  the 
PSC.  If  in  addition,  the  upper  electrode  of  the  tunnel 
junction  is  biased  at  a  dc  voltage  V'  relative  to  the 
lower  one  (via  the  superconducting  side-tab  in  Fig.  5.1a), 
then  the  potential  difference  Aug  across  the  junction 
will  consist  of  the  sum  of  a  constant  term  eV'  and  a 

•  i  Qf 

nonequilibrium  oscillatory  part  <*  e  "  .  Then,  since 

the  phase  difference  AP  across  the  junction  evolves  as 

fiAG  =  —  2 A vi  ,  AO  will  consist  of  a  linear  term  and  an 
s 

oscillatory  term.  Because  of  the  nonlinear  mixing 
properties  of  the  Josephson  current  through  the  junction, 
a  time-average  supercurrent  will  flow  given  the  usual 
conditions  for  an  ac  Josephson  step,  namely  V'  =  nfift/2e 
=  nV.  Taking  the  amplitude  of  the  nonequilibrium 
oscillations  from  Section  4.2,  we  find  that  the  resulting 
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Fig.  5.1  Detection  scheme  for  charge  imbalance  waves 

near  a  P SC.  a)  Physical  picture,  b)  Schematic 
c)  Equivalent  of  b)  ,  emphasizing  the  similarity 
with  the  configuration  for  investigating 
coupled  microbridqes . 
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step  size  is  at  best  several  percent  of  the  tunnel  junction 
critical  current  for  any  reasonable  set  of  parameters. 

This  is  small,  but  it  should  still  be  observable.  Further¬ 
more,  a  closely-spaced  array  of  such  tunnel  junction  probes 
would  indicate  directly  the  decay  length  of  the  charge 
imbalance  waves,  and  thus  distinguish  this  effect  from 
that  due  to  any  possible  microwave  coupling  between  the 
PSC  and  the  tunnel  junctions. 

A  number  of  other  issues  should  be  addressed  in 

connection  with  this  proposed  experiment.  First,  the 

current  in  the  tunnel  junction  will  include  not  only 

a  supercurrent  but  also  a  quasiparticle  current  (which 

corresponds  to  a  shifted  potential  due  to  the  presence 

of  a  time-averaged  charge  imbalance  in  the  strip) ,  which 

if  large  might  mask  the  desired  effect.  This  can  be 

minimized  by  making  the  upper  electrode  out  of  a  material 

with  a  hicher  T  than  that  of  Ihc  strip,  which  will  be 
c 

just  below  its  own  T  .  Secondly,  the  phase  difference  AO 
is  dependent  not  only  on  time,  but  on  the  position  along 
the  strip.  In  analogy  to  the  time  domain,  this  spatial 
dependence  will  consist  of  the  sum  of  a  linear  part  due 
to  the  dc  supercurrent  in  the  strip,  and  an  oscillatory 
part  due  to  the  waves.  If  the  tunnel  junction  probe 
is  larger  than  either  the  wavelength  of  the  charge 
imbalance  waves  or  the  coherence  length  f,  (T)  (assuming 


the  presence  of  a  near-critical  supercurrent  in  the 
strip) ,  then  the  size  of  the  desired  effect  will  be 
attentuated.  Fortunately,  a  tunnel  junction  -  1  um  wide 
or  smaller  (within  the  capabilities  of  standard  photoresist 
and  liftoff  techniques)  seems  to  satisfy  these  require¬ 
ments  for  reasonable  materials,  temperatures,  and  voltages. 

A  final  point  is  that  it  has  been  assumed  that  the 
tunnel  junction  is  acting  as  a  passive  probe  of  the  strip, 
although  it  will  be  injecting  some  current  into  the  strip. 
Ideally  this  current  will  be  much  smaller  than  the  currents 
otherwise  flowing  in  the  strip.  If  it  is  too  small, 
however,  then  the  desired  effects  would  likely  be  washed 
out  by  noise  in  the  tunnel  junction  characteristics. 

In  the  opposite  extreme  of  large  injected  currents,  one 
can  determine  the  effect  of  the  tunnel  junction  on  the 
strip  by  reversing  the  situation,  and  using  the  I-V 
characteristics  of  the  PSC  to  probe  the  disequilibrium 
caused  by  the  tunnel  junction. 

The  experimental  configuration  we  have  been  discussing, 
as  indicated  in  Fig.  5.1b,  is  topologically  equivalent 
to  the  configuration  in  Fig.  5.1c  used  to  investigate 

the  interactions  between  two  closely-spaced  microbridges. 

6  5  6  7 

These  experiments  '  and  others  on  arrays  of  micro- 
6  6 

bridges  have  shown  that  under  certain  circumstances 
the  Josephson  oscillations  of  closely-spaced  microbridges 
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will  synchronize,  i.e.  the  voltages  across  the  individual 
bridges  tend  to  lock  together.  By  analogy  with  the 
tunnel  junction  detection  scheme,  our  approach  leads  us 
to  believe  that  these  effects  may  be  due  in  part  to 
coupling  via  charge  imbalance  waves  emitted  by  the 
adjacent  bridges.  Artemenko,  Volkov,  and  Zaitsev'’7 
have  made  apparently  similar  predictions  in  discussing 
collective  excitations  in  large  arrays  of  proximity 
effect  bridqes,  and  Lindelof6 ‘  has  made  similar  suggestions. 
A  better  understanding  of  the  criteria  for  coupling 
between  microbridges  should  assist  in  the  develooment 
of  arrays  of  microbridges  with  potential  applications 
as  detectors  and  transmitters  of  radiation. 

In  experiments  investiaating  microbridge  coupling, 
the  critical  currents  of  the  two  microbridges  typically 
have  comparable  values,  so  that  each  has  a  significant 
effect  on  the  other,  and  the  treatment  becomes  somewhat 
more  complicated  than  the  tunnel  junction  probe  discussed 
earlier.  We  can  present  this  problem  in  a  tractable  form, 
however,  through  the  use  of  our  transmission  line  equiva¬ 
lent  circuit,  in  particular  the  formulas  for  the  finite 
lengths  of  transmission  line  developed  in  Appendix  III. 
Consider  the  simplified  problem  consisting  of  two  PSC's 
located  a  distance  21  apart  on  a  uniform  strip  (Fig.  5.2a), 
represented  in  our  transmission-line  picture  by  Fig.  5.2b. 


Fig.  5.2  Two  interacting  PSC's. 

a)  Superconducting  strip  with  two  weak  spots 
(indicated  by  dashed  lines)  which  act  to  nucleate 
PSC's. 

b)  Transmission-line  equivalent  of  (a) . 

c)  Lumped-impedance  equivalent  (using  the  pi- 
section  of  Eq.  (£3.13)  )  of  the  transmission 
line  in  (b) . 
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We  can  relate  the  currents  and  potentials  at  the  two 
Josephson  elements  by  replacing  the  length  of  the  strip 
between  the  two  with  its  lumped  impedance  equivalent, 
either  the  pi  equivalent  (Fig.  5.2c)  or  the  tee  equivalent 
(Fig.  5.5) .  The  impedances  here  are  the  same  as  those 
of  Fqs.  (A3. 13)  and  (A3. 14).  In  general,  the  currents 
at  one  of  the  Josephson  elements  will  modulate  the 
voltaqe  across  the  other,  and  vice  versa. 

Some  simple  consequences  of  these  pictures  can  be 
seen  even  without  working  through  the  details  of  the 
interaction.  Consider  a  string  of  identical  FSC's, 
each  adiacent  pair  separated  by  the  distance  21,  repre¬ 
sented  in  terms  of  the  tee  section  equivalent  [Fq.  (A3. 14)] 
in  Fig.  5.3a.  By  symmetry,  the  time-averaged  voltage 
V  across  each  Josephson  element  is  the  same,  as  is  the 

average  normal  current  I  at  each  element.  Therefore, 

n 

in  Fig.  5.3a,  the  average  current  passing  through  the 
impedances  7^  must  be  7.ero,  and  the  voltage  across  each 
PSC  can  be  expressed  as 

V  =  2ZJ  (w«0,2I)I  =  27.  (m=0)  T  tanh(K(0)C) 

1  n  o  n 

=  2ARI  tanh(8/A)  (5.1) 

using  Fq.  (A3.1^)f  where  R  is  the  resistance  ner  unit 

length  and  A  =  1/*T3C>  =  K  *  (0)  is  the  dc  decay  length. 

54 

Within  the  context  of  the  BBT  model,  Tinkham  considered 
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the  similar  problem  of  an  ideal  uniform  filament  of 
length  L,  with  n  identical  PSC's  distributed  evenly 
along  its  length.  This  corresponds  to  a  separation 
of  2 l  =  L/n  between  adjacent  PSC's,  with  a  distance 
i  =  L/2n  between  the  end  of  the  filament  and  the  initial 
(or  final)  PSC.  By  symmetry  the  nonequilibrium  potential 
<t>  =  0  half-way  between  each  pair  of  adjacent  PSC's, 
and  if  we  take  the  boundary  condition  at  the  end  to  be 
♦  -  0  as  well  (the  "banks"  are  maintained  in  equilibrium) , 
then  the  voltage  across  each  PSC  is  equal  to  that  in 
Eq.  (5.1),  and  the  total  voltage  across  the  entire 
filament  is 

VT(n)  =  2AnRIntanh [L/ ( 2nA) )  .  (5.2) 

This  is  the  same  as  Eq.  (1)  of  ref.  54,  derived  from 
direct  consideration  of  the  equations  for  the  SBT  model. 

Another  problem  that  can  be  treated  rather  simply 
within  the  transmission-line  picture  is  that  of  a  single 
PSC  in  a  filament  of  finite  length,  represented  in  Fig. 
5.3b  by  a  tee  section  equivalent  for  the  length  S,  on 
either  side  of  a  Josephson  element.  By  direct  calcula¬ 
tion  we  can  show  that 

ZJU)  +  [Zjm  ||7.£(*)  ]  =  Zotanh(Ke)  =  7.J(2*)f  (5.3) 

(+  and  ||  refer  to  the  series  and  parallel  combination 
of  the  impedances) ,  so  that  we  can  represent  the  problem 
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by  Fig.  5.3c,  which  shows  a  Josephson  element  shunted 
by  an  impedance 

Z  =  2Z'  (2C)  =  2Z  tanh(KC).  (5.4) 

s  X  o 

(Note  that  in  dc  V  =  Z  I  reduces  to  the  result  in 

s  n 

Eq.  (5.2)  for  n=l ,  as  it  should.)  We  can  then  proceed 

to  solve  the  problem  using  the  methods  of  Appendix  II 

and  Chapter  IV.  In  Fig.  5.4  we  have  plotted  the  resulting 

I-V  curve  and  its  derivative  (dl/dV  vs.  V)  for  the  same 

parameters  as  curve  b  of  Fig.  4.8  (linear  f(AP),  ujcto=2.2, 

t„/t  =11)  and  usina  21  =  A,  The  I-V  curve  is  not  par- 
t  o 

ticularly  striking,  but  its  derivative  shows  a  series  of 
sharp  features.  Their  physical  origin  can  be  seen  from 
the  following  discussion. 

In  a  PSC  in  a  long  filament,  charge  imbalance  waves 
propagate  out  from  the  center  and  decay  in  a  distance 
of  order  A.  If  the  filament  is  rot  so  long,  the  waves 
reach  the  end  before  having  decayed  completely  to  zero, 
and  partially  reflect  back  due  to  the  impedance  mismatch 
at  the  boundary.  They  can  then  interact  back  on  the 
Josephson  element  in  the  center.  In  general  they  will  be 
phase-shifted  relative  to  the  waves  then  being  generated, 
with  the  greatest  effect  occuring  when  they  are  either 

O 

in  phase  or  180  out  of  phase.  These  resonance  conditions 
correspond  to  a  finite  number  of  wavelengths  or  half¬ 
wavelengths  contained  in  the  length  2A,  and  lead  to 


«mMM> 


V(2eWTi) 

Fig.  5.4  Electrical  characteristics  of  a  PSC  in  a 

strip  of  length  21  =  A,  with  the  same  parameters 
as  for  curve  b  of  Fig.  4.8.  a)  V  vs.  I. 
b)  dl/dV  vs.  v. 
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variations  in  the  magnitude  of  the  shunt  impedance  in 
Eq.  (5.4),  which  in  turn  cause  the  structure  in  Fig.  5.4. 

Returnina  to  the  interaction  of  two  PSC ' s ,  consider 
the  currents  and  Dotentials  using  the  tee  section  equi¬ 
valent  in  Fig.  5.5.  Their  Fourier  transforms  are  related 


V1  =  (Zo  +  ZP  Inl  +  Z2(Inl  "  I„2) 


V2  =  (Zo  +  Zi):[n2  +  Z2(In2  “  W 


(5.5) 


Using  these  together  with  the  current-phase  relations 
for  the  two  Josephson  elements,  the  problem  can  be 
solved  numerically.  We  can,  however  gain  a  better 
understanding  of  the  phenomena  involved  by  assuming 
that  the  two  oscillators  are  synchronized  at  the  same 
fundamental  frequency  Q  =  2eV/F\  (although  they  may  be 
out  of  phase) ,  so  that  the  currents  and  potentials  are 
all  periodic  in  time.  We  can  then  follow  through  with 
a  Fourier  series  method  similar  to  that  used  in  Chapter 
IV  and  Appendix  II.  Usinq  the  linear-periodic  super¬ 
current-phase  relation  of  Fq.  (A2.3),  we  obtain  the 
Fourier  coefficients 


'pi'01  =  In2<°>  ’  V/'VZi> 


(5.6) 
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Fig.  5.5  Lumped-el ement  equivalent  for  two  interactino 

PSC's,  separated  by  a  distance  2£,  using  the 
tee  section  [Eq.  (A3 .14) ] . 
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where 

ZT  =  Zq  +  Zj  +  Z'2  =  Zq[1  +  coth  ( 2K£ )  ]  (5.8) 

and  the  constant  t  is  the  time  lag  between  the  phase- 
slip  events  at  the  two  Josephson  elements.  The  ex¬ 
pression  for  I  _(m0)  takes  the  same  form  but  with  I  . 

nl  cl 

and  I  2  interchanged  and  e  *m^o  replaced  with  e+^m  to. 

The  total  current  at  each  element  is  equal  to 

[as  in  Eq.  (4.33)] 

00 

+  W2 

m=-°° 

(5.9) 

00 

h  mYi  I»2I"OT  +  W2 

m=-co 

For  the  present  case,  since  the  PSC's  are  fed  by  the  same 

constant  current,  I^=I2=I ,  and  this  condition  will 

determine  the  value  of  t  .  If  such  a  value  of  t  exists, 

o  o 

then  a  solution  with  the  two  oscillations  synchronized 
("voltage  locking")  is  possible.  If  I  j=Ic2  the  two 
oscillations  will  always  be  synchronized ,  but  if  they 
differ  by  even  a  small  amount  then  this  voltage  locking 
may  be  restricted  to  certain  currents  and  separations. 

We  may  obtain  an  estimate  of  the  strength  of  this  voltage 
locking  by  considering  Eq.  (5.7)  as  the  sum  of  two  terms, 
and  examining  the  second  term  which  contains  the  parameter 
t  which  can  be  adjusted  to  achieve  coupling.  Taking  for 


tfranift.  'i.j.jL.. 


SiH  % £i-„. 
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simplicity  m=l  and  Ic^=Ic-,,  the  absolute  value  of  this 
term  reduces  to  the  expression 


2  V 


F(V,n  =  I  /  |  4  ti  sinh  1 2K  (f!)  K]  - 


c  o 


I  Z  (ft) 
c  o 


cosh [2K (ft) C  ; 


(5.10) 


which  is  plotted  in  Fig.  5.6  as  a  function  of  V  for  a 
separation  2S  =  0.5 A  and  the  parameters  u>  to  =  Mclg  =  20 
[see  Fq.  (4.7)1.  Where  this  function  is  large,  voltage 
locking  should  occur  over  a  wider  range  of  critical 
currents  I  ^  and  Ic->-  The  periodic  peaks  (after  the 
first  on<')  correspond  to  a  resonance  condition  similar 
to  that  discussed  earlier  for  the  single  PSC;  the  two 
junctions  create  an  effective  resonant  cavity  between 
them. 


This  type  of  multiple-peak  structure  in  the  strength 
of  voltage  locking  is  reminiscent  of  qualitatively 
similar  structure  seen  by  Jillie  et  al.*’"’  in  the  behavior 
of  their  pairs  of  In  microbridges.  Although  the  experi¬ 
mental  situation  was  somewhat  more  complicated  than  our 
idealized  problem,  and  the  quantities  being  compared 
are  actually  different,  our  preliminary  results  are 
encouraging,  and  further  work  may  be  expected  to  clarify 
the  picture. 


coupled  microbridqes . 


5 . 2  Unexplained  Features  of  Microbridges  in  Magnetic  Fields 

We  emphasized  in  Chapter  II  the  importance  of  the 
localized  weak  spot  in  extending  the  range  of  useful  data 
undistorted  by  heating,  essentially  because  of  the 
reduction  in  the  local  critical  current  relative  to  that 
of  the  uniform  strip.  As  we  noted  in  Section  2.5,  this 
decrease  in  I  was  strongly  dependent  on  the  magnitude 
and  the  direction  of  the  applied  magnetic  field  (within 
the  plane  of  the  film) .  For  H  perpendicular  to  the 
current  density  J,  the  relative  "weakness"  of  the  weak 
spot  corresponded  to  a  simple  decrease  in  the  local 
cross-sectional  area.  For  H  ||J,  however,  it  was  far 
weaker  than  this  would  indicate.  We  speculated  in 
Section  2.5  that  this  difference  may  have  been  due  to 
the  presence  of  a  field  component  perpendicular  to  the 
surface  at  the  weak  spot,  because  of  its  geometry  (see 
Fig.  2.2).  Other  orientations  would  presumably  cause 
effects  between  these  two  limits,  although  these  measure¬ 
ments  were  not  carried  out,  in  oart  because  the  changes  in 
orientation  required  rotating  the  sample  at  room  tempera¬ 
ture  (the  magnet  rotation  was  used  to  orient  the  field 
parallel  to  the  film) . 

The  deoression  of  I  for  H  II J  consisted  of  two 

c 

basic  effects.  First,  the  local  critical  temperature 

T  (H)  at  the  weak  spot,  as  determined  by  the  temperature 
c 
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at  which  1^  went  to  zero,  was  depressed  below  the  value 

for  the  uniform  strio  (as  indicated  by  the  other  PSC's). 

2 

This  "extra"  denression  of  T  (H)  increased  as  H  ,  or 

c 

what  is  equivalent,  the  critical  field  curves  for  both 
the  uniform  strio  and  the  weak  spot  fit  the  form 

HC||(T)  =  1.44  «cp(0)  (1-T/TJ  (5.11) 

but  for  different  values  of  (0)  ,  with  that  for  the 

weak  spot  being  typically  10%  less.  It  is  not  clear 

whether  this  behavior  can  be  reconciled  with  the  suggestion 

about  a  perpendicular  component  of  field,  inasmuch  as 

H  |  (T)  (rather  than  its  scruare)  should  be  linear  in 

cj_ 

temperature.  In  this  regard,  Tinkham  has  made  an  alter- 
6  9 

native  suggestion  that  the  properties  of  the  weak  spot 

might  be  related  to  an  enhanced  value  of  the  parallel 

magnetic  field  at  the  weak  spot,  due  to  the  partial 

Meissner  effect  in  the  thin  film. 

Secondly,  even  after  its  onset,  the  critical 

current  of  the  weak  spot  for  H  ||J  increased  substantially 

more  slowly  as  the  temperature  was  lowered  than  that  of 

the  uniform  film.  The  Ginzburg-Landau  expression  for  the 

critical  current  of  a  uniform  film  in  a  oarallel  field 
6  8 

can  be  written  as 

IC(T,H)  =  Ic(T,H«0)U  -  H2/Hc^ (T) 1 3/2 
“  (1+t2)^  [l-t2-h2(l+t2) ]3/2 


(5.12) 


163 


where  t  =  T/Tc>  h  =  H/H c 1 1 ( 0 ) #  and  in  obtaining  the  latter 
expression  the  usual  two-fluid  temperature  dependences  of 


H  ^  and  A  have  been  used.  When  re-exDressed  in  terms  of 

T  (H)  ,  I  varies  as  [T  (H)-T]3/2,  and 
c  c  c 


dl 


2/3 


dT 


T=Tc (H) 


«  (1+h2) 1/6 (1-h2)^ 


(5.13) 


2/3 

This  was  compared  with  experiment  by  clotting  I 


vs  . 


c 

T  near  T  (H)  ,  and  measuring  the  slope  of  the  straight 
line  drawn  through  the  points.  For  both  H  ||j  and  H  j_  J , 
the  data  seemed  to  fit  a  straight  line  reasonably, 
although  for  the  latter  case  for  large  fields  there 
was  sometimes  some  minor  irregularity  or  even  irrepro- 
ducibilitv  of  Ic  as  a  function  of  T,  perhaps  associated 
with  trapped  flux.  (For  comparison,  the  behavior  I  <* 

[T  (H)-T]  has  been  observed  for  thin  films  in  a  perpen¬ 
dicular  field.  "^)  in  Fig.  5.7,  these  measured  slopes 
for  Sn  bridge  #29B  are  plotted  as  a  function  of  H  for 
the  two  orientations,  together  with  Eq.  (5.13)  normalized 
to  the  zero-field  behavior.  Although  the  results  for 
both  orientations  differ  from  the  theoretical  curve, 
those  for  H  ||J  differ  far  more.  The  reduced  value  of 
the  slope  for  H  ||J  reflects  the  comment  made  above, 
that  the  rate  of  increase  of  Ic  with  temperature  below 

T  (H)  is  reduced  below  the  value  for  the  uniform  film, 
c 

This,  taken  together  with  the  localized  depression  of 


fiMMi 
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Fig.  5.7  Effect  of  maqnetic  field  H  on  I  for  Sn 

2/3  , 

bridge  #29B  with  a  weak  spot:  d(Ic  )/3T 
for  T  near  m  (H)  ,  vs.  H  for  the  parallel  and 
perpendicular  orientations.  The  theory 
curve  is  from  Eo.  (5.13). 
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T  (II) ,  assures  that  the  local  value  of  I  is  well  below 
c  c 

that  for  the  uniform  film,  particularly  in  large  fields. 

The  particular  mechanism  for  this  effect  remains  to  be 
elucidated . 

Another  class  of  features  which  we  sometimes  observed 

in  bridnes  with  magnetically  weakened  spots  is  illustrated 

in  Fig.  5.8.  As  the  temperature  was  lowered  to  significantly 

below  T  (II)  ,  a  low  resistance  "foot"  seemed  to  develop 

for  currents  just  above  the  critical  current.  (This  is  the 

reason  for  some  of  the  anomalously  small  values  of  R  in 

Fin.  2.8.)  The  part  of  the  curve  for  somewhat  higher 

voltages  seemed  not  to  change  very  much  over  the  same 

temperature  interval,  and  our  supposition  is  that  this 

upeor  region  represents  the  usual  PRC  behavior.  In  the 

sample  for  which  the  data  in  Fig.  5.8  was  taken,  this 

"foot"  disappeared  for  large  fields,  and  for  small  fields 

the  correspond  inn  temperature  regime  was  dominated  by 

heating  effects.  The  exact  form  of  this  anomalous  low- 

voltage  behavior,  which  was  most  prominent  in  those 

samples  with  the  "weakest"  weak  spots,  varied  from 

sample  to  sample,  so  that  it  may  reflect  the  particular 

geometry  of  the  weak  spot.  It  is  speculated  that  perhaps 

these  features  may  be  effectively  a  short  microbridge 

phenomenon.  Qualitatively  similar  features  were  observed 

by  Octavio  in  the  I -V  character i sties  o f  Rn  variable- 

7 1 

•h  v  Me-.  .  microbridges  in  zero  maunetic  iield. 
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Fig.  5.8  I-V  character!  st  ics  for  Sn  bridue  #29R 
for  H  =  190  Oe  (  ||  J)  ,  showinq  the  development  of  a 
low-voltage  foot .  The  temperatures  corresponding  to 
the  curves  are:  A  3.620  K ;  B  3.606;  C  3.592;  p  3.579; 
F  3.564;  F  3.543;  G  3.497;  H  3.450;  I  3.415;  J  3.3S0; 

K  3.345;  L  3.313;  M  3.285;  N  3.258;  n  3.236;  P  3.212; 

Q  3.182;  R  3. 163. 


( vry )  I 


We  have  also  nade  some  preliminary  measurements 
on  short  Sn  two-dimensional  microbridnes  in  parallel 
magnetic  fields.  These  samples  were  made  on  a  uniform 

O 

(1000  A  thick)  Sn  film  by  cutting  with  the  sharp  point  of 

a  diamond  knife,  lifting  the  knife,  moving  a  short 

distance  further,  and  then  continuino  the  cut.  Tin 

bridge  #35S,  for  which  we  present  some  I-V  curves  in 

Fig.  5.9,  was  approximately  2  pm  long  by  2  pm  wide. 

These  curves  correspond  to  a  large  range  of  fields 

for  a  constant  temperature  of  about  3.2  K. 

We  will  focus  on  the  set  of  structures,  circled  in 

Fig.  5.9,  which  may  represent  some  form  of  "gap  structure" 

7  2 

occuring  at  V  =  2A(T,H)/e.  For  this  sample  these 
features  were  absent  in  zero  field  (although  other 
samples  did  show  comparable  gap  structure  in  zero  field) , 
and  were  also  lacking  in  large  fields,  for  which  a  low- 
voltage  foot  structure  became  prominent.  In  the  inset 
of  Fig.  5.9  we  have  plotted  the  square  of  the  voltage 
at  which  these  features  occur,  versus  the  square  of  the 
magnetic  field.  Taking  Hc||(T)  =  545  Oe  and  fitting  to 
the  relation 


A2 (T,H)  +  H2  =  1 
A2(T)  HcJ(T) 

we  obtain  the  value  A(T)  =  420  neV. 


(5.14) 

This  compares 


favorably  with  the  estimate  of  402  peV  obtained  by 


Fig.  5.9  I -V  characteristics  of  short  Sn  bridge  #35S 

(for  T  =  3.2  K) ,  showing  the  development  of  possible 

gap  structure  (indicated  by  small  dashed  circles). 

The  values  of  the  magnetic  field  are:  A  600  Oe;  B  580 

C  560;  D  r>5  0;  E  540;  F  530;  G  520;  H  510;  I  500;  J  490 

K  480;  L  470;  M  460;  N  450;  O  440;  P  420;  Q  400;  R  380 

S  360;  T  330  Oe .  The  inset  shows  a  straight-line  fit 

for  a  plot  of  the  sauare  of  the  voltage  at  which 

2 

this  structure  apoears  vs.  H  [see  Fo .  (5.14)]. 
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using 

A (T)  =  1.74  A ( 0 ) (1-T/T  )*  (5.15) 

together  with  the  values  Tc  =  3.81  K  and  y  =  3.20  K. 

If  these  features  were  in  fact  due  to  gap  structure, 
a  study  of  their  behavior  over  a  wide  range  of  tempera¬ 
tures  might  help  provide  some  insight  into  the  origins 
of  this  kind  of  structure  in  zero  field. 

The  behavior  of  the  small  microbridges  which  we 
examined  tended  to  differ  substantially  from  sample  to 
sample,  a  problem  that  was  also  present  in  the  low- 
voltage  behavior  of  long  bridges  with  weak  spots.  This 
may  be  due  in  part  to  the  difficulty  of  fabricating 
reproducible  geometries  on  this  small  scale  with  the 
diamond-knife  technology.  A  change  to  some  sort  of 
microlithographic  fabrication  therefore  seems  desirable. 
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CHAPTER  SIX:  SUMMARY  AND  CONCLUSIONS 

The  seemingly  simple  problem  of  the  current-induced 
breakdown  of  superconductivity  in  a  ouasi-one-dimensional 
filament  has  turned  out  to  have  a  great  deal  of  structure 
and  richness,  and  its  investigation  has  provided  major 
insights  into  a  whole  range  of  related  phenomena.  As  the 
current  through  the  filament  is  increased,  the  voltage 
increases  in  a  series  of  steps  alternating  with  plateaus 
of  anproximately  constant  differential  resistance.  As  was 
demonstrated  by  Skocool ,  Beasley,  and  Tinkham,*  this 
behavior  corresponds  to  the  successive  onset  of  localized 
dissipative  units  known  as  phase-slip  centers.  In  the 
research  which  we  have  reported  on  here,  we  have  continued 
the  earlier  work  of  SBT,  extending  and  confirming  their 
model  of  PSC ' s  in  superconducting  filaments,  in  both  its 
experimental  and  its  theoretical  aspects.  After  reviewing 
briefly  the  BBT  model  and  discussing  some  of  the  questions 
which  this  earlier  work  left  unanswered,  we  summarize  the 
present  work,  emphasizing  the  progress  that  has  been  made 
in  the  understanding  of  PSC's.  Finally,  we  conclude  with 
a  brief  discussion  of  some  of  the  major  implications  of 
our  work  for  future  research  in  the  field. 

Within  the  context  of  the  SBT  model,  a  PSC  is  a 
localized  dynamic  nonequilibrium  state,  which  retains  its 
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ability  to  carry  a  time-averaged  supercurrent  while 
supporting  a  dissipative  normal  current.  It  was  assumed 
that  the  problem  could  be  separated  into  the  phase-slip 
process  per  se ,  restricted  to  a  core  region  on  the  scale 
of  the  coherence  length  £,  and  the  diffusion  of  non¬ 
equilibrium  quasiparticles,  occuring  on  the  larger  scale 
of  A,  the  diffusion  length  corresponding  to  an  appropriate 
quasiparticle  relaxation  time  (as  originally  proposed, 
the  normal-state  electron-phonon  time  !„) .  SBT  used  a 

Ti 

particular  form  for  the  oscillations  in  the  core,  but 
all  that  it  essential  is  that  the  time-averaged  super¬ 
current  I  is  depressed  at  the  core,  so  that  an  inter¬ 
conversion  between  supercurrent  and  normal  current 
occurs  over  the  scale  of  A.  The  time-averaged  voltage  V 
across  the  PSC  is  then  proportional  to  I  at  the  core, 
with  the  proportionality  constant  being  the  normal 

resistance  R  of  a  length  2A  of  the  filament.  If  I 
n  s 

at  the  core  is  a  constant  fraction  of  the  critical  current 
Ic,  then  the  differential  resistance  dV/dl  of  the  PSC 
has  this  constant  value  R  ,  in  agreement  with  the  observa¬ 
tion  of  a  constant  slope  in  the  plateau  region  of  the 
I -V  curves.  The  voltage  profile  over  the  scale  of  A 
consistent  with  this  picture  involves  a  sudden  jump  in 
the  condensate  potential  u"s  at  the  core,  and  an  expo¬ 
nential  decay  of  the  auasiparticle  potential  un  to  the 
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value  of  us  over  the  length  A,  as  was  confirmed  by 
19 

Dolan  and  Jackel  with  normal  and  superconducting 
probes  near  a  PSC. 

Notwithstanding  a  substantial  amount  of  experi¬ 
mental  evidence,  the  SBT  model  was  largely  unsupported 
by  theory  at  the  time  that  it  was  proposed.  The  nature 
of  the  diffusion  length  and  the  corresponding  relaxation 
time  were  somewhat  uncertain,  as  was  the  connection 
between  the  PSC  and  other  experiments  involving  non¬ 
equilibrium  potential  differences  in  superconductos , 
such  as  the  quasiparticle  injection  experiment  of  Clarke.*'* 
Another  question  was  whether  the  model  really  required 
the  strict  separation  of  length  scales,  inasmuch  as  some 
of  the  strongest  experimental  support  came  from  tin 
samples  in  a  range  where  A  ®  3f>.  SBT  left  open  the  issue 
of  time-dependent  processes  outside  the  core,  dealing  only 
with  time  averages  in  the  region.  They  were  also  unspecific 
with  regard  to  the  shape  of  the  voltage  step  associated 
with  the  onset  of  the  PSC,  which  varied  experimentally 
from  rounded  through  steep  to  hysteretic.  Although  the 

importance  of  dissipative  heatinq  in  producing  hysteresis 

1  46 

was  emphasized  by  Skocpol  et  al.  '  the  possibility 
remained  of  some  intrinsic  hysteresis  in  the  regime  very 
close  to  T  .  The  question  of  interactions  between 
multiple  PSCs  was  also  not  explicitly  pursued  in  the 


original  SBT  paper,  although  Tinkham  has  since  described 
the  dc  interactions  which  follow  from  the  model.'’4 
Finally,  it  was  unclear  whether  the  SBT  picture  of  PSC's 
in  a  superconducting  filament  could  be  extended  to 
include  the  presence  of  oa i r-breakinq  perturbations 
such  as  a  magnetic  field.  We  have  tried  to  respond  to 
many  of  these  questions  in  our  research,  and  to  clarify 
the  issues  for  those  problems  we  could  not  resolve. 

In  the  course  of  developing  the  theoretical  back¬ 
ground  for  the  phase-slip  center,  we  have  emphasized  the 
centrality  of  the  concept  of  charge  imbalance  (or  branch 
imbalance)  to  all  phenomena  involving  the  existence  of 
a  nonequilibrium  potential  difference  in  a  superconductor . 
Basically,  this  is  because  a  charge  imbalance  is  the 
net  excess  charge  density  in  the  quasiparticles,  and  to 
maintain  overall  neutrality  the  condensate  must  respond 
by  shifting  its  own  chemical  potential.  Charge  imbalance 
can  be  generated  not  only  by  quasiparticle  injection,  but 
also  by  the  net  influx  of  normal  current,  which  transports 
quasiparticle  charge.  By  working  within  the  context  of 

recent  generalized  two-fluid  theories  of  nonequilibrium 
24 

superconductors,  which  are  rigorously  valid  in  the 
limit  of  small  A/k_T,  we  have  obtained  equations  describing 

n 

the  dynamics  of  charge  imbalance  in  nonequilibrium 
superconductors.  One  consequence  of  this  approach  is 


the  prediction,  still  experimentally  unverified,  that  a 
homogeneous  charge  imbalance  will  decay  to  zero  expo¬ 
nentially  with  a  relaxat’.-  time  x^,.  rather  than  the 

"branch  imbalance  relaxation  time"  x~*  which  was  intro- 

Q* 

9 

duced  by  Tinkham  and  Clarke  for  the  dc  quasiparticle 
tunnel  injection  experiment.  ^hese  equations  are  exactly 
equivalent  to  those  describing  electrical  signals  propa¬ 
gating  down  a  classical  transmission  line,  where  the 
two  channels  of  the  line  represent  the  condensate  and 
quasiparticle  channels  for  the  flow  of  current,  and 
the  shunt  conductance  between  the  two  lines  represents 
conversion  processes  via  charge  imbalance  decay.  These 
equations  and  the  associated  transmission  line  form  a 
unified  framework  for  dealing  with  the  whole  gamut  of 
phenomena  involving  charge  imbalance.  The  charge  imbalance 
waves  which  propaqate  down  the  superconducting  filament 
have  a  decay  length  whose  low-frequency  limit  is  =  /Dx^* 

This  quasiparticle  diffusion  length  thus  characterizes 
the  exponential  decay  of  normal  current  on  the  supercon¬ 
ducting  side  of  a  current-carrying  normal -superconduct i ng 
boundary.  For  hiah  frequencies,  these  charge  imbalance 
waves  correspond  to  the  propagating  collective  mode 

observed  in  superconducting  fluctuations  by  Carlson 
23 

and  Goldman. 

We  have  also  used  this  transmission-line  equivalent 
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to  develop  a  new  model  for  a  phase-slip  center.  The 
Josephson  oscillations  connected  with  the  phase-slip 
process  in  the  core  act  to  generate  charge  imbalance 
waves  which  oropaqate  out  in  both  directions  alonq 
the  filament.  Continuinq  for  simplicity  the  separation 
of  length  scales  contained  in  the  SBT  model,  we  describe 
these  Joseohson  oscillations  in  terms  of  a  sunercur rent- 
phase  relation  for  an  ideal  Josephson  element  of  negligible 
extent,  located  on  the  condensate  channel  of  our  trans¬ 
mission  line.  The  effect  of  the  transmission  line  on 
the  core  is  equivalent  to  shunting  the  Josephson  element 
with  an  impedance  equal  to  twice  the  characteristic 
impedance  of  the  line.  The  resulting  time-averaged 
behavior  reduces  essentially  to  that  of  SBT,  except 
that  now  the  length  A  is  unambiguously  identified  with 

A r...  The  I-V  characteristics  calculated  within  the 

w 

model  show  a  low-voltage  onset  whose  shape  depends  on 
the  characteristic  times  in  the  problem,  and  a  high- 
voltage  asymptote  with  constant  differential  resistance 
proportional  to  A0* .  The  onset  becomes  hysteretic 
when  the  effective  ac  resistance  shunting  the  Josephson 
element  becomes  smaller  than  the  dc  resistance.  Pre¬ 
liminary  experimental  evidence  on  this  issue  is  incon¬ 
clusive,  but  suggests  qualitative  and  some  rough  quanti¬ 
tative  agreement  in  tin  microstrips  and  whiskers.  With 
pairs  of  closely-spaced  PSC's,  the  model  predicts  that 


178 


the  charge  imbalance  waves  emitted  by  each  junction  will 
modulate  the  voltage  driving  the  other,  leading  to  the 
possibility  of  "voltage  locking"  or  synchronization  of 
the  Josephson  oscillations.  Other  effects  similar  to 
resonance  are  predicted  in  situations  where  discontinuities 
in  the  wave  properties  cause  reflections  of  the  charge 
imbalance  waves  back  onto  the  Josephson  element  which 
had  earlier  emitted  them.  All  these  effects  can  be 
treated  very  conveniently  using  standard  formulas  for 
lumped-element  equivalents  of  finite  lengths  of  the 
transmission  line.  Again,  these  predictions  are  quali¬ 
tatively  similar  to  certain  effects  reported  in  the 
literature  for  closely-spaced  microbridges,  but  more 
work  is  required  to  provide  a  rigorous  test  of  the  model. 

Our  experimental  work  on  tin  microstrips  in  parallel 

magnetic  fields  has  provided  additional  evidence  that 

the  basic  concepts  contained  in  the  SBT  model  of  PSC's 

remain  valid  well  beyond  the  context  in  which  they 

were  originally  proposed.  The  qualitative  behavior 

of  the  I-V  curves  remains  the  same,  but  the  characteristic 

resistance  level  R  decreases  substantially  as  the  field 

n 

is  increased.  The  approximation  of  the  separation  of 
length  scales  is  then  even  poorer  than  in  zero  field, 
but  if  we  extract  a  quasiparticle  relaxation  time  from 
Rn  in  the  usual  way,  we  obtain  general  agreement  with  a 
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formula  derived  by  Schmid  and  Schon  for  the  "transverse 
mode"  relaxation  time  in  the  presence  of  pair-breaking 
(produced  by  the  magnetic  field) .  Since  this  transverse 
mode  time  is  effectively  a  generalized  charge  imbalance 
relaxation  time,  this  agreement  provides  support  for  our 
theoretical  analysis  of  PSC’s.  The  clearest  agreement 
with  theory  was  obtained  form  our  samples  containing 
deliberate  weak  spots  which  depressed  Ic  locally, 

isolating  a  single  PSC  and  decreasing  heat  dissipation 

2  -Jfc 

«  Ic  R  .  The  value  of  Rn  diverged  weakly  as  [Tc(H)-T] 

at  the  critical  temperature  of  the  uniform  strip,  rather 

than  at  that  of  the  weak  spot.  This  helps  confirm  the 

contention  that  the  dc  properties  of  PSC's  are  relatively 

independent  of  the  dynamics  in  the  core;  the  key  is  the 

quasiparticle  relaxation  over  the  larger  scale  of  A. 

This  divergence  was  not  apparent  in  much  of  the  earlier 

work  with  uniform  strips  in  zero  field,  in  part  because 

of  the  effects  of  heating  in  enhancing  the  apparent 

resistance  level  away  from  T  .  We  were  able  to  remove 

these  effects  for  moderate  amounts  of  heating  using  the 

local  heating  model  of  Skocpol  et  al .  '  ,  and  thereby 

expand  the  temperature  range  near  T  (H)  where  we  could 

accurately  infer  the  relaxation  time.  The  mechanism  for 

the  reduction  of  I  at  the  weak  spots  remains  uncertain, 

although  the  magnetic  field  clearly  plays  an  important 

part. 


1  «i  ■ 
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Our  investigations  into  suoerconduct  i  ng  phase- 
slip  centers  have  helped  to  clarify  the  relationships 
between  the  I-V  characteristics  and  the  nonequilibrium 
phenomena  associated  with  charge  imbalance  and  heatinq. 

At  the  same  time,  they  have  introduced  a  series  of  new 
questions,  which  provide  an  agenda  for  future  research 
in  the  field.  The  most  important  of  these  relate  to  the 
ac  properties  of  PSC's,  for  which  experimental  evidence 
for  the  theoretical  model  we  have  developed  is  largely 
lacking.  Careful  measurement  of  hysteresis  close  to  Tc 
in  a  variety  of  materials  may  provide  an  indirect  test, 
although  the  presence  of  heating  may  complicate  the  inter¬ 
pretation.  More  clearcut  evidence  might  be  obtained  by 
the  observation  of  resonance-related  effects  in  a  single 
PSC  and  vol tage- lock i ng  behavior  in  accordance  with  the 
model  in  closely-spaced  microbridges.  Also  of  significant 
interest  is  whether  our  model  can  be  extended  to  include 
the  effects  of  magnetic  fields  and  incident  microwave 
radiation  on  PSC's  and  microbridges.  In  conclusion, 
the  study  of  the  T-V  characteristics  of  superconducting 
filaments  has  helped  to  provide  insinht  into  a  wide 
range  of  nonequilibrium  phenomena  in  superconductors , 
and  will  likely  continue  to  do  so  in  the  future. 


APPENDIX  ONE:  ELECTRONIC  'rpMPrRATURE  REP ELATION  SYSTF‘1 

The  main  purpose  of  this  appendix  is  to  provide  a 
brief  guide  to  the  operation  of  our  home-made  electronic 
temperature  regulator,  in  case  there  should  be  a  use 
for  it  in  the  future.  In  my  experiments,  for  those  runs 
in  which  the  sample  was  located  inside  a  vacuum  can,  it 
was  used  to  provide  electronic  feedback  to  requlate  the 
temperature  of  the  copper  block  on  which  the  sample  was 
mounted.  It  usually  worked  ouite  well,  although  as  it 
turned  out  it  was  generally  sufficient  to  run  directly  in 
the  liquid  He  bath,  without  having  to  worry  about  whether 
the  vacuum  can  was  actually  sealed. 

The  circuit  functioned  basically  as  a  simple  thermostat 
system  with  electronic  feedback.  The  signal  from  a 
temperature  sensor  was  amplified,  compared  with  a  thermo¬ 
stat  setting,  and  a  correction  signal  was  applied  to 
adjust  the  power  delivered  to  a  heater.  A  schematic  of 
the  circuit  is  shown  in  Fig.  A1 . 1 .  In  more  detail,  the 
thermometer  signal  was  the  voltaae  across  a  Tryocal 
doped  germanium  resistance  thermometer  (or  earlier  an 
A1 len- Brad  ley  resistor)  with  4.2  K  resistance  of  order 
1  kohm  and  a  negative  temperature  coefficient,  biased  with 
a  constant  current  of  10  uA  or  less.  It  was  mounted  on 
the  same  copper  block  as  was  the  sample  and  also  the  heater, 
which  consisted  of  a  length  of  alloy  resistance  wire 
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Resistonce  Heoter 
Thermometer 


('vlOO  n)  wrapped  around  the  block  and  heat-sunk  with 
GF-7031  varnish.  The  thermometer  voltage  was  measured 
with  a  Keithley  model  160  digital  voltmeter,  and  the 
analog  output  of  the  DVM  was  used  to  amDlify  the  signal 
to  the  1  volt  range  (usually  by  a  factor  of  100) .  The 
difference  between  this  and  the  thermostat  voltage  was 
amplified  and  added  to  the  raw  heater  bias  voltage, 
which  was  then  used  to  supply  power  for  the  heater. 

Because  the  thermometer  had  a  negative  temperature 
coefficient,  the  circuit  was  set  up  so  that  an  increase 
in  the  thermometer  voltage  caused  a  proportional  increase 
in  the  heater  current. 

The  outline  of  the  operating  procedure  is  as  follows: 
First,  with  the  feedback  switch  open,  the  bath  was  pumped 
down  to  somewhat  below  the  operating  range  of  temperatures, 
and  the  raw  heater  bias  was  turned  up  to  yield  approximately 
the  desired  temperature.  The  temperature-setting  knob 
was  then  turned  to  null  the  amplified  correction  signal, 
and  the  feedback  was  then  switched  in.  The  feedback  would 
then  act  to  stabilize  the  temperature  at  the  given  value, 
even  while  the  bath  temperature  was  drifting  somewhat. 

The  stability  was  better  than  required  for  these  experi¬ 
ments,  permitting  a  drift  in  the  temperature  of  the 
copper  block  of  less  than  about  10  mK  per  hour.  Small 
adjustments  in  temperature  were  made  by  adjusting  the 
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thermometer  setting,  and  sometimes  also  the  heater 
bias  if  the  correction  voltage  got  too  large.  An  attempt 
to  make  a  large  change  of  temperature  in  this  way  some¬ 
times  led  to  unstable  oscillations.  A  tendency  towards 
oscillation  was  sometimes  a  problem  even  for  a  fixed 
setting,  particularly  for  the  larger  feedback  gains. 

This  may  have  been  due  to  a  time  lag  in  the  thermal 
response  of  the  cooper  block,  causing  the  heater  correction 
voltage  to  be  out  of  phase  with  the  temperature.  A 
systematic  study  of  this  was  not  made,  but  the  insertion 
of  a  capacitor  to  slow  down  the  electrical  response  of 
the  system  usually  got  rid  of  those  osc i 1 1  at  ions . 

Operated  in  the  way  described  above,  the  electronic 
temperature  regulation  system  provided  fairly  reliable 
constant  temperatures  in  the  range  from  about  1.5  K  to 


over  4.2  K . 


APPENDIX  TWO:  GENrRALIZED  RFSISTIVELY  SHUNTED  JUNCTION  MODEL 

We  treat  in  this  appendix  the  resistively  shunted 
junction  (RSJ)  model,  generalized  to  allow  for  a  non- 
sinusoidal  supercurrent-phase  relation. ®  In  particular, 
we  consider  an  ideal  Josephson  element  which  can  be 
described  by  the  supercurrent-phase  relation 

Is  =  Icf(A0)  (A2.1) 

where  I  is  the  critical  current  and  f  is  some  periodic 
function  of  the  phase  difference  A0  across  the  junction 
with  amplitude  one.  We  will  emphasize  the  classic 
sinusoidal  relation 

f ( A  6 )  =  sin  ( A0 )  (A2.2) 

and  the  linear-periodic  relation 

f  ( A0)  =  mod  (  A0/2tt)  (A2.3) 

(mod  or  modulus  is  the  remainder  after  the  integral  part 
is  subtracted  off),  but  unless  otherwise  noted,  the 
treatment  will  be  valid  for  a  general  f(A0).  As  indicated 
in  Fig.  A2.1,  this  junction  is  shunted  by  a  simple 
resistance  R.  This  parallel  combination  is  fed  with  a 
constant  current  bias  I,  and  the  voltage  drop  V  =  I^R 
will  be  a  function  of  time.  By  the  Josephson  relation 

V  =  A0  (A2.4) 
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Fig.  A2.1  The  generalized  resistively-shunted  junction 
(GRSJ)  model,  a)  Schematic,  b)  Two  examples 
of  periodic  functions  f ( A  © )  which  are  permitted 
in  this  model  for  the  supercurrent-phase 
relation  [Fes.  (A2.2)  and  (A2.3)]. 
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we  have 

1  •  1,  *  ■„  -  V'4"1  +  K5  4e  •  <fl2-5> 

This  differential  eauation  can  be  integrated  in  the  form 


/Ae  (t )  . 

d ( AO  ' )  =  / 

46(0)  >T  -  J° 


t  2el  R  2el  Rt 

-^C-  dt*  =  — jp —  .  (A2.6) 


For  the  dc  I-V  relation,  we  can  integrate  over  an 
entire  period,  and  using  the  fact  that 


_  .  2tt  2eV 


(A2.7) 


(f!  and  t  are  the  Josephson  anqular  frequency  and  period) 
J 

we  obtain 

2n 


V  =  2*1  R 


d  (AO) 


0  (f  -  f(A0)l 

"c 


(A2. 8) 


For  any  current-phase  relation,  this  starts  out  with 
V=0  at  Ic,  rises  with  vertical  slope  (taking  current  as 
the  horizontal  axis)  and  bends  over  to  meet  asymptotically 
the  voltage  biased  line 


V  =  R(  I  -  <f>I  ) 


(A2.9) 


where 


21 J  0 


<f>  =  /  f  ( A  6 ) d  ( A  9 )  .  (A2.10) 

Fq.  (A2.9)  can  also  be  written  in  the  form 


1  -  S  +  <f>Ic 


( A2 . 1 1 ) 


where  the  second  term,  the  "excess  current",  corresponds 


a  Fourier-series  method.  Although  this  leaves  the 
solution  in  the  form  of  an  infinite  sum,  it  is  more 
generally  applicable  to  cases  where  the  element  shunting 
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I/Ic 

Fig.  A2 . 2  V(I)  for  sinusoidal  f(A0).  The  dashed  lino 
is  the  high-voltage  asymptote.  The  inset 
shows  the  time  dependence  of  1^  for  I  =  1.1  1 
corresponding  to  the  small  circle  on  the 
V(l)  curve. 


V ( I )  for  linear-periodic  f  (  A 0 )  .  The  dashed 
line  is  the  high-voltage  asymptote.  The  inset 
shows  the  time  dependence  of  Ig  for  I  =  1.1  Ic 
corresponding  to  the  small  circle  on  the 
V ( I )  curve. 


Given  any  relation  f(A0),  Eqs .  (A2.18)  and  (A2.20) 

provide  a  formal  solution  to  the  problem,  which  can 
be  evaluated  on  a  computer  by  iteration  (although  a 
solution  may  not  always  exist  for  a  complex  Zg) .  The 
I-V  relation  can  be  expressed  in  the  form 


I 


10.' 


«  I  (0)  +  I  (0) 

n  s 

\K> 

=  ^  c  +  1  |f  (At’(Of)  )  +  f  (.'•  (0~)  )  1/2  (A?. 21) 

m=-"'  ' 

where  the  second  term  .Allows  toi  the  discontuity  ot 
f ( AO ( t ) )  at  t  =  0 . 

If  we  use  tlie  linear-periodic  current -phase  relation 


of  Eq.  (A2.3),  an  analytic  expression  for  the  Fourier 
coefficients  can  be  obtained.  Pubstitut  mu  Fq .  (A2.20) 

into  Fa.  (A2.18),  we  obtain 


which  yields  the  explicit  solution  for  c 


x-4-  /(l  + 

2nim  2  mm 


Z  imd) 
s 


(A2 . 23) 


For  the  resistive  shunt,  we  simply  set  Z  (mO)  R . 
Provided  that  a  sufficient  number  of  harmonics  are  taken, 
this  reproduces  the  results  of  Fas.  (A2.ll)  and  (A.’.lr>). 
The  more  qeneral  form  with  the  shunt  impedance  is  used 


in  Chapter  TV. 


AFPFNDI X  THREE :  FINITE  LENGTH  AND  ARTIFICIAL 

TRANSMISSION  LINES 

In  this  appendix  we  treat  the  transmission  line  of 
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Chapter  III  for  cases  involvinq  finite  lengths,  where 
the  possibility  of  reflected  waves  cannot  be  ignored. 

We  then  show  how  this  information  can  be  used  to  design 
an  artif.eial  transmission  line,  which  might  serve  in 
analog  computations  of  some  of  the  problems  in  Chapters 
IV  and  V  where  numerical  computations  become  difficult. 

Consider  the  section  of  transmission  line  in  Fig. 

A3. la  of  length  2 C ,  between  x  =  -C  and  x  =  +S.  As  in 
Chapter  IIT,  its  continuous  electrical  elements  can  be 
characterized  by  a  resistance  R  per  unit  length  in  the 
upper  channel,  an  inductance  I.  per  unit  length  in  the 
lower  channel  ,  and  a  conductance  C.  and  capacitance  C 
per  unit  length  between  the  two  channels.  From  Kirchoff's 
laws,  the  equations  governing  the  evolution  of  the 
currents  and  potentials  are 

=  -RI  +  LI  (A3.  1) 

3x  n  s 

31 

n  =  -g<J>  -  C*  ,  (A3.  2) 

3x 

where  as  before  I  and  1^  are  the  currents  in  the  upper 
and  lower  channels  respectively  and  <t>  is  the  potential 
difference  between  the  two  channels.  We  will  be  dealing 
exclusively  here  with  the  case  of  constant  total  current, 
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Fig.  A3.1  Equivalents  for  finite  length  transmission  1 

a)  Schematic  of  transmission-line  section  of  length  2?. 

b)  Pi-section  equivalent  of  (a) . 

c)  'T’ee-section  equivalent  of  (a)  . 


r 


195 


so  that  I 

s 

=  -1 

n 

and 

Fa .  ( A3 . 1 )  can  be 

rewritten 

2*  = 
3x 

-R I 

n 

LI 

n 

• 

(A3. 3) 

Under  this 

cond  i 

t  ion 

we  may  also  place 

the  inductance 

in  Fig.  A3 

.la  in 

the 

upner  channel  and 

leave  the  lower 

channel  as 

a  "common 

ground";  Fqs.  (A3 

.2)  and  (A3. 3) 

are  still 

va  lid. 

(I 

n  the  more  general 

case  where  1^0, 

this  equivalence  and  the  following  treatment  do  not 
apply;  however,  a  modified  approach  can  be  developed 
which  will  allow  for  time-varying  total  currents.) 

From  Fgs.  (A3. 2)  and  (A3. 3) ,  both  $  and  I  satisfy 


the  equation 
3  J  f 


3x ' 


=  LCf  +  ( RC+GL)  f  +  RGf 


with  dispersion  relation 

ik  K  =  * (R+iwL) (fi+iwC) 


(A3. 4) 


(A3. 5) 


Consider  I  for  a  particular  frequency  component: 


T  .  lkx  ,  _  -lkx.  iu't 

I  ( x  ,  t )  =  ( Ae  +  Be  )e 

n 


(A3 .6) 


Here  A  and  R  are  arbitrary  constants,  and  the  first 
term  corresponds  to  a  wave  propagating  to  the  left, 
the  second  term  to  the  right.  Since  we  will  be  dealing 
only  with  a  single  freouency  component,  we  will  for 
simplicity  suppress  the  time  dependence: 

I  (x)  =  Aelkx  +  Be-1 kx  .  (A3. 7) 

n 


1 


19  6 


From  Ecr.  (A3. 3)  we  have 

1^-  =  -  (R+iuiL)  I  (A3. 8) 

3x  n 

or 

*  _  E+ia)L  ,»_ikx  „_-ikx> 

„  ,  .  ikx  ,  ^  -ikx.  ,  „ ,  n. 

=  Zq  (-Ae  +  Be  )  ,  (A3. 9) 

where 

Z  (w)  =  /(E+iuiL)  /(C  +  iwC)  (A3. 10) 

o 

is  the  characteristic  impedance  of  the  line.  If  there 
is  only  a  wave  propagating  to  the  right  then  i>  =  z0Inf 
while  a  wave  propagating  to  the  left  will  give  4>  =  -ZQIn. 
When  both  are  present  the  more  complete  relations  of 
Fqs .  (A3. 7)  and  (A3. 9)  are  necessary. 

The  transmission  line  section  of  Fiq.  A3. la  can 
be  treated  as  a  four-terminal  linear  passive  device, 
or  if  the  inductance  is  in  the  upper  line,  an  effective 
three-terminal  device.  We  can  therefore  look  for  a 


simple  equivalent  circuit  with  lumped  impedances  as  in 
Figs.  A3. lb  and  c,  called  for  obvious  reasons  the  "pi" 
abd  "tee"  equivalents.  (The  two  impedances  labeled 


<Mi) 
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Considerinq  the  pi  eouivalent,  we  introduce  I  °  as  the 


effective  current  passing  through  the  impedance  7.^. 
By  Kirchoff's  laws,  we  have  three  exoressions  for  I 


Xn  =  In(-,)  ' 

In°  =  In(t)  + 

in°  =  m-n  -  m)  ]/z2  . 


(A3. 12) 


We  can  solve  these  three  equations  for  the  three  unknowns 

Z .  ,  Z_,  and  I  ° : 

12  n 


Z,  =  Z  coth(KC) 

1  o 

Z~  =  27,  sinh  ( KC )  cosh  ( K  { )  =  7  sinh(2KC) 

2  o  o 


(A3. 13) 


I  °  =  (A+B) /cosh (KC ) 
n 

The  impedances  are,  as  thev  should  be,  independent  of 
the  amplitudes  of  the  waves.  For  the  tee  equivalent, 
a  similar  derivation  is  possible,  but  it  is  sufficient 
to  note  that  the  two  are  connected  by  the  standard 
"wye-delta"  (Y-A)  transformation: 


7»  Z 

Jzjrij  *  zotanh(Ktl 

7i 

Iz^zT  "  Vsinh(2KE> 

1  z 


(A3. 14) 


In  the  limit  where  the  length  of  the  transmission 
line  becomes  long,  the  two  ends  become  decoupled  with 
input  and  output  impedances  equal  to  Z  : 

Z  ^ ,  7,  j '  -*■  ZQ  ;  2  2  ■*  00  ;  Z?'  -*•  0  as  C  ■+  «*.  (A3. 15) 


1.  9R 


(The  hyperbolic  functions  in  the  impedances  general ly 

have  oscillatory  parts  as  well,  hut  in  the  limit  of 

large  8  the  leadincj  exponential  term  will  dominate  the 

sinusoidal  oscillations.)  In  the  other  limit  of  small 

8  ,  we  have 

Z  , 

Z1  "  K?  =  1  (G+iu,C) 

Z~  •*  2Z  K8  =  (R+iuL)  ( 2  e ) 

2  o 

,8-0  (A3. 16) 

Z!  ■*  7.  K 8  =  ( R+iwL)  8 
1  o 

Z  . 

Z2  *  2Kl  =  1  (G  +  i“'C)  (21)  l"1 

As  one  would  expect,  this  reduces  the  section  to  the 
basic  elements  of  the  transmission  line:  for  the  small 
length  28,  a  series  impedance  (R+iwL) 28  and  a  shunt 
admittance  (G+iwC)28. 

The  transmission  line  can  of  course  be  constructed 
of  these  infinitesimal  elements,  but  it  can  also  be 
"constructed"  of  the  finite  length  equivalents  connected 
in  series.  A  semi-infinite  line  has  impedance  Zq,  since 
there  are  no  reflections.  If  we  attach  a  pi  section 
equivalent  to  this,  cS  in  Fig.  A3. 2a,  we  obtain  an 
impedance 

Z  =  Z1  ||  IZ2  +  (Zx  ||Zo)  1  =  ZQ  (A3. 17) 

as  can  be  shown  by  direct  calculation  (the  ||  and  + 
represent  parallel  and  series  connection  respectively). 


Fig.  A3. 2  Reconstruction  of  transmission  line  from 
pi-section  eauivalents.  a)  Insertion  of 
Di-section  does  not  alter  the  characteristic 
impedance  of  the  line.  b)  Representation 
of  the  transmission  line  as  a  series  of 
pi-section  equivalents,  c)  An  example  of 
a  reconstruction  via  approximate  pi-section 
equivalents,  as  discussed  in  the  text. 
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We  can  therefore  attach  an  additional  pi  section,  or  as 
many  as  we  want,  and  still  have  an  impedance  Zp . 

An  infinite  series  of  these  (Fig. A3. 2b)  is  therefore  an 
equivalent  for  the  entire  transmission  line. 

Unfortunately,  for  finite  leneths  these  lumped 

impedances  cannot  easily  be  constructed  out  of  the 

usual  resistors,  capacitors,  and  inductors.  However, 

we  can  use  the  small  length  approx imat  ions  of  F.q.  (A3. 16), 

and  estimate  how  great  the  errors  will  be.  Substituting 

these  back  into  the  definitions  of  the  impedances ,  we 

can  invert  Fqs.  (A3. 13)  and  (A3. 14)  to  obtain  an  effective 

propagation  constant  K'  and  characteristic  impedance  2d 

in  terms  of  the  ideal  parameters  K  and  Z  : 

o 

K’  =  ^  sinh"1(KC)  '  K 1 1  +  -...1 

J  (A3. 18) 

Z  ’  =  Z  /,'\  +  (KCT7  '  Z  [1  -  ~-^-’  +  ...l 
o  o  o  2 

In  the  non-ideal  line  thus  obtained,  partial  reflections 
may  occur  at  the  connections  between  elements,  and  this 
causes  the  deviation  from  ideality  in  the  line  parameters. 
In  the  limit  where  (  is  much  smaller  than  both  the  decay 
length  and  the  wavelength,  Fq .  (A3. 18)  indicates  that 
these  deviations  may  be  neglected. 

To  take  a  particular  example,  let  us  measure  length 
in  terms  of  the  number  of  repeated  units,  and  let  each 
unit  consist  of  the  following  elements,  as  shown  in  Fig. 


A3. 2c: 
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R  =  10  Q  ;  G  =  1  mq_1  (G-1  =  1  k«) 

L  =  100  mH  ;  C  =  0.1  uF 
The  characteristic  times  of  the  line  are  then 


(A3. 19) 


L  ,  n- 5 
-  =  10  sec 


C  i  n-4 
rz  =  10  sec 


(A3. 20) 


and  the  low-frequency  decay  length  and  high-frequency 
are 

(0)  =  1//RG  =  10  units 


velocity  are 
,  -1 


(A3. 21) 


v  =  1//LC  =  3  *  10  units/sec 


-1 


From  Eq.  (A3. 18),  since  we  have  l  =  (0)/20,  this 

analog  transmission  line  should  be  suitable  for  wavelengths 
down  to  approximately  10  units,  with  errors  of  not  more 
than  a  few  percent.  It  should  be  possible  to  construct 
such  a  line,  couple  it  to  one  or  more  analog  Josephson 
junctions,  and  use  the  resulting  system  to  simulate 
resonance  or  synchronization  effects  in  phase-slip 
centers  (see  Section  5.1). 
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